全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3464 3
2012-11-02
library(SparseM)
library(quantreg)
RD1=read.table("18-18.txt",header=TRUE)
RD1
Fixed.wired.Internetsubscriptionsper100inhabitants=RD1[1:78,
Fixedtelephonelinesper100inhabitants=RD1[1:78,2:2]
population_density=RD1[1:78,3:3]
GDP_PER_CAPITA=RD1[1:78,4:4]
fit1=rq(Fixed.wired.Internetsubscriptionsper100inhabitants~Fixedtelephonelinesper100inhabitants+population_density+GDP_PER_CAPITA,tau=c(0.05,0.25,0.5,0.75,0.95),data=RD1)
fit1


> confint(fit1,level=0.95)
Error in UseMethod("vcov") :
no applicable method for 'vcov' applied to an object of class "rqs"


我哪裡打錯了嗎?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-11-2 08:24:43
summary(fit1)没有结果吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-2 16:13:15
> summary(fit1,se="boot")

Call: rq(formula = Fixed.wired.Internetsubscriptionsper100inhabitants ~
Fixedtelephonelinesper100inhabitants + population_density +
GDP_PER_CAPITA, tau = c(0.05, 0.25, 0.5, 0.75, 0.95),
data = RD1)

tau: [1] 0.05

Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) -2.28359 1.38271 -1.65153 0.10287
Fixedtelephonelinesper100inhabitants 0.23708 0.09440 2.51146 0.01421
population_density 0.00086 0.00127 0.67946 0.49897
GDP_PER_CAPITA 0.00007 0.00009 0.74140 0.46080

Call: rq(formula = Fixed.wired.Internetsubscriptionsper100inhabitants ~
Fixedtelephonelinesper100inhabitants + population_density +
GDP_PER_CAPITA, tau = c(0.05, 0.25, 0.5, 0.75, 0.95),
data = RD1)

tau: [1] 0.25

Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) -1.91703 1.35011 -1.41990 0.15984
Fixedtelephonelinesper100inhabitants 0.26346 0.11604 2.27053 0.02609
population_density 0.00055 0.00143 0.38272 0.70302
GDP_PER_CAPITA 0.00019 0.00011 1.75389 0.08359

Call: rq(formula = Fixed.wired.Internetsubscriptionsper100inhabitants ~
Fixedtelephonelinesper100inhabitants + population_density +
GDP_PER_CAPITA, tau = c(0.05, 0.25, 0.5, 0.75, 0.95),
data = RD1)

tau: [1] 0.5

Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) -1.45516 0.77641 -1.87422 0.06485
Fixedtelephonelinesper100inhabitants 0.36948 0.10434 3.54092 0.00069
population_density 0.00022 0.00121 0.18305 0.85526
GDP_PER_CAPITA 0.00024 0.00011 2.15242 0.03463

Call: rq(formula = Fixed.wired.Internetsubscriptionsper100inhabitants ~
Fixedtelephonelinesper100inhabitants + population_density +
GDP_PER_CAPITA, tau = c(0.05, 0.25, 0.5, 0.75, 0.95),
data = RD1)

tau: [1] 0.75

Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) -0.59898 0.28577 -2.09603 0.03950
Fixedtelephonelinesper100inhabitants 0.40807 0.10233 3.98780 0.00016
population_density 0.00151 0.00114 1.32090 0.19061
GDP_PER_CAPITA 0.00024 0.00011 2.20931 0.03025

Call: rq(formula = Fixed.wired.Internetsubscriptionsper100inhabitants ~
Fixedtelephonelinesper100inhabitants + population_density +
GDP_PER_CAPITA, tau = c(0.05, 0.25, 0.5, 0.75, 0.95),
data = RD1)

tau: [1] 0.95

Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) -0.30633 2.28814 -0.13388 0.89386
Fixedtelephonelinesper100inhabitants 0.57010 0.16291 3.49956 0.00079
population_density 0.00198 0.00119 1.65866 0.10142
GDP_PER_CAPITA 0.00016 0.00019 0.86980 0.38722
> confint(fit1)
Error in UseMethod("vcov") :
no applicable method for 'vcov' applied to an object of class "rqs"
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-2 16:14:03
要怎麼跑出confidence interval ??
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群