使用Betareg生成的NaN

Question

当尝试对下面的数据集执行Beta回归时，我得到以下错误。谁能告诉我如何处理这个错误？

提前感谢

Error in quantile.default(x$residuals) : 
  missing values and NaN's not allowed if 'na.rm' is FALSE
In addition: Warning message:
In sqrt(v * (1 - hatvalues(object))) : NaNs produced

代码

model.beta = betareg(Percentage ~ Kuikenweek,data = dieet)
summary(model.beta)

数据

Jaar    Kuikenweek  Percentage  Locatie
2013    Week 1  0,4238  Markenje
2013    Week -1 0,8248  Markenje
2013    Week 2  0,5927  Markenje
2013    Week -2 0,3208  Markenje
2013    Week 3  0,8317  Markenje
2013    Week -3 0,7314  Markenje
2014    Week 1  0,4624  Markenje
2014    Week 2  0,5224  Markenje
2014    Week -2 0,3528  Markenje
2014    Week 3  0,7739  Markenje
2014    Week -3 0,1264  Markenje
2015    Week -1 0,4516  Markenje
2015    Week -2 0,4306  Markenje
2015    Week -3 0,5243  Markenje
2013    Week 1  0,5117  Scheelhoek
2013    Week -1 0,1490  Scheelhoek
2013    Week 2  0,3342  Scheelhoek
2013    Week -2 0,5069  Scheelhoek
2013    Week 3  0,7101  Scheelhoek
2013    Week -3 0,8019  Scheelhoek
2013    Week 4  0,7569  Scheelhoek
2013    Week 5  0,9096  Scheelhoek
2013    Week 6  0,9013  Scheelhoek
2014    Week 1  0,5403  Scheelhoek
2014    Week -1 0,2264  Scheelhoek
2014    Week 2  0,5775  Scheelhoek
2014    Week 3  0,7366  Scheelhoek
2014    Week -3 0,4086  Scheelhoek
2014    Week 4  0,5750  Scheelhoek
2014    Week 5  0,6623  Scheelhoek
2015    Week 1  0,7758  Scheelhoek
2015    Week -1 0,1517  Scheelhoek
2015    Week 2  0,7587  Scheelhoek
2015    Week -2 0,3247  Scheelhoek
2015    Week 3  0,7463  Scheelhoek
2015    Week -3 0,6325  Scheelhoek

Answer 1

我在这里重现你的错误：

# here are the data
devtools::source_gist("169bfa3a6c709fd2fd31c5bfa46648ee")
library(betareg)
# fix percentages 
dieet$Percentage <- gsub(",",".",dieet$Percentage)
dieet$Percentage <- as.numeric(dieet$Percentage)
model.beta = betareg(Percentage ~ Kuikenweek,data = dieet)
summary(model.beta)
## Error in quantile.default(x$residuals) : 
##  missing values and NaN's not allowed if 'na.rm' is FALSE
## In addition: Warning message:
## In sqrt(v * (1 - hatvalues(object))) : NaNs produced

出现这个问题是因为默认的残差类型"sweighted2“不适用于您的模型。所以：

summary(model.beta, type = "deviance")

至少能给你总结一下。有几种类型的残差可用，哪一种最适合您的目的，这更多的是一个统计问题。