search for: lm4

I have R under Linux Kubuntu Hardy Heron 8.04 LTS. I am trying to install package 'lme4' and get the following error: ---- * installing *source* package ?lme4? ... ** libs gcc -std=gnu99 -I/usr/share/R/include -I/usr/share/R/include -I"/usr/lib/R/library/Matrix/include" -I"/usr/lib/R/library/stats/include" -fpic -g -O2 -c init.c -o init.o gcc -std=gnu99

...sed so far is Does sports predict GPA? > lm1<-lm(GPA~sports) > summary(lm1) Does sports predict income? > lm2<-lm(income~sports) > summary(lm2) Does GPA predict income? > lm3<-lm(income~GPA) > summary(lm3) Does sports predict income after accounting for GPA? > lm4<-lm(income~GPA+sports) > summary(lm4) Can someone let me know is the above is correct? I am not sure if to keep all four regressions or only 1 and 4. Also, I need to plot the data on a graph or table. Can anyone suggest how to do this? Thank you in advance :) -- View this message in cont...

Dear all, First, I would like to thank you for your immense work. My question is about a frequent topic which I am not able to solve - even after hours of search in the mailing lisy. I would like to analyse random-effects (and fixed-effects)models of longitudinal / panel data with sampling weights. I have an unbalanced panel of different individuals in 5 years and income data as well as their

