Here ben
I have a tutorial on how to do it
http://stevemosher.wordpress.com/step-8-the-sample-package/
On Sun, Dec 18, 2011 at 11:49 PM, Petr PIKAL wrote:
> Hi
>
> If I remember correctly I do
>
> start clear R -vanilla session
> copy my functions to it.
> run package.skeleton("some.name")
Dear R-users,
I would like to determine the probability of event at specific time using
cox model fit. On the development sample data I am able to get the
probability of a event at time point(t).
I need probability score of a event at specific time, using scoring scoring
dataset which will have o
Dear R users,
I have a very simple query.
I am using bootstrap (library(boot)). I want to access the std. error in
the results and further generate confidence interval(CI) for n no. of
samples which will give me n values for std. error and CI.
*
*
b1
ORDINARY NONPARAMETRIC BOOTSTRAP
Call:
boot(
Dear R users,
I have a problem. I would like to solve the following:
I have
pL = 1/(1+e^(-b0+b1))
pM = 1/(1+e^(-b0))
pH = 1/(1+e^(-b0-b1))
My target function is
TF= mean(pL,pM,pH) which must equal 0.5%
My non-linear constraint is
nl.Const = 1-(pM/pH), which must equal 20%, and would like t
Hi
After good advices i post this problem here again. Now the attached sambledata
is in csv form and it is a part of my real data. I have tried different
na.action operations but the result stays the same...
My code
> library(mlogit)
> library(foreign)
> z<-odbcConnectExcel("D:\\SAMBLEDATA.csv
Hi there,
This seems like it should be simple. I have a data frame of climate data
sampled every 10 min. I want to average the entire data frame into 30
min values (i.e., one value for each half hour). Functions like
running.mean give me a moving average but I want to reduce the size of
the
On 19/12/2011 09:12, Vikram Bahure wrote:
Dear R users,
I have a very simple query.
I am using bootstrap (library(boot)). I want to access the std. error in
the results and further generate confidence interval(CI) for n no. of
samples which will give me n values for std. error and CI.
*
*
b1
O
Hi,
19.12.2011 11:12, Ville Iiskola wrote:
Error in if (abs(x - oldx)< ftol) { : missing value where TRUE/FALSE
needed
The reason for this error is in the row 563. There the choice has
value 1 and Ie has missing value. If the choice has value 0 and Ie
has missing values, then there is no erro
Dear all
I tried to find some packages (or programs) for image analysis and
especially fractal dimension image analysis but so far I had not success.
It shall be used for particle surface layer analysis from TEM images.
Any suggestions?
Best regards
Petr
_
Michael Griffiths wrote
>
> Dear R users,
>
> I have a problem. I would like to solve the following:
> I have
>
> pL = 1/(1+e^(-b0+b1))
> pM = 1/(1+e^(-b0))
> pH = 1/(1+e^(-b0-b1))
>
> My target function is
>
> TF= mean(pL,pM,pH) which must equal 0.5%
>
> My non-linear constraint is
>
> nl.
Hi Keith,
You do need to reorganize your data. cor() will work on any number of variables
as long as they are columns in a matrix or data frame.
There are a lot of ways to reorganize data, of various power and complexity.
Here's one simple way:
> library(ecodist)
> WW_Sample_table <- with(WW_Sam
Hi,
I've always used FRAGSTATS, but it looks like the SDMTools package will
do it within R.
Sarah
On Mon, Dec 19, 2011 at 4:56 AM, Petr PIKAL wrote:
> Dear all
>
> I tried to find some packages (or programs) for image analysis and
> especially fractal dimension image analysis but so far I had n
Em 18-12-2011 18:55, Rolf Turner escreveu:
> On 19/12/11 04:29, Uwe Ligges wrote:
>>
>>
>> On 18.12.2011 12:58, Ronaldo Reis Júnior wrote:
>
>
>>> Why the write.xls dont find the object a inside a function?
>>
>> Because at least that part of the function is poorly written.
>
> Surely this should
Does this work for you:
> x <- read.table(text = " date time Voltage LwTempDownelling LwDownwelling
> LwDownwelling_min LwDownwelling_max LwTempUpwelling
+ 1 2011-11-01 00:00:00 2.73244717.30 30.0
14.0 39.5 17.83
+ 2 2011-11-01 00:10:00 2.731534
Hi,
Thank you David, I resolved the problem with ?text.
Komine
--
View this message in context:
http://r.789695.n4.nabble.com/Label-in-oblique-orientation-tp4211944p4213465.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-pro
Hi all,
I know this a general question, not specific for any R package, even so I
hope someone may give me his/her opinion on this.
I have a set of 20 candidate models in a binomial GLM. The global model has
52 estimable parameters and sample size is made of about 1500 observations.
The global mo
Hello
Is there any package for rendering or raytracing?
--
View this message in context:
http://r.789695.n4.nabble.com/rendering-or-raytracing-tp4213944p4213944.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mail
Hi
How to perform hausman test when a model is created using lmer ?
Also please any one give link to mail it to r-sig-mixed-models mailing
list. I'm not able to find it
Thanks in advance
Arun
--
View this message in context:
http://r.789695.n4.nabble.com/Hausman-test-for-lmer-model-
Hi ,
I have issue obtaining p.value in urersTest(x, type = c("DF-GLS",
"P-test"),lag.max = 4,doplot = FALSE)@test$p.value
urersTest(...)@test indicates correct results but upon obtaining p.value as
this the result is NULL. Which seems not to be correct .
Any help on this ?
Mike
--
View this m
Jana Makedonska wrote on 12/14/2011 07:19:23 PM:
> Dear R-users,
>
>
> I am a relatively new R user and I have a simple question. Is there a
> command attached to a R function, which performs matrix pooling?
>
>
> Thanks in advance,
> Jana
>
> --
>
> --
> "The world is full of mysteries. L
Sean Baumgarten wrote on 12/14/2011 06:38:08 PM:
> Hello,
>
> I have a data frame with hourly or sub-hourly weather records that span
> several years, and from that data frame I'm trying to select only the
> records taken closest to noon for each day. Here's what I've done so
far:
>
> #Add a co
On 19/12/2011 7:54 AM, skan wrote:
Hello
Is there any package for rendering or raytracing?
The rgl package uses OpenGL for rendering solid objects, but it doesn't
do raytracing. A nice addition would be to output a scene in a format
that some external ray-tracer (e.g. POV-ray) could render.
Hi, I'm sure this is simple, but I haven't been able to find this in TFM,
say I have some data in R like this (pasted here:
http://pastebin.com/raw.php?i=sjS9Zkup):
> head(df)
gender age smokes diseaseY
1 female 65 ever control 0.18
2 female 77 never control 0.12
3 male 40
Hmm, almost but my real data frame is much larger (128 columns).
The problem is with this line
newData <- do.call(rbind, tapply(x[4:8], x$interval, colMeans))
when I do
tapply(x[4:8], x$interval, colMeans)
all works but when I do
tapply(x[4:7], x$interval, colMeans)
I get
Error in tapply(
Almost except
tapply(x[4:8], x$interval, colMeans)
works but with a larger data frame I have problems, even
tapply(x[4:7], x$interval, colMeans)
gets the "arguments must have the same length" error. But they do! Any
ideas?
Thanks for your help Jim
On 12/19/2011 2:00 PM, jim holtman wrote:
Do
On Dec 19, 2011, at 9:09 AM, Brent Pedersen wrote:
Hi, I'm sure this is simple, but I haven't been able to find this in
TFM,
say I have some data in R like this (pasted here:
http://pastebin.com/raw.php?i=sjS9Zkup):
One of the reason this is not in TFM is that these are questions that
sho
On Dec 19, 2011, at 3:34 AM, aajit75 wrote:
Dear R-users,
I would like to determine the probability of event at specific time
using
cox model fit. On the development sample data I am able to get the
probability of a event at time point(t).
I need probability score of a event at specific tim
On Dec 19, 2011, at 15:09 , Brent Pedersen wrote:
> Hi, I'm sure this is simple, but I haven't been able to find this in TFM,
It's not _that_ simple. You likely need TFtextbook rather than TFM. Most (but
not all) will go into at least some detail of coding categorical variables
using dummy var
On 19.12.2011 13:30, Ronaldo Reis Júnior wrote:
Em 18-12-2011 18:55, Rolf Turner escreveu:
On 19/12/11 04:29, Uwe Ligges wrote:
On 18.12.2011 12:58, Ronaldo Reis Júnior wrote:
Why the write.xls dont find the object a inside a function?
Because at least that part of the function is poo
n.via...@libero.itg wrote on 12/16/2011 05:02:24 AM:
> Dear all,
> I'm having problems with the tick of my graph. I'mpcombining lines and
> barplot.
> For my I'm using the function axis combined with the function pretty to
have
> more efficient tick, but all my tick (for example, 300 as my max
Hi,
I'm a newbie to the world of HMMs and HMMs in R. I've had a look at
the hmm package and the RHmm package but I couldn't see anything
straightforward on how a labelled sequential dataset with observed
values and underlying states might be used to construct and train a
HMM based on that data and
Claus,
Given that you have observed the hidden states they are not really hidden
anymore ... so maybe an HMM is not what you are looking for.
Assuming that you had such super observation powers using something like:
by(obs,hid,table)
will get you the state dependent counts of the categories of
Suppose I have the following:
x1<-as.vector(rnorm(10))
x2<-as.vector(rnorm(10))
x3<-as.vector(rnorm(10))
x4<-as.vector(rnorm(10))
x5<-as.vector(rnorm(10))
x6<-as.vector(rnorm(10))
x7<-as.vector(rnorm(10))
x8<-as.vector(rnorm(10))
x9<-as.vector(rnorm(10))
x10<-as.vector(rnorm(10))
I would like the
Hi
Then OK, I mean raytracing directly with some R package.
Or maybe some optics package
I've seen there exist an optics package called planar but it's only for
reflection and transmission at planar interfaces.
--
View this message in context:
http://r.789695.n4.nabble.com/rendering-or-raytracin
I'd like to run an R code automatically on Windows 7, say once every 2
hours. Could anyone tell me how to manage it?
Thanks in advance!
Best regards,
Ken
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.et
Not sure what you mean by mean of each position, so here are two
possibilities: rowMeans and colMeans:
> x1<-as.vector(rnorm(10))
> x2<-as.vector(rnorm(10))
> x3<-as.vector(rnorm(10))
> x4<-as.vector(rnorm(10))
> x5<-as.vector(rnorm(10))
> x6<-as.vector(rnorm(10))
> x7<-as.vector(rnorm(10))
> x8<-
That is what the task scheduler is for. Just create a ".bat" file
that invokes R and schedule it every two hours. On Win XP it is
Control Panel/Scheduled Tasks.
On Mon, Dec 19, 2011 at 11:10 AM, Kenneth Zhang wrote:
> I'd like to run an R code automatically on Windows 7, say once every 2
> hour
On Dec 19, 2011, at 11:00 AM, Pablo wrote:
Suppose I have the following:
x1<-as.vector(rnorm(10))
x2<-as.vector(rnorm(10))
x3<-as.vector(rnorm(10))
x4<-as.vector(rnorm(10))
x5<-as.vector(rnorm(10))
x6<-as.vector(rnorm(10))
x7<-as.vector(rnorm(10))
x8<-as.vector(rnorm(10))
x9<-as.vector(rnorm(1
Dear All
I have a constrained optimisation problem, I want to maximise the following
function
t(weights) %*% CovarianceMatrix %*% weights
for the weights,
subject to constraints on each element within the weights & the weights
vector summing to 1.
i.e.
weights = (x1, x2, x3), where x1 is w
Dear all,
I am not sure if this mail is for R-help or should be sent to R-devel
instead, and therefore post to both.
While using nlrob from package 'robustbase', I ran into the following
problem:
For psi-functions that can become zero (e.g. psi.bisquare), weights in
the internal call to nl
Dear Petr,
Well, It is not R but works well:
Fractal dimension and lacunarity a plugin for ImageJ.
Best regards,
Andrés AM
*
*
2011/12/19 Sarah Goslee
> Hi,
>
> I've always used FRAGSTATS, but it looks like the SDMTools package will
> do it within R.
>
> Sarah
>
> On Mon, Dec 19, 2011 at 4:5
Greetings.
I've written a convenience function for multicollinearity diagnosis.
I'd like to report to the user the formula that is used in a
regression. I get output like this:
> mcDiagnose(m1)
[1] "The following auxiliary models are being estimated and returned in a list:"
[1] "`x1` ~ ."
formul
Does applying formula() to the output of lm() (or terms())
do what you want?
> d <- data.frame(y=1:10, x1=log(1:10), x2=sqrt(1:10), x3=1/(1:10))
> formula(lm(formula("y~."), data=d)) # dot expanded
y ~ x1 + x2 + x3
> lm(formula("y~."), data=d)$call[[2]] # your original
formula("y~.")
>
Hi,
package planar is concerned with the full electromagnetic problem at
planar interfaces, which is not very useful for raytracing. The cda
package includes a small demo interfacing either rgl or povray (via
system call) to visualize 3D clusters of metallic particles. A more
general interface to
Komine wrote:
Hi everybody,
I have a matrix with 3 columns (Date, MeanArea and SdArea). I want to
draw a figure showing the variable MeanArea in terms of the Date. But
instead to use the variable SdArea as bar error, I want to use “polygon
error”. I use this code but the output does not seem goo
TOPIC
My question regards the philosophy behind how R implements corrections to
chi-square statistical tests. At least in recent versions (I'm using 2.13.1
(2011-07-08) on OSX 10.6.8.), the chisq.test function applies the Yates
continuity correction for 2 by 2 contingency tables. But when used a
On 12/13/2011 02:15 AM, jalfaro wrote:
Hi Jim,
Thanks so much for your help. I have read several of your responses on this
mailing list and they have helped me out quite a bit as I have gotten more
and more used to R.
I am still a little confused here by your response.
I think you understood my
try the lsei function from the limSolve package.
On Mon, Dec 19, 2011 at 2:32 PM, Russell2 wrote:
> Dear All
>
> I have a constrained optimisation problem, I want to maximise the following
> function
>
> t(weights) %*% CovarianceMatrix %*% weights
>
> for the weights,
>
> subject to constraints
>> I tried to be clever by trying get(paste(paste("x", 1:6, sep=""),
>> collapse="+")) but it didn't work.
Maybe you meant something like this:
sapply(paste("x", 1:6, sep=""), function(x)mean(get(x)))
On Mon, Dec 19, 2011 at 2:30 PM, David Winsemius wrote:
>
> On Dec 19, 2011, at 11:00 AM, Pab
I don't know if this code match what you want. You should use tcltk library.
code:
library(tcltk)
z <- function () { cat("Hello you!\n"); .id <<- tcl("after", 1000, z)}
z()
If you want to cancel you can use:
.id <<- tcl("after", 1000, z)
tcl("after", "info", .id) # To get info about this schedul
Dear all,
Does any one know how to calculate the risk (death rate) at some time point
(say 5 yrs) based on Kaplan-Meier curve in coxph model in R? I know how to fit
the survival data with coxph. But need to calculate risk at some time point.
Many thanks.
Bill
_
On Dec 19, 2011, at 10:05 PM, Bill Hyman wrote:
Dear all,
Does any one know how to calculate the risk (death rate) at some
time point (say 5 yrs) based on Kaplan-Meier curve in coxph model in
R? I know how to fit the survival data with coxph. But need to
calculate risk at some time point
Hello, all.
I have come across a problem. Previously, when I recoded my DV, added the
variable to my dataset and ran the multinom, I was just fine. But I am doing
it again and I am getting this error.
R> Poland$trust <- trust
Error in `$<-.data.frame`(`*tmp*`, "trust", value = c(NA, NA, NA, N
Dear R users,
I am getting following error while using boot.ci. I have int.inc function
with 2 values. I am generating CI for the sample estimate.
*> med <- function(x,i) median(x[i])*
*> b1 <- boot(int.inc,med,2)*
*> ci.out <- boot.ci(b1,conf = c(0.95),type = c("bca"))*
*Error: cannot
Hi,
I'm very new to working with sparse matrices and would like to know how I
can column permute a sparse matrix. Here is a small example:
> M1 <-
> spMatrix(nrow=5,ncol=6,i=sample(5,15,replace=TRUE),j=sample(6,15,replace=TRUE),x=round_any(rnorm(15,2),0.001))
> M1
5 x 6 sparse Matrix of class "d
This looks much like the kinds of problems that function 'solve.QP' in
package 'quadprog' can handle. But if you want people to help you,
follow the posting guide and 'provide commented, minimal,
self-contained, reproducible code'. You need to show people what you
tried and how it failed (so
56 matches
Mail list logo
