RE: [R] how to use try()

Just check the class() of the lo. You could try something like this... for(i in 1:20){ lo - try(nls(y~y0+a/(1+(x/x0)^b),start=list(y0=0.1,a=a0,x0=x00,b=-8.1))) if (class(lo)==try-error) beta[i]=NA else beta[i] - lo$m$getPars()[4] } Hth Andrej _ Andrej Kveder, M.A. researcher

[R] bootstrap estimates for lme

? Are there any special bootstrap tools for that? thanks for all the help best Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 [[alternative HTML version deleted

[R] memory optimization and use of recursion

, 512MbRam. Thanks for the suggestions... Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 [[alternative HTML version deleted

RE: [R] 2 questions regarding base-n and identifing digits + the source problem description

think there was a thread on the list just a short while ago. If anybody has a hint to follow, i would be more then glad to try it out... Thanks again Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440

[R] 2 questions regarding base-n and identifing digits

and store them as a vector. I could do it through converting top string, parsing the string and converting back. Is there another way. It would look something like that: a-1234 a [1] 1234 b-foo(a) b [,1] [,2] [,3] [,4] [1,]1234 Thanks! Andrej _ Andrej Kveder, M.A

RE: [R] r editors

statement in R, as was suggested by Gabor. I have nothing but the best to tell about the TextPad. Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 -Original Message

[R] FW: error predicting values from the LME

to be indispensible for using lme. hope this helps. spencer graves Andrej Kveder wrote: HI all, I might add some more information in order to possibly solve my problem. I'm really stuck and no obvious solutions do the trick. I'm using R 1.7.1 on Windows 2000 with the packages regurarly updated

[R] FW: error predicting values from the LME

(Intercept) 0.385620605 (Intr) v11 0.003147431 -0.048 Residual 0.450012367 Number of Observations: 729 Number of Groups: 50 If this give you some more insight to my problem. I would reallly appreciate any suggestion. Thanks Andrej -Original Message- From: Andrej Kveder [mailto:[EMAIL

[R] predicting values from the LME

? Thanks for all the suggestions. Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 [[alternative HTML version deleted

[R] glmmPQL probelm

. Thanks Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 [[alternative HTML version deleted]] __ [EMAIL PROTECTED

RE: [R] within group variance of the coeficients in LME

:[EMAIL PROTECTED] Behalf Of Douglas Bates Sent: Monday, June 30, 2003 6:09 PM To: J.R. Lockwood Cc: Harold Doran; R-Help; Andrej Kveder Subject: Re: [R] within group variance of the coeficients in LME J.R. Lockwood [EMAIL PROTECTED] writes: Dear listers, I can't find the variance

RE: [R] within group variance of the coeficients in LME

looking at the best possible solution to do it quantatively. Andrej -Original Message- From: Harold Doran [mailto:[EMAIL PROTECTED] Sent: Wednesday, July 02, 2003 3:20 PM To: Andrej Kveder; Douglas Bates; J.R. Lockwood Cc: R-Help Subject: RE: [R] within group variance of the coeficients

[R] within group variance of the coeficients in LME

to them. Does anybody have an idea how to get to them? I would really appreciate any suggestion. Andrej _ Andrej Kveder, M.A. researcher Institute of Medical Sciences SRS SASA; Novi trg 2, SI-1000 Ljubljana, Slovenia phone: +386 1 47 06 440 fax: +386 1 42 61 493 [[alternative

[R] simulating correlated vars

#proportion r12-0.4 #correlation x1-rbinom(n,1,p) x1-(x1-p)/sqrt(p*(1-p)) x2-rnorm(n,0,1) r-array(1,c(2,2)) r[1,2]-r12 r[2,1]-r12 r.ch-chol(r) d.1-cbind(x1,x2) d.1-d%*%r.ch d.1[,1]-round((d.1[,1]*sqrt(p*(1-p)))+p) #d.1[,1]-round(d.1[,1]+p)*sqrt(p*(1-p)) cor(d.1) mean(d.1[,1]) Andrej _ Andrej

[R] a question regarding s-plus libraries and R

others have not. For the matter of beeing more exact I'm interested in the dealing with missing survey data and Joe Shafer has the routines in S posted on his web site. Could you advise how I could get them to work in R? Thank you very much for the raply. Best regards Andrej _ Andrej Kveder