...s) for(itn in 1:nitns){ g<-rbinom(nsims,1,0.5) b<-rnorm(nsims,0,1)*10 rn<-rnorm(nsims,0,1)*10 a<-b*r+rn*(1-r^2)^0.5 a<-round(a)+50 a<-a-g*5 b<-round(b)+50 abs.2<-function(x) ifelse(x<1,1,x) b<-abs.2(b) c<-b-a p<-c/b lm1<-lm(a~g) lm2<-lm(c~g) lm3<-lm(p~g) lm4<-lm(a~b+g) gr<-c(g,g) occasion<-rep(0:1,c(nsims,nsims)) occ<-occasion-0.5 ppd<-c(b,a) h<-rep(0,nsims) mb<-mean(b) bppd<-b-mb bappd<-c(h,bppd) occgr<-occ*gr subject<-c(1:nsims) sub<-c(subject,subject) library(nlme) lm5<-lme(ppd~occ+gr+occgr,random=~1|sub) lm5f&...

Hello, I am using a simple linear model and I would like to get an AIC value. I came across both AIC() and extractAIC() and I am not sure which is best to use. I assumed that I should use AIC for a glm and extractAIC() for lm, but if I run my model in glm the AIC value is the same if I use AIC() on an lm object. What might be going on? Did I interpret these functions incorrectly? Thanks,

Dear friends. I use nls() and encounter the following puzzling problem: I have a function f(a,b,c,x), I have a data vector of x and a vectory y of realized value of f. Case1 I tried to estimate c with (a=0.3, b=0.5) fixed: nls(y~f(a,b,c,x), control=list(maxiter = 100000, minFactor=0.5 ^2048),start=list(c=0.5)). The error message is: "number of iterations exceeded maximum of

...e plm package to examine a longitudinal dataset. Each individual in this dataset has a longitudinal weight for the probability that he or she remains in the sample. Unfortunately, I have not found an argument to use weights in the plm function? I tried ?weights=? like in standard lm or in nlme or lm4 but it does not work. I asked the maintainer but I have not received a reply yet. Therefore, I would like to ask you whether I can use plm with longitudinal weights or if I have to switch to other packages dealing with longitudinal data? Thank you in advance! Best regards, Arne -- View th...

la libreria lm4 cambio con la versón 3 de R. Fijate en las versiones Prof. Julio Di Rienzo Estadística y Biometría FCA- U.N. Córdoba IBS-RARG President http://sites.google.com/site/juliodirienzo "Biometry, the active pursuit of biological knowledge by quantitative methods." (R.A. Fisher, 1948) 2013/10...

Estimados integrantes de la lista. Disculpas por posteo cruzado. Estoy ajustando un modelo con lmer (lm4). La variable respuesta es un índice (ADI) que se midió en 3 áreas diferentes en 4 estaciones climáticas diferentes, así mis efectos fijos son area y estaciones climáticas. Cada área tiene 12 puntos de muestreo (Punto). En cada punto de muestreo el índice se midió repetidamente a diferentes hora...

...los procesos abiertos por el sistema, si se > abren 4 rsession y los 2 procesadores + los 2 hilos se ponen en uso del > 100%. > > > Saludos > > > El 24 de octubre de 2013 12:36, Julio Alejandro Di Rienzo < > dirienzo.julio@gmail.com> escribió: > > la libreria lm4 cambio con la versón 3 de R. Fijate en las versiones >> >> Prof. Julio Di Rienzo >> Estadística y Biometría >> FCA- U.N. Córdoba >> IBS-RARG President >> http://sites.google.com/site/juliodirienzo >> "Biometry, the active pursuit of biological >>...

...", ncpus=4L Pero el caso es que si miro los procesos abiertos por el sistema, si se abren 4 rsession y los 2 procesadores + los 2 hilos se ponen en uso del 100%. Saludos El 24 de octubre de 2013 12:36, Julio Alejandro Di Rienzo < dirienzo.julio@gmail.com> escribió: > la libreria lm4 cambio con la versón 3 de R. Fijate en las versiones > > Prof. Julio Di Rienzo > Estadística y Biometría > FCA- U.N. Córdoba > IBS-RARG President > http://sites.google.com/site/juliodirienzo > "Biometry, the active pursuit of biological > knowledge by quantitative metho...

Hola. Algo raro me pasa con la librería boot y con lme4. En una misma máquina, core i5. Un bootstrap R=500 utilizando los 4 hilos de procesamiento tarda 1 hora con xubuntu 12.04 y R 2.14 y 3 con xubuntu 13.10 y R 3.0.1 Alguna idea de que puede estar pasando?? [[alternative HTML version deleted]]

...happen to be all true or false. # # God I wish I had taken more then a Trig course. Where was I? # # Yes that phantom 'i' you see is that this is all in a big loop # for 6 possible models # lm1 = all regressors w/ intercept # lm2 = lm1 stepwise removal # lm3 = all regressors wo/ intercept # lm4 = lm3 stepwise removal # lm5 = Hand Tuned # lm6 = lm5 stepwise removal ################################################### toPurge=c() for(k in names(newmcReg[[i]])) { print (paste("check to see if",k,"is a useless regressors for model",i)) if(sum(newmcReg[[i]][k][,1])==length...

Hola a todos. Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo mixto. Tengo varias variables en mi data.frame, unas son a nivel individual y otras a nivel de comarcas. Listo algunas. ingre_6 : Ingresos (nivel individual) iscs_a : un indicador sintético resumen de otras variables, calculado mediante componentes principales. sau_com : superficie agraria útil de cada una

Dear arun and all R users, I will first of all try to simply define my issue.. I have data in the following format Year Discharge dd/mm/yyyy x .. … … … There are some NA values in the discharge which I would like to predict by using “predict command”. I cant figure out the way to write the coding for that. Could you please help me on that??? I have also ,written

...1,mean,na.rm=T))/2 sf12.1 <- rep(sf12[,1],4) QOL.numbers.sm <-apply(QOL.m[,5:12],1,mean,na.rm=T) worry.numbers.sm <-apply(worry.m[,5:12],1,mean,na.rm=T) trust.numbers.sm <-apply(trust.m[,5:12],1,mean,na.rm=T) comfort.numbers.sm <-apply(comfort.m[,5:12],1,mean,na.rm=T) lm4 <- lm(deviat.logo.sm ~ stai.score+lot.score+educ+denial.score+sex+truste.score+minrisk.numbers.sm+exagg.numbers.sm+QOL.numbers.sm +worry.numbers.sm+trust.numbers.sm+comfort.numbers.sm+mss.m.score+mss.b.score) anova(lm4) #lm5 <- lm(deviat.logo.sm ~ stai.score+lot.score+educ+sex+minrisk.numbe...

.....p [020] 00 00 00 00 00 10 00 94 D1 D3 7D 74 75 AC 39 1A ........ ..}tu.9. [030] BA B8 A7 F8 45 7A AE 9B 7B FD 01 DF D3 61 9C 29 ....Ez.. {....a.) [040] AE AB 10 58 0A 25 F8 02 8A EA 98 FE BB 09 5D CA ...X.%.. ......]. [050] 66 A3 0A 10 64 54 A7 AE E2 4C 6D 34 20 B9 CC F0 f...dT.. .Lm4 ... [060] E1 2F 80 6B B4 3E DE 5E 1F E6 55 F3 D4 78 58 0F ./.k.>.^ ..U..xX. [070] B5 8D 6C 13 EF 5C 24 0D 63 49 70 0E D9 56 A8 72 ..l..\$. cIp..V.r [080] FB 60 D8 EA 82 41 51 D7 D8 E0 F2 2E 4E A5 AB 44 .`...AQ. ....N..D [090] 06 00 00 01 00 00 00 77 00 7A 00 FF FF 00 00 31 .........