Re: [R-br] GLM

2013-12-13 Por tôpico ana paula coelho madeira
Entendi, muito obrigada! Mas e quanto aos fatores que não aparecem na análise? Tenho 8 tempos e 4 cultivares e quando peço o summary do modelo o tempo 0 e a cultivar C1 não aparecem na análise. O que ocorreu? Como avaliar a significância desses fatores? Segue a saída: NP1 <- glm (Y ~ T+C, family

Re: [R-br] Otimização (boost + simplex)

2013-12-13 Por tôpico Ivan Bezerra Allaman
Debug a função para que tenha um diagnóstico melhor. library(debug) mtrace(simplex) mtrace(simplex1) O erro está acontecendo dentro a função simplex1 que é interna da função simplex. O erro ocorre em um while. Então, só você mesmo pra saber o porque do erro. (s,f,p) Allaman _

[R-br] Otimização (boost + simplex)

2013-12-13 Por tôpico Rodolfo Timoteo da Silva
Senhores, ao tentar otimizar a partir do arquivo abaixo: /library(boot)// //z = c(1,1,1,1,1,1,1)// //res2 = matrix(c(1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,1,1,1,1,1,0,0,0,1,1,1,1,1),7,7)// //simplex(a=z,A1=NULL,b1=NULL,A2=rbind(res2),b2=c(18,12,15,19,14,16,11),A3=

Re: [R-br] GLM

2013-12-13 Por tôpico Paulo Justiniano
Nao! o ponto "." siginifica todas as variaveis no data-frame alem da resposta (Y) Se o data frame contem apenas Y, T, e C entao as especificaqcoes sao iguais On Fri, 13 Dec 2013, Fátima Lima Paula wrote: Ana, o modelo com o ponto é o modelo nulo. O outro está ajustando pelo T e C. Não apare

Re: [R-br] GLM

2013-12-13 Por tôpico Fátima Lima Paula
Ana, o modelo com o ponto é o modelo nulo. O outro está ajustando pelo T e C. Não aparecem as categorias que são as de referência. Abs Em 13 de dezembro de 2013 13:26, ana paula coelho madeira < apcmade...@hotmail.com> escreveu: > Boa tarde pessoal. > > Tenho dados de contagem em função do tempo

[R-br] GLM

2013-12-13 Por tôpico ana paula coelho madeira
Boa tarde pessoal. Tenho dados de contagem em função do tempo e da cultivar. Tenho duvida na maneira como declaro o modelo. Qual é a diferença entre os modelos NP1 <- glm (Y ~ T+C, family = poisson, data=dados) e NP2 <- glm (Y ~ ., family = poisson, data=dados). Outra coisa, quando rodo o NP