Buenos días compañeros, al hacer el salto del software SAS a R, me he perdido unas cosas. He buscado información pero no he resuelto la duda, me he quedado demasiado abrumado por la cantidad de información al respecto. Tengo un conjunto de datos vibb: (Case factor fijo, source factor aleatorio y vib variable) ID case source vib 1 S 1 13.1 2 S 1 13.2 3 S 2 16.3 4 S 2 15.8 5 S 3 13.7 6 S 3 14.3 7 S 4 15.7 8 S 4 15.8 9 S 5 13.5 10 S 5 12.5 11 A 1 15.0 12 A 1 14.8 13 A 2 15.7 14 A 2 16.4 15 A 3 13.9 16 A 3 14.3 17 A 4 13.7 18 A 4 14.2 19 A 5 13.4 20 A 5 13.8 21 P 1 14.0 22 P 1 14.3 23 P 2 17.2 24 P 2 16.7 25 P 3 12.4 26 P 3 12.3 27 P 4 14.4 28 P 4 13.9 29 P 5 13.2 30 P 5 13.1 Quiero tomar el modelo vib~case*source, con source como factor aleatorio y hacer la separación de medias para cada factor (Tukey y LSD). En SAS programaria: PROC GLM; CLASS case source; MODEL vib=case|source; RANDOM source source*case/TEST; MEANS case/tukey e=source*case; MEANS source/tukey e=source*case; RUN; Aparte tengo que tener en cuenta el error de tipo III (Type III), cuando la matriz no esta balanceda, no en este caso, pero si que lo tendré en casos futuros. Muchas gracias por la atención. Ya no ser donde buscar Marçal
Hola Marcal, En R, el procedimiento de SAS (GLM) se llama igual "glm()" pero es una función. Empieza por help(glm) y podrás ajustar el modelo tal y como estás indicando para SAS. Y como estás en la transición de SAS a R, seguramente te hará falta la ayuda de este libro: http://www.amazon.com/SAS-SPSS-Users-Statistics-Computing/dp/0387094172/ref=sr_1_1?ie=UTF8&qid=1298747173&sr=8-1 Que afortunadamente para unos y lamentablemente para otros si buscas por ahí, podrás conseguir el pdf. Saludos, Carlos Ortega www.qualityexcellence.es 2011/2/25 <marcal.plans@upc.edu>> Buenos días compañeros, > > al hacer el salto del software SAS a R, me he perdido unas cosas. He > buscado información pero no he resuelto la duda, me he quedado demasiado > abrumado por la cantidad de información al respecto. > > Tengo un conjunto de datos vibb: (Case factor fijo, source factor aleatorio > y vib variable) > ID case source vib > 1 S 1 13.1 > 2 S 1 13.2 > 3 S 2 16.3 > 4 S 2 15.8 > 5 S 3 13.7 > 6 S 3 14.3 > 7 S 4 15.7 > 8 S 4 15.8 > 9 S 5 13.5 > 10 S 5 12.5 > 11 A 1 15.0 > 12 A 1 14.8 > 13 A 2 15.7 > 14 A 2 16.4 > 15 A 3 13.9 > 16 A 3 14.3 > 17 A 4 13.7 > 18 A 4 14.2 > 19 A 5 13.4 > 20 A 5 13.8 > 21 P 1 14.0 > 22 P 1 14.3 > 23 P 2 17.2 > 24 P 2 16.7 > 25 P 3 12.4 > 26 P 3 12.3 > 27 P 4 14.4 > 28 P 4 13.9 > 29 P 5 13.2 > 30 P 5 13.1 > > Quiero tomar el modelo vib~case*source, con source como factor aleatorio y > hacer la separación de medias para cada factor (Tukey y LSD). En SAS > programaria: > > PROC GLM; > CLASS case source; > MODEL vib=case|source; > RANDOM source source*case/TEST; > MEANS case/tukey e=source*case; > MEANS source/tukey e=source*case; > RUN; > > Aparte tengo que tener en cuenta el error de tipo III (Type III), cuando la > matriz no esta balanceda, no en este caso, pero si que lo tendré en casos > futuros. > > Muchas gracias por la atención. Ya no ser donde buscar > > Marçal > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >[[alternative HTML version deleted]]
Marçal, creo que lo siguiente, contesta en parte a tu pregunta: require(lme4) datos=read.table(file="datos.txt",header = TRUE) #adjunto el fichero de datos al mensaje modelo=lmer(vib~case+(case|source),data=datos) summary(modelo) Un saludo. Olivier -- ____________________________________ Olivier G. Nuñez Email: onunez en iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es ____________________________________ ------------ próxima parte ------------ An embedded and charset-unspecified text was scrubbed... Name: datos.txt URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20110226/dec6a518/attachment.txt> ------------ próxima parte ------------ El 25/02/2011, a las 18:36, marcal.plans en upc.edu escribió:> Buenos días compañeros, > > al hacer el salto del software SAS a R, me he perdido unas cosas. > He buscado información pero no he resuelto la duda, me he quedado > demasiado abrumado por la cantidad de información al respecto. > > Tengo un conjunto de datos vibb: (Case factor fijo, source factor > aleatorio y vib variable) > ID case source vib > 1 S 1 13.1 > 2 S 1 13.2 > 3 S 2 16.3 > 4 S 2 15.8 > 5 S 3 13.7 > 6 S 3 14.3 > 7 S 4 15.7 > 8 S 4 15.8 > 9 S 5 13.5 > 10 S 5 12.5 > 11 A 1 15.0 > 12 A 1 14.8 > 13 A 2 15.7 > 14 A 2 16.4 > 15 A 3 13.9 > 16 A 3 14.3 > 17 A 4 13.7 > 18 A 4 14.2 > 19 A 5 13.4 > 20 A 5 13.8 > 21 P 1 14.0 > 22 P 1 14.3 > 23 P 2 17.2 > 24 P 2 16.7 > 25 P 3 12.4 > 26 P 3 12.3 > 27 P 4 14.4 > 28 P 4 13.9 > 29 P 5 13.2 > 30 P 5 13.1 > > Quiero tomar el modelo vib~case*source, con source como factor > aleatorio y hacer la separación de medias para cada factor (Tukey y > LSD). En SAS programaria: > > PROC GLM; > CLASS case source; > MODEL vib=case|source; > RANDOM source source*case/TEST; > MEANS case/tukey e=source*case; > MEANS source/tukey e=source*case; > RUN; > > Aparte tengo que tener en cuenta el error de tipo III (Type III), > cuando la matriz no esta balanceda, no en este caso, pero si que lo > tendré en casos futuros. > > Muchas gracias por la atención. Ya no ser donde buscar > > Marçal > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es
.... Marçal,
una cosa más, para que la columna source sea tratada como un factor
(tu codificación de source es ambigua para R),
pon más bien:
datos=read.table(file="datos.txt",colClasses=c
("factor","factor","factor","numeric"),header
= TRUE)
y
anova(modelo)
para el análisis de varianza.
Un saludo. Olivier
--
____________________________________
Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es
____________________________________
El 26/02/2011, a las 20:31, Olivier Nuñez escribió:
> Marçal,
>
> creo que lo siguiente, contesta en parte a tu pregunta:
>
> require(lme4)
> datos=read.table(file="datos.txt",header = TRUE) #adjunto el
> fichero de datos al mensaje
> modelo=lmer(vib~case+(case|source),data=datos)
> summary(modelo)
>
> Un saludo. Olivier
>