關於Probit Analysis R的警告消息

Question

我正在嘗試對一些死亡率測試進行概率分析以獲得LT50。 該數據集是針對一只昆蟲的。 在不同的時間，狀態要么是死的（1）要么是活的（0）。

> mbugdata
   state time
1      0   10
2      0   20
3      1   30
4      1   40
5      1   50
6      1   60
7      1   70
8      1   80
9      1   90
10     1  100

y<-cbind(mbugdata$state)
x<-cbind(mbugdata$time)
myprobit <- glm(state~time,family=binomial(link="probit"),data=mbugdata)

但是我不斷收到這些警告消息

Warning messages:
1: glm.fit: algorithm did not converge 
2: glm.fit: fitted probabilities numerically 0 or 1 occurred 

有人可以幫我嗎？ （我是R新手）

Answer 1

這不是R問題，而是統計問題。 您有過度安裝，分離等問題。

您應該修改統計模型規格以避免該錯誤。

如果您可以提供數據副本，那么我可以舉例說明，但是最好在CrossValidated上寫一篇文章，詳細介紹您的數據和研究問題。

還有一些罐頭方法可以避免此問題，例如使用bestglm軟件包，其中涉及自動選擇統計模型。

install.packages("bestglm")
require(bestglm)

這是一個用法示例 。