[R] Consulta sobre GLM-log linear

Rui Barradas ruipbarradas at sapo.pt
Fri Jun 8 12:30:52 CEST 2012


Ola, outra vez.

Desculpa, está errado!

É

a[[2]]
a[["Deviance"]]  # claro!

com a[[1]] ou a[["Df"]] graus de liberdade. Os p-values estão correctos.


Rui Barradas


Em 08-06-2012 11:14, Rui Barradas escreveu:
> Ola,
>
> Eu sou português, não sei escrever castellano/espanhol mas compreendo.
> Aqui é melhor em inglês.
> Para responder à questão,
>
>
>
> modelo <- glm(...etc...)
> a <- anova(modelo, test="Chisq")
> a
> length(a)  # 'a' tem 5 elementos
> sapply(a, print)  # ver o objecto 'a'
> a[[4]]  # valores da estatistica, coluna 4
> a[["Resid. Dev"]] # o mesmo
>
> a[[5]]    # p-values
> a[["Pr(>Chi)"]]  # o mesmo
>
>
> Espero que ajude,
>
> Rui Barradas
>
> Em 08-06-2012 06:54, Luis Fernando García Hernández escreveu:
>> Estimados amigos,
>>
>> Estoy familiarizándome con los modelos lineales generalizados en R. Estoy
>> interesado en realizar un análisis lig linear y me gustaría saber cuáles
>> son  o como extraer los valores correspondientes al chi cuadrado en el
>> análisis para cada grupo y para las interacciones. Desde ya muchas
>> gracias
>> y disculpas si la pregunta es muy básica, adjunto los comandos que estoy
>> utilizando. Si hiciera falta el archivo de texto puedo adjuntarlo.
>>
>> Desde ya mcuhas gracias
>>
>> library(MuMIn)
>>
>> archivo<-read.table("C:\\Users\\R\\arch.txt",header=T)
>>
>> attach(C.narino)
>>
>> glm1=glm(Frec~Sexo*Hora*Comportamiento, family=poisson)
>>
>> summary(glm1)
>>
>> anova(glm1,test="Chisq")
>>
>> dredge(glm1, rank = "AIC")
>> .
>>
>>     [[alternative HTML version deleted]]
>>
>>
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list