similar to: about lme

Hi R-users: I've got a problem with lme. In Rail data,I try to model the next one lme(travel ~ Rail, data = Rail, random = ~ Rail | Rail) I want travel = Rail(i) + Rail(j) + epsilon(i,j) say, an effect fixed for every Rail and other for Rail(j) random, R says fewer observations or false convergence. Why can't I model an effect fixed and other random? Thanks in advance. My e-mail is

Hi R users: In R for windows. is there any package to work, in experimental design? With functions like in S-PLUS , fact.design and the other ones. How to implement in R, all theory about experimental design, I'm talking about ssType3, fact.design and the theory of orthogonal design, combinatory...... Thanks in advance and please excuse me, my "european english". My e-mail is

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

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 <-

Hi all I am trying to generate a normal unbalanced data to estimate the coefficients of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to estimate the variance components and their standard errors. Further, I am trying to use the likelihood ratio test to test H0: sigma^2_b = 0 (random effects variance component), and the t-test to test H0:mu=0 (intercept of the model Yij = mu

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

I am trying to use varcomp to obtain the variance partitioning across different nested levels of random effects (say x,y and z). I get the three variance components (for each of my along with an additional one called 'within' from varcomp output. I am using the 'scale total variance to 1' option and though the within component is small, it does form a part of what explains the

Hi try str(varcompobject) to see structure of this object. You can extract parts by standard R means. Cheers Petr > -----Original Message----- > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Sharada > Ramadass > Sent: Tuesday, August 8, 2017 3:33 PM > To: r-help at r-project.org > Subject: [R] how to extract individual values from varcomp? > >

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

Dear all I need to run an anova from a factorial model y_{ijk}=\alpha_i+\beta_j+(\alpha\beta)_{ij}+e_{ijk} and calculate type II and III sums of square, but I have an empty cells, so anova function from package car fail. (I believe) y<-c(7,13,6,10,8,11,8,3,7,5,65) a<-as.factor(c(1,1,2,2,3,3,3,1,1,1,2)) b<-as.factor( c(rep(1,7),rep(2,4)) ) table(b,a) # cell (2,3) is empty

Hi, I'm still having problems putting the variance components of my model in to a data frame, it is a continuation of this discussion, http://r.789695.n4.nabble.com/ANOVA-problem-td4609062.html, but now focussed on the problem of extracting variance components. I have got my mixed effects model now /narrow$line<-as.factor(narrow$line) rg.lmer <- lapply(split(narrow,

Hello, I am trying to use varcomp to decompose the variance across multiple nested levels on a lme object. I am able to successfully do this and when I view the varcomp object I can see the individual values / estimates for the variance at different levels. However, I want to be able to extract each of them separately, as I need to build a confidence interval using bootstrapping on the sample

Hello, I have some data covering contaminant concentrations in fish over a time period of ~35 years. Each year, multiple samples of fish were taken (with varying sample sizes each year). Ultimately, I want an estimation of the variance between years, and the variance within years + random effects. I used a linear mixed model to estimate these variances, but after reading a number of different

estimados, siempre habia usado data.frame, pero estoy tratando de aprender a usar data.table porque normalmente tengo que manejar una buena cantidad de datos ahora necesito hacer un calculo simple, para el que antes usaba ciclos FOR, pero quiero aprovechar las ventajas de data.table mi preguna es como podria calcular, por ejemplo, la columna 2 de un data.table llamada DT si el valor de la

Par evitar el warning.... > DT=data.table(col1=1:10) > DT[,col2:=diff(log(c(1,col1)))] > DT col1 col2 1: 1 0.0000000 2: 2 0.6931472 3: 3 0.4054651 4: 4 0.2876821 5: 5 0.2231436 6: 6 0.1823216 7: 7 0.1541507 8: 8 0.1335314 9: 9 0.1177830 10: 10 0.1053605 > ----- Mensaje original ----- De: "Carlos Ortega" <cof en

Hello R experts, I haven't found anything like the S function 'varcomp' as described in W.N. Venables & B.D. Ripley's 'Modern Applied Statistics ...' for R; does something comparable exist for R, or is planned for future releases? More generally, are there libraries with post-anova test procedures, like Student-Newman-Keuls? Or do those of you who frequently use

Hello: I need to prepare some detailed stats from asterisk, and I'm asked to show data I don't know how to obtain it: It's the 'final' number (don't know what's its name) In the stats I have to show the caller_id (I have it), the called_id (I have it) and the final number that actually accepted the call. In extensions.conf file, I try to pass the call to several

Hola. En un cierto directorio tengo archivos excel y de otros tipos. Cada uno de ellos corresponde a una encuesta a una empresa X. Me gustaria importar los datos de ciertas celdas de ciertas hojas de todos los archivos excel en una sola pasada. ¿Hay alguna forma de hacerlo? Yo se que la variable Z está siempre en la celda X de la hora tanto. Me gustaría ingresarlo todo en una matriz en que

Hola, buenas tardes: Necesitaba un poco de socorro con xtable, si fuese posible. Mi problema es que quiero 'formatear' específicamente algunas celdas de una tabla generada con xtable y no soy capaz. Detallo el problema de ejemplo: vv<-rbind(c(1,1,0,1),c(0,0,1,1)) vv[2,3]<-<-paste("\\textbf{",vv[2,3],"}",sep="") xtable(vv) > xtable(vv) % latex

Buenas tardes, Necesito contar cuántas comas hay en cada celda. Los datos tienen esta forma: Pros Comas Opinión 1 calidad, 1 Opinión 2 calidad, diseño, color, 3 Opinión 3 precio, accesibilidad, sienta bien, luminoso, 4 La variable Pros es la que tengo y la variable Comas es la