Questao 107.5 Seriemas e Carcaras

classic Clássica list Lista threaded Em Árvore
7 mensagens Opções
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Questao 107.5 Seriemas e Carcaras

Alan Eduardo de Barros
Aqui estou tendo problemas ao determinar o p-value de p.ce

Esta dando 35%.Mas nao consigo encontrar o erro. Segue o codigo

#P ajeitar os dados
aves.c <- read.csv("aves_cerrado.csv", header=TRUE, sep=";", as.is=T)
aves.c$urubu[is.na(aves.c$urubu)==T] <- 0
aves.c$carcara[is.na(aves.c$carcara)==T] <- 0
aves.c$seriema[is.na(aves.c$seriema)==T] <- 0
aves.c$fisionomia[aves.c$fisionomia=="ce"] <- "Ce"
aves.c$fisionomia <- factor(aves.c$fisionomia, levels=c("CL","CC","Ce"))


#E o exercicio em si

aves_ce <- subset(aves, fisionomia == "Ce", select = c("fisionomia","carcara","seriema"))
aves_cc <- subset(aves, fisionomia == "CC", select = c("fisionomia","carcara","seriema"))
aves_cl<- subset(aves, fisionomia == "CL", select = c("fisionomia","carcara","seriema"))
mod.ce<- lm( seriema~carcara,data=aves_ce)
mod.cc<- lm( seriema~carcara,data=aves_cc)
mod.cl<- lm( seriema~carcara,data=aves_cl)

summary(mod.ce)
coef(mod.ce)
p.ce= 0.2446
coef.ce= -0.14116

summary(mod.cc)
coef(mod.cc)
p.cc= 0.6454
coef.cc=-0.08262712

summary(mod.cl)
coef(mod.cl)
p.cl= 0.1045
coef.cl=-0.4507576  
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

Diogro.
Não copie valores numéricos do console para o código, use indexação.
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

Alan Eduardo de Barros
Valeu!  

Deu um pouco de trabalho p conseguir ajustar, mas deu certo!

Dicas para quem estiver trabalhando nisso.

 p pede um só valor (descarta o intercepto e mantem a variavel de interesse, no meu caso carcara)

 para coeficientes sao 2 valores (mantem intercepto e a variavel de interesse)

Abraços!!!
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

GCImprota
Empaquei nos mesmos 35% e estou tentando indexar o p pelo "anova(mod.ce)$"Pr(>F)", mas não faço ideia de porque não está funcionando. Estou começando a questionar se o problema foram os meus "mods", mas acho que do jeito que eu fiz, dá no mesmo que fazer subsets para Ce, CC e CL:

aves <- read.csv("aves_cerrado.csv", sep=";", as.is=TRUE, header= TRUE)

head(aves)
table(aves$fisionomia)
aves$fisionomia[aves$fisionomia=="ce"] <- "Ce"
aves[is.na(aves)==TRUE] <- 0
aves$fisionomia <- as.factor(aves$fisionomia)

mod.ce <- lm(aves$seriema~aves$carcara, data= aves[aves$fisionomia=="Ce",])
mod.cc <- lm(aves$seriema~aves$carcara, data= aves[aves$fisionomia=="CC",])
mod.cl <- lm(aves$seriema~aves$carcara, data= aves[aves$fisionomia=="CL",])

p.ce <- anova(mod.ce)$"Pr(>F)"
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

GCImprota
Olhando de novo o código do Alan, resolvi tentar simplesmente "seriema~carcara" em lm() (ao invés de lm(aves$seriema~aves$carcara....) e deu certo. Não entendi direito o porque, mas agora rodou certo os meu valores de P.
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

Dylan Padilla
Esta mensagem foi atualizada em .
Em resposta à esta mensagem postada por GCImprota
unica coisa que falta é [1] depois do "Pr(>F)"

o "Pr(>F)" é um objeto com dois valores, dai você tem que escolher o primeiro desses valores. Ou seja, anova(seu modelo)$"Pr(>F)"[1]
Responder | Em Árvore
Abra essa mensagem na visão em árvore
|

Re: Questao 107.5 Seriemas e Carcaras

Diogro.
Em resposta à esta mensagem postada por GCImprota
Não estava funcionando pq quando vc usava as colunas do data.frame original na formula a indexação por fisionomia que vc fez no argumento data era ignorada. Quando vc passa o nome das colunas ele usa o argumento data pra procurar as colunas, e ai a indexação funciona.