[R-br] (sem assunto)

[Luís Fernando Pirondi Junior]

Olá,

escrevi o seguinte script:

> A<-read.table("A.txt", header=TRUE)
> A
  Dias Massa
1   38   370
2  121  4200
3  208 13200
4  291 17100
5  381 18000
> plot(A$Dias,A$Massa,xlab="Dias",ylab="Massa", col="blue")
> modelo<-glm(Massa~Dias, data=A, family=binomial)
Mensagens de aviso perdidas:
In eval(expr, envir, enclos) : #sucessos não-inteiro em um glm binomial!
> summary(modelo)

Call:
glm(formula = Massa ~ Dias, family = binomial, data = A)

Deviance Residuals:
       1         2         3         4         5
-0.19282  -0.03793   0.15578   0.09105  -0.11083

Coefficients:
             Estimate Std. Error z value Pr(>|z|)
(Intercept) -3.737847   4.166859  -0.897     0.37
Dias         0.006566   0.014028   0.468     0.64

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 0.339765  on 4  degrees of freedom
Residual deviance: 0.083459  on 3  degrees of freedom
AIC: 5.1496

Number of Fisher Scoring iterations: 6

E obtive o gráfico em anexo.

Minha intenção era fazer uma regressão logística com os valores da Tabela
"A", no entanto a curva não parece muito ajustada aos parâmetros. O que
pode estar errado ?

-- 
*Luís Fernando Pirondi Junior*
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/bc376416/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: Logística.jpg
Tipo: image/jpeg
Tamanho: 31647 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140321/bc376416/attachment.jpg>