以下のコード予測を実行しようとすると、間隔引数がなしまたは信頼されるべきであるというエラーが表示されるため、分位回帰を使用して予測間隔を取得する方法があるかどうか疑問に思っていますか?
注: rq を lm に置き換えると、予測間隔を見つけることができますが、rq を使用して見つけたいと考えています。
#mycode:
library(quantreg)
mydata <- read.csv("C:\\Users\\gokul\\Desktop\\Book4.csv")
attach(mydata)
summary(mydata)
mod = rq(y~q) # y is my dependent and q is my independent variables
summary(mod)
predict(mod, data.frame(q), interval = "prediction",level = 0.10)
CSV ファイルのスニペット:
dput(head(mydata,10))
structure(list(y = c(71143L, 68061L, 66603L, 66907L, 69073L,
72901L, 77521L, 81728L, 84842L, 87877L), q = c(71416.79329, 68003.59226,
66533.66142, 66620.44529, 68640.60953, 72945.13676, 77743.82153,
81604.52442, 84887.47483, 87904.33486)), row.names = c(NA, 10L
), class = "data.frame")