Re: [R] R2WinBUGS results not different with different runs

Gregor Gorjanc wrote: toby909 at gmail.com writes: Is this a specialty with R2WinBUGS? Does it have something to do with the seed value? Isnt the seed value reset everytime I restart winbugs? I always have the same seed if I start WinBUGS multiple times. So you get exactly the same

[R] R2WinBUGS results not different with different runs

Hi All I dont know if anyone else has noticed the same thing, but with 2 subsequent runs of the same syntax, I am getting exactly the same results. I was expecting that results differ slighlty, say in the 4th or 5th decimal place. Is this a specialty with R2WinBUGS? Does it have something to do

[R] R2WinBUGS more updates after model did not converge

After running a model for a while and seeing that it did not converge yet, how can I continue to run, ie not starting anew, the model? I know if I manually/interactively use winbugs, this is possible anytime, but how can I do this in r2winbugs, so that my existing sim$sims.array and other

[R] R2WinBUGS awkward to use

Hi All Does anyone know if I can avoid to use the write.model() function below? I dont want to do this. Can't bugs() do that automatically for me just by specifying the 4th argument 'model'? Just I like I am also using the 'inits' object! If I use 'model' in the same way as I use 'inits' I am

[R] avoiding timconsuming for loop renaming identifiers

Hi All I was wondering if I can avoid a time-consuming for loop on my 60 obs dataset. school_id y 8 9.87 8 8.89 8 7.89 8 8.88 20 6.78 20 9.99 20 8.79 31 10.1 31 11 There are, say, 143 different schools

[R] constraining residual variance to zero in lme or lmer

Hi All I havent seen that lme or lmer allows to specify certain constraints (although I heard one could program ones own covariance structure). Specific constraints, including one on the residual variance can be specified in sas proc mixed. In my model I want to fix the residual variance to

[R] opening a file from within a zipfile that is online

