similar to: type II and III Sum square whit empty cells

Hi all, I have some difficulties with the lme function and so this is my problem. Supoose i have the following model y_(ijk)=beta_j + e_i + epsilon_(ijk) where beta_j are fixed effects, e_i is a random effect and epsilon_(ijk) is the error. If i want to estimate a such model, i execute >lme(y~vec.J , random~1 |vec .I ) where y is the vector of my data, vec.J is a factor object

HI, I would like to extract the variance components estimation in lme function like a.fit<-lme(distance~age, data=aaa, random=~day/subject) There should be three variances \sigma_day, \sigma_{day %in% subject } and \sigma_e. I can extract the \sigma_e using something like a.fit$var. However, I cannot manage to extract the first two variance components. I can only see the results in

Greetings, The following question has come up in an off-list discussion. Is it possible to construct a regular expression 'rex' out of two given regular expressions 'rex1' and 'rex2', such that a character string X matches 'rex' if and only if X matches 'rex1' AND X does not match 'rex2'? The desired end result can be achieved by logically combining

I'm fitting a GLMM to some questionnaire data. The structure is J individuals, nested within I areas, all of whom answer the same K (ordinal) questions. The model I'm using is based on so-called continuation ratios, so that it can be fitted using the lme4 package. The lmer function fits the model just fine, but using mcmcsamp to judge the variability of the parameter estimates produces

Hello everyone, I wanted to ask if there is an R-package to fit the following Poisson regression model log(\lambda_{ijk}) = \phi_{i} + \alpha_{j} + \beta_{k} i=1,\cdots,N (subjects) j=0,1 (two levels) k=0,1 (two levels) treating the \phi_{i} as nuinsance parameters. Thank you very much -- -Tony [[alternative HTML version deleted]]

Dear useRs, I've just encountered a serious problem involving the F-test being carried out in aov() and anova(). In the provided example, aov() is not making the correct F-test for an hypothesis involving the expected mean square (EMS) of a factor divided by the EMS of another factor (i.e., instead of the error EMS). Here is the example: Expected Mean Square

Hi all, I have a random effect model that can be written as y_{ij} = \beta +\alpha_i+ \epslion_{ij} where \alpha_i ~ N(0,\sigma^2_effect) and \espilon_{ij} ~ N(0,\sigma^2_error) and i compute res<-lme(data~1,data=Data,random=~1 | veci) The estimation of \sigma_error is given by res$sigma but I don't understand how i can find the estimation of \sigma _effect with the object res (I

Hello, I'm trying to fit a special kind of proportional odds model from: Whitehead et al. (2001). Meta-analysis of ordinal outcome using individual patient data. Statistics in medicine 20: 2243-2260. (model 2) The data are as follows: library(nlme) library(geepack) library(Design) library(MASS) options(contrasts=c("contr.SAS","contr.poly")) counts <-

Dear R users: I've a problem with lme function, when I want to model an unbalanced two-way anova, with 2 random factors say t and b. My two models are: model1- y(ijk) = beta+b(i)+t(j)+epsilon(ijk) model2- y(ijk)= beta+b(i)+t(j)+b:t(ij)+epsilon(ijk) beta overall mean effect The data.frame is X t b med celda 1 1 10 1 1 1 12 1 1 1 11 1 1 2 13 2 1

Dear All, I wonder if there is an efficient way to fit the generalized linear mixed model for multivariate outcomes. More specifically, Suppose that for a given subject i and at a given time j we observe a multivariate outcome Yij = (Y_ij1, Y_ij2, ..., Y_ijK). where Y_ijk is a binomial(n_ijk, p_ijk). One way to jointly model the data is to use the following specification: g(p_ijk) =

Hello, I have a simulated data structure in which students are nested within teachers, and with each student are associated two test scores. There are 20 classrooms and 25 students per classroom, for a total of 500 students and two scores per student. Here are the first 10 lines of my dataframe "d": studid tchid Y time 1 1 1 -1.0833222 0 2 1 1

Hi R-helpers, I'm new to time series modelling, but my requirement seems to fall just outside the capabilities of the arima function in R. I'd like to fit an ARMA model where the variance of the disturbances is a function of some exogenous variable. So something like: Y_t = a_0 + a_1 * Y_(t-1) +...+ a_p * Y_(t-p) + b_1 * e_(t-1) +...+ b_q * e_(t-q) + e_t, where e_t ~ N(0, sigma^2_t),

Dear R-users Coming from a proc mixed (SAS) background I am trying to get into the use of (n)lme. In this connection, I have some (presumably stupid) questions which I am sure someone out there can answer: 1) With proc mixed it is easy to get a hold on the estimated variance parameters as they can be put out into a SAS data set. How do I do the same with lme-objects? For example, I can see the

Hi friends, I'm beginning in R and I have simple question. I have this piece of my program and how you see, that's ok (whit > num<- 0.002) num<-0.002 # ok, but not when I change whit num<-0... ? factor1<-1; while(1) { if (num*factor1<1) factor1<-factor1*10 else { print("out ok!!"); break; } } [1] "out ok!!" but when I change (whit >

Dear List, I'm just teaching myself semi-parametric techniques. Apologies in advance for the long post. I've got observational data and a longitudinal, semi-parametric model that I want to fit in GAM (or potentially something equivalent), and I'm not sure how to do it. I'm posting this to ask whether it is possible to do what I want to do using "canned" commands

Hello... i have a reciept printer (SAMSUNG SRP-250) and i can't print correclty whit wine, some one know a driver for this kind of printers? i tried whit some drivers(whit CUPS) but this printing garbage ( $$$$$$$$$%%%%%%%%%%%%%%%&&&&&&&&&&&&&&&&). if i execute " cat file > /dep/lp0 " in the terminal that works

Hi List, I installed Samba 3.0.0beta3 on RH 8. I'm trying to make a shared directory so I can access from M$ machines using NT users ( from my PDC ). I got to join the Domain and I can see the users from my domain (in the Linux), I can even login as one... I used the RTFM "Samba-HOWTO-Collection" to do this... But I can't open the shared directory, not from any M$ client (or

Hello i am new to Wine and it is a very nice program how ever i have been trying to install a program name *Cain & Abel and every time i try to install it give me this error "C:\user\root\Temp\GLF1502.tmp could not be opened. Please check that your disk is not full and that you have access to the destination directory". I am not a big windows User but i remember that this kind of

I install WineHQ and I can't understand what they write me! Can you help me?? I use Linux Ubuntu 8.10 - the Intrepid Ibex [Image: http://img440.imageshack.us/img440/7590/wine.jpg ]

Hi. I'm entering in Samba's world and I need help with this. An Nt4.0 client was upgraded to XP PRO, and Win does not connect to domain (before upgrade worded fine), windows registry was modificated as samba howto explains. This is my smb.conf file of samba v. 2.0.7. Thanks, thanks, thanks. Dr Jekyll [global] ############## workgroup = X security = user status = yes