search for: andrejk

Dear listers, I have two questions: (1) Is there a way in R to change the base-n of the calculations. I wnat to run some calculations either in binary (base-2) or base-4. Is there a way to specify that in R - to chnage from the decimal? (2) I also want to extract the digits from a larger number and store them as a vector. I could do it through converting top string, parsing the string and

...dDev Corr (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:andrejk@zrc-sazu.si] Sent: Monday, September 29, 2003 7:05 PM To: R-Help Subject: predicting values from the LME Dear listers, I experinced a problem prdicting the values using the LME with multilevel data. I have NA's in my dependent variable and the model is fitted only on the completed cases. I w...

Dear listers, I can't find the variance or se of the coefficients in a multilevel model using lme. I want to calculate a Chi square test statistics for the variability of the coefficients across levels. I have a simple 2-level problem, where I want to check weather a certain covariate varies across level 2 units. Pinheiro Bates suggest just looking at the intervals or doing a rather

Hi! I am a relatively new user of R and I use it to prepare my dissertation. I have come to some very usefull and specific libraries written for S-PLUS 4 and would like to use them in R. Is that possible? I just found out that one of these libraries has already been transfered to R, while 3 others have not. For the matter of beeing more exact I''m interested in the dealing with missing

Dear listers, I must first say that all the suggestions and help I got from the list so far was of great help to my work. I approach you with another question. Is it possible to generate an automatic specification of the formula for lme (or other models in R)? Let me clarify. I'm using models in a simulation run and can have variable number of variables in the data. I would like to generate

...---- Message: 19 Date: Thu, 26 Jun 2003 10:54:50 -0400 (EDT) From: "J.R. Lockwood" <lockwood at rand.org> Subject: RE: [R] within group variance of the coeficients in LME To: Harold Doran <hdoran at nasdc.org> Cc: R-Help <r-help at stat.math.ethz.ch>, Andrej Kveder <andrejk at zrc-sazu.si> Message-ID: <Pine.LNX.4.33.0306261037560.6707-100000 at penguin.rand.org> Content-Type: TEXT/PLAIN; charset=US-ASCII > > Dear listers, > > I can't find the variance or se of the coefficients in a multilevel model > using lme. > The component...

Hi! I'm pretty raw when working with the R models (linear or not). I'm wondering has anybody worked with the NLME library and dichotomous outcomes. I have a binary outcome variable that I woul like to model in a nested (multilevel) model. I started to fit a logistic model to a NLS function, but could not suceed. I know there are better ways to do it in R with either the LRM or GLM wih

Dear listers, I experinced a problem prdicting the values using the LME with multilevel data. I have NA's in my dependent variable and the model is fitted only on the completed cases. I want to estimate the predicted values for the rest of the data (those cases with missing dep. variable) I extracted a subset from the original file containing the variables used in the model as well as the

Hi listers, In light with the recent discussion on the optimizing the use of memory in straneous proceudres i present you m problem, and hope to some additional ideas. I'm running a simualtion that in each step uses quite an amount of memory (but not exceedingly) - just to give you an idea - I create a pseudo population (n=1000, m=3) run lme and lm model and multiply impute (M=5) and do the

Dear listers, I would like to get the bootstrap estimates form my lme model. I have an HLM (multillevel) 2-level model with the dichotomous outcome. I used glmmPQL procedure. However I have a problem since I have a rather unbalanced proportion (90-99% of events, i.e. ones and only 1-10% of nonevents - zeros) although sample sizes are not that small between 500 and 1000 I get pretty weird and

Dear listers, First let me appologize if the same mail arrives multiple times. Recently I had some probelms sending my e-mails to the list. I encountered a problem when running glmmPQL procuedure doing multilevel modeling with a dichotomous outcome. Those are the two error messages I usually get: Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3)