Hi Reading the help for ?unz I was wondering if I can read data into R from within an zipfile that is on some website, like maybe: dtaa = read.table(unz(http://www.ats.ucla.edu/stat/examples/alsm/alsm.zip,Ch01pr19.dat;)) Thanks for letting me know if you came acros such a thing before. Toby

[R] fixed effects anova in lme lmer

Can lme or lmer fit a plain regular fixed effects anova? Ie a model without a random effect, or have there be at least one random effect in order for these functions to work? Trying to run such, (1) without specifying a random effect produces an error, (2) specifying that there is no random

Re: [R] how to specify starting values in varIdent() of lme()

Douglas Bates wrote: On 6/1/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: I was reading the help but just did not get how to specify starting values for varIdent() of the lme() function, although I managed to do it for corSymm(). Do I specify the values just as they are printed out in an

[R] how to specify starting values in varIdent() of lme()

I was reading the help but just did not get how to specify starting values for varIdent() of the lme() function, although I managed to do it for corSymm(). Do I specify the values just as they are printed out in an output, like c(1, 1.3473, 1.0195). Or do I need to take the residual and

[R] parallel processing an lme model

Hi All Had anyone of you seen if it is possible to split a large lme() job to be processed by multiple cpus/computers? I am just at the very beginning to understand related things, but does the lme() use solution finding functions like nlm() mle(), and I-dont-know-what-else of the standard R

Re: [R] covariance question which has nothing to do with R

while my other program is running. The reference I mentioned previously addresses exactly this. Snijders and Bosker's Multilevel Analysis book on page 31 and 33, section 3.6.2 and 363 discuss this. When you say that the Xs are correlated then you would need to say according to which

Re: [R] in unix opening data object created under win

On Unix R's version is 2.3.1 and on PC its 2.4.1. I dont have the rights to install newer version of R on Unix. I tried different upload methods. No one worked. On Unix it looks as follows (dots to hide my userid): load(/afs/ir/users/../project/ps/data/dtaa) head(dtaa)

[R] in unix opening data object created under win

Hi All I am saving a dataframe in my MS-Win R with save(). Then I copy it onto my personal AFS space. Then I start R and run it with emacs and load() the data. It loads only 2 lines: head() shows only two lines nrow() als say it has only 2 lines, I get error message, when trying to use this data

Re: [R] Running an R script without running R

This gives you a GUI that runs R without you having to run it yourself on your machine: http://rss.acs.unt.edu/cgi-bin/R/Rprog You could upload your data to a webspace if you have one, and catch it from there in that service they provide. You could provide to whomever you like such a service

[R] saving datafreame object problem

Do I miss here something? dtaa = read.table(http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat;, sep=,) head(dtaa) # shows the data as it should be save(dtaa,dtaa,file=c:/dtaa) d = load(c:/dtaa) head(d) # all data is lost, it only shows [1] dtaa dtaa Thanks for

Re: [R] can I get same results using lme and gls?

Douglas Bates wrote: On 5/20/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls

[R] can I get same results using lme and gls?

Hi All I was wondering how to get the same results with gls and lme. In my lme, the design matrix for the random effects is (should be) a identity matrix and therefore G should add up with R to produce the R matrix that gls would report (V=ZGZ'+R). Added complexity is that I have 3 levels, so

[R] gls() error

Hi All How can I fit a repeated measures analysis using gls? I want to start with a unstructured correlation structure, as if the the measures at the occations are not longitudinal (no AR) but plainly multivariate (corSymm). My data (ignore the prox_pup and gender, occ means occasion):

Re: [R] plotting multilevel / lme lines

Yesterday, I tried to do exactly this, too. Below is my approach. Unfortunately, I did not find a textbook example against which I could verify my code. Hints, simplifications, and verification highly appreciated!!! ### random slope! model v = lme(POPULAR ~ SEX + TEXP, data=dta, random = ~ SEX

[R] translating sas proc mixed to lme()

Hi All I am trying to translate a proc mixed into a lme() syntax. It seems that I was able to do it for part of the model, but a few things are still different. It is a 2-level bivariate model (some call it a pseudo-3-level model). PROC MIXED DATA=psdata.bivar COVTEST METHOD = ml; CLASS

[R] controling omitted category in factor()

Hi All How can I specify which category R should omit when running a linear model with categorical predictors? I saw it omits the first category by default, but I like to have the 3rd category omitted. Thanks for your hints. Toby __

[R] long character string problem

Hi All I am having 2 very long character strings (550chars) and I want to put them as expressions together with c(). The problem is that I also get these double-quotes, as seen below in 'fct'. How can I remove these double-quotes? I tried as.name() but it did not work (because of size?). These

[R] piecing together statements (macro?)

Hi All I am pretty new to R but saw stata and sas's macro facilities and am looking for how such things work in R. I am trying to piece together a series of statements: n = 5 #want to have it dynamic with respect to n for (j in 1:n) { eval(paste(x, j, =x[, j, ], sep=)) } I want the

[R] autoload libraries at startup

Hi All I was wondering if there is a way I can specify in R that it should load libraries automatically at startup, so that I do not have to manually issue the command. Thanks Toby __ R-help@stat.math.ethz.ch mailing list

[R] long character string problem

Hi All I am having 2 very long character strings (550chars) and I want to put them as expressions together with c(). The problem is that I also get these double-quotes, as seen below in 'fct'. How can I remove these double-quotes? I tried as.name() but it did not work (because of size?). These

Re: [R] long character string problem

Yes parse(text=) is what I was looking for. thanks. (I also fixed the missing operator). I will feed this entire expression to the deriv() function. T Stephen Tucker wrote: I think you are looking for fct - c(parse(text=cum1),parse(text=cum2)) although you need to include operators

[R] D() and deriv() accessing components of returned expression

Hi All I like the D() function since it directly gives me the derivative. D() however does not take multiple arguments like c(x1, x2), but deriv() does, but deriv() I am having trouble accessing the actual derivative function. So strange, its sitting direct in front of me but I cant access it:

[R] does function deriv take a vector of expressions?

Hi All I was just wondering if the function deriv also takes a vector of expressions If someone knows I would appreciate letting me know. Thanks Toby __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do