我使用以下数据:
Dose Cells Dicentrics
0.1 328 2
0.2 723 1
0.5 513 4
1.0 369 4
1.5 280 7
2.0 361 19
3.0 238 26
我正在使用 read.table
和 attach
加载表 . 我被建议使用模型Dicentrics = Cells(b1 b2剂量b3剂量^ 2) . 因此,我适合一个glm
fit1 <- glm(Dicentrics ~ Cells+Cells*Dose+Cells*(Dose^2)-1, family=poisson(link="identity"))
现在,我想将其可视化并将其绘制在Dicentric / Cells~Dose Plot之上 .
我想在这个情节上绘制我的预测线 .
当我使用 predict(fit1)
时,我明白了
predict(fit1, type="response") # predicted values
1 2 3 4 5
1.124833 2.521284 2.677220 4.906738 7.802772
6 7
17.662902 26.304251