similar to: lattice

This is only partly a question about R, as I am not quite sure about the underlying statistical theory either. I have fitted a non-linear mixed-effects model with nlme. In the fixed part of the model I have a factor with three levels as explanatory variable. I would like to use Tukey HSD or a similar test to test for differences between these three levels. I have two grouping factors:

Dear all, I have tried to calculate a GLMM fit with lmer (lme4) and glmmPQL (MASS), I also used glm for comparison. I am getting very different results from different functions, and I suspect that the problem is with our dataset rather than the functions, but I would appreciate help in deciding whether my suspicions are right. If indeed we are attempting the wrong type of analysis, some

Dear all, I am confused. I have encountered some strange behaviour of gls > data(co2) > co2.y <- aggregate(co2,1,mean) > co2.y.data <- data.frame(co2=as.numeric(co2.y),year=seq(1959-1980,along=co2.y)) > co2.1.gls <- gls(co2~year+I(year^2), co2.y.data) > co2.2.gls <- update(CO2.1.gls, corr=corAR1()) > summary(CO2.2.gls) > plot(CO2.2.gls) plot shows standardized

JSS, which is at http://www.jstatsoft.org, publishes statistical software with manuals, as well as papers describing, comparing, and evaluating statistical software. Abstracts of all contributions accepted in JSS are published in printed form in JCGS. Although we publish software written in any language, we especially welcome software written in R. The general idea is that packages written in R

JSS, which is at http://www.jstatsoft.org, publishes statistical software with manuals, as well as papers describing, comparing, and evaluating statistical software. Abstracts of all contributions accepted in JSS are published in printed form in JCGS. Although we publish software written in any language, we especially welcome software written in R. The general idea is that packages written in R

Hi, I cant log into my Fedora core 2 box from another linux machine or an OSX machine It worked the very first time I tried then never since. I have downl?oaded the newest openssh sources, and done ./configure make make install as root but seems to still use the old openssh. I stopped and restarted the sshd any Ideas? cheers Dan I get this message with the -v flag on [daniel:~] dan%

Could someone explain the meaming of this error message from augPred: > augPred(area3.pen.nlme, primary=~day) Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels 1,2,3 not allowed for block > predict.nlme(area3.pen.nlme) does not produce an error. area3.pen.nlme was created with: > area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),

I am learning R, and although I have looked in the documentation, I may be asking something obvious. Sorry, if that is the case. In a split-plot design if there is a missing subunit summary gives me a table with two rows for the same factor, one in the error within section and one in the section using error between units. With no data missing the table is "normal". How does one

Hello everyone, is there a package/packages for factor analysis, particularly PCA? thanks, Katja Katja L??ytynoja Taitoniekantie 9 A 218 40 740 Jyv??skyl?? Finland tel.+35814 608058 cell.+35850 336 0174 kaloytyn at jyu.fi

I am new to R (and S) and I hope someone on the list can help me out. I find the following behaviour surprising, is it a bug or a feature? (looks like a bug to me) R Version 1.1.0, Windows (x86, Win32) June 15, 2000 OS Windows 95 cat.test <- function(){ summary(c(1:20)) cat("test\n") } summary produces no output, just "test" appears on RGUI's

Dear all, I am pleased to announce the CRAN release of a new package called 'KFAS' - Kalman filter and smoother. The package KFAS contains functions of multivariate Kalman filter, smoother, simulation smoother and forecasting. It uses univariate approach algorithm (aka sequential processing), which is faster than normal method, and it also allows mean square prediction error matrix Ft to

Dear all, I am pleased to announce the CRAN release of a new package called 'KFAS' - Kalman filter and smoother. The package KFAS contains functions of multivariate Kalman filter, smoother, simulation smoother and forecasting. It uses univariate approach algorithm (aka sequential processing), which is faster than normal method, and it also allows mean square prediction error matrix Ft to

Dear all, I'd like to discuss about a possible bug in function StructTS of stats package. It seems that the function returns wrong value of the log-likelihood, as the added constant to the relevant part of the log-likelihood is misspecified. Here is an simple example: > data(Nile) > fit <- StructTS(Nile, type = "level") > fit$loglik [1] -367.5194 When computing the

I installed wine on mandriva 2008.1 spring I installes a software called BeInSync. the install was succesfull but when I try to start the program I get a message that the computer needs to be restarted. restarting wine q booting the computer diidn't solve the problem. It seems as if the application don't recognize the restart and waits for the restart after the installtion. anyone has

Hello I am having this compilation error while running 'configure'. It fails to detect the xpm.h file which in fact is where it is supposed be. I am running RH70 and xpm&xpm-dev rpms are installed. Codeweavers FAQ offers a solution to this problem: http://wine.codeweavers.com/fom-meta/cache/38.html ... "I changed 'configure' line 980 from: if mkdir conftestdir;

Vaya, polémica. ¿R puro es R base? Porque van ya por 6618 paquetes en CRAN http://cran.r-project.org/web/packages/index.html Y nadie sabe bien cuántos hay en github. Por cierto, ¿alguno conoce un código para hacer una estimación? El 06/05/2015 a las 16:09, Freddy Omar López Quintero escribió: 2015-05-06 11:00 GMT-03:00 Pedro Concejero Cerezo <pedro.concejerocerezo en

Hello, I''ve setuped an ethernet bridge on a debian sarge 3.1 with l7-filter + ipp2 shapper rules The server is a supermicro p4sci + celeron pentium 4 base 3ghz + 512Mb + 2 ethernet e1000 One interface is connected to a cisco catalyst switch The other interface is connected directly to a CMTS (a sort of router for cable modem) configured as bridge too. More than 20Mbps of bandwith

Y con fread.. directamente sobre el ".rel"... > datIn <- fread("DRUG-AE.rel") > dim(datIn) [1] 6821 8 Saludos, Carlos Ortega www.qualityexcellence.es El 29 de septiembre de 2015, 18:45, Pedro Concejero Cerezo < pedro.concejerocerezo en telefonica.com> escribió: > No tiene cabecera. Tendrás que ponerle nombre a las variables > Prueba este código,

Hola, Las presentaciones de la sesión sobre ecuaciones estructurales (SEM) así como los vídeos -como ya anunció Carlos ayer- están subidos a la página del Grupo de Madrid: http://r-es.org/Grupo+de+Inter%C3%A9s+Local+de+Madrid+-+GIL+Madrid Saludos! -- Pedro Concejero BigData - Analytics @ Telefónica I+D<http://www.tid.es/es/Paginas/default.aspx> Part of Telefónica

Hola de nuevo, el archivo .rel esta en una carpeta comprimida .zip, yo la descomprimo con win rar. Como no sea ese el motivo por el que no me da el fichero completo? El 30 de septiembre de 2015, 10:30, MªLuz Morales <mlzmrls en gmail.com> escribió: > Pues soy gafe entonces, > no lo entiendo: > > > datIn <- fread("C:\\Users\\iphealthMariluz\\Documents\\Proyecto