r - 绘制使用 bbmle 包创建的非线性模型的残差

Question

我正在尝试获取我用 bbmle 构建的非线性模型的残差图，但不知道如何处理此任务。bble 包有一些关于残差的注释，但无法绘制直方图之类的东西。任何帮助将不胜感激

Answer 1

残差（）函数似乎工作。然后你可以做任何你想做的事情——直方图、qqplot、残差散点图与预测值（predict() 也有一个方法）。例如：

set.seed(1002)
lymax <- c(0,2)
lhalf <- 0
x <- runif(200)
g <- factor(rep(c("a","b"),each=100))
y <- rnbinom(200,mu=exp(lymax[g])/(1+x/exp(lhalf)),size=2)
dat <- data.frame(y,g,x)

fit3 <- mle2(y~dnbinom(mu=exp(lymax)/(1+x/exp(lhalf)),size=exp(logk)),
    parameters=list(lymax~g),
    start=list(lymax=0,lhalf=0,logk=0),
    data=dat)


par(mfrow=c(2,2))
hist(residuals(fit3))
qqnorm(residuals(fit3))
hist(residuals(fit3, type="response"))
qqnorm(residuals(fit3, type="response"))

还是我错过了重点？