[R-br] dúvidas sobre utilização de disitribuição beta

[Nicolay Cunha]

Olá a todos,

Tenho uma variável dependente (Y) que varia de 0 a 1 e quero testar em
relação a duas variáveis categóricas (X1 com três níveis e X2 com 7
níveis). A categórica X2 transformo em números para ganhar graus de
liberdade.

Os dados não atendem nenhum tipo de premissa de linearidade, portanto tive
que buscar outras alternativas para analisa-los.
Ao estudar sobre a melhor forma de ajustar esses dados, encontrei o artigo (
http://cran.r-project.org/web/packages/betareg/vignettes/betareg.pdf) sobre
o tema, e os exemplos dados aparentam encaixar com o meu caso.


No entanto meu Y tem valores inteiros (0 e 1) e a distribuiçãoa aparenta
não aceitar isso, para tal transformei todos os 0s em 0.0001 e 1s em
0.9999.  Ainda tenho dúvidas sobre a interpretação dos dados (inclusive o
papel do phi), utilização e escolha da melhor forma de gerar o modelo da
forma que estou propondo, logo gostaria de sugestões e orientações sobre
como proceder com essas análises.


Abaixo um CMR com um exemplo do que estou fazendo onde aparece uma mensagem
de aviso perdida ao estimar o modelo.


*"In betareg.fit(X, Y, Z, weights, offset, link, link.phi, type, control) :
  no valid starting value for precision parameter found, using 1 instead"*

#############################################################################

X1<-factor(rep(c("A", "B", "C"),50)) ## categoria 1
X2<- rep(sample(1:7), len=150 ) ## categoria 2
Y<-rep(sample(0:100)/100, len=150) ## minha variável dependente em proporção
Y[Y=='1']<-0.99999 # transformo 1 em 0.99999 pois a distribuição beta não
aceita 0 e 1, somente o intervalo entre eles
Y[Y=='0']<-0.00001 # transformo 0 em 0.00001

library(betareg)
beta1<-betareg(Y~X1+X2, link="loglog")
coef(beta1)
summary(beta1)

beta2<-betareg(Y~X1+X2|X2, link="loglog") ## adicionando X2 como um
regressor adicional
coef(beta2)
summary(beta2)

AIC(beta1, beta2,  k = log(nrow(data.frame(X1,X2,Y))))

##############################################################################



Agradeço a ajuda!
Nicolay







-- 
Nicolay Leme da Cunha

Biólogo, Mestre, Doutorando em Ecologia e Conservação
Universidade Federal de Mato Grosso do Sul, 79070-900
Campo Grande, MS, Brasil
E-mail: nicolaycunha em gmail.com
lattes.cnpq.br/5916316648872099




-- 
Nicolay Leme da Cunha

Biólogo, Mestre, Doutorando em Ecologia e Conservação
Universidade Federal de Mato Grosso do Sul, 79070-900
Campo Grande, MS, Brasil
E-mail: nicolaycunha em gmail.com
lattes.cnpq.br/5916316648872099
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121018/66797bb6/attachment.html>