我想将以下模型设置为仅考虑响应变量“Pres”的上 25%,而不是整个数据集。我在想 xmin 和 xmax 线需要改变吗?
aa2<-gam(Pres~s(Mud,bs="ps", k=15),family=Gamma(link=log),data=Antho)
xmin <- ceiling(min(Antho$Mud[Antho$Bin==1]))
xmax <- floor(max(Antho$Mud[Antho$Bin==1]))
Mudnew <- seq(from=xmin, to=xmax, by=0.1)
pred.dat <- data.frame(Mudnew)
names(pred.dat) <- "Mud"
pred.aa2 <- data.frame(predict.gam(aa2, pred.dat, se.fit=TRUE, type="response"))
pred.aa2.comb <- data.frame(pred.dat, pred.aa2)
names(pred.aa2.comb)
plot(fit ~ Mud, data=pred.aa2.comb, type="l", lwd=2, col=1, ylab="Density per 0.0132 m2", xlab="Mud content (%)")
感谢您的帮助
--
谢谢,但我想做后者,即明确地模拟“上四分位数”..这合理吗?我正在使用 mgcv 仅供参考。