作者ching0629 (Syameroke)
看板Statistics
标题Re: [程式] R分量回归结果显示
时间Thu Jan 3 20:56:47 2013
※ 引述《Ennail5566 (小妹)》之铭言：
: [软体程式类别]:
: R
: [程式问题]:分量回归结果输出
: [软体熟悉度]:
: 低(1~3个月)
: [问题叙述]:
: 目前我可以将我想输出分量回归的表显示出，
: 如:
: 0.01分量
: 变数 value std. error .....
: x1 . . .....
: x2 . . .....
: x3 . . .....
: 但我希望单独显示所有分量个别的系数
: 如:
: x1变数
: 分量 value
: 0.01 3
: 0.05 -2
: 0.10 -1
: 0.25 -10
: . .
: . .
: . .
: 请问我该输入什麽指令
: 感谢帮忙
A=c(0.01,0.05,0.1,0.25,0.5,0.75,0.9,0.95,0.99)
X1变数=matrix(NA,nrow=length(A),ncol=2)
for (分量 in 1:length(A)) {
X1变数[分量,1]=A[分量]
X1变数[分量,2]=rq(y~x,A[分量])$coefficients[2]
}
colnames(X1变数)=c("分量","value")
其中A是分量的总数(我想你应该看的懂)
第五行的y与x摆上自己想摆的,但x1(或是该自变数)要放第一个
由於你给你资讯很少,我甚至不知道你的自变项共有多少个,所以你就按照自己需要的呈现
吧
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 27.105.2.50
※ 编辑: ching0629 来自: 27.105.2.50 (01/03 20:59)
1F：→ ching0629:资料储存在"X1变数"这个物件内 01/03 21:00
2F：推 Ennail5566:谢谢 再请教一个问题 如果想输出标准误 01/03 21:31
3F：→ Ennail5566:该将coefficients改成什麽 01/03 21:31
4F：→ ching0629:手上没R,你用ls(rq(y~x,0.5))去试试看里面有哪些 01/04 11:01
5F：→ ching0629:应该会看到标准误是哪一个名字 01/04 11:01
6F：推 q6261901:std 但输出出来是null 我再研究看看好了 谢谢 01/04 12:45
7F：→ ching0629:我查了一下,rq(y~x,0.5)里面似乎没有标准误,而 01/04 15:36
8F：→ ching0629:summary(rq(y~x,0.5))里面也只有上下界,没有标准误 01/04 15:36
9F：→ ching0629:下界是这样叫的: 01/04 15:37
10F：→ ching0629:summary(rq(y~x,0.5))$coefficients[2,2] 01/04 15:38
11F：→ ching0629:这是X1的下界,你自己在for回圈中计算标准误吧 01/04 15:38
12F：推 q6261901:谢谢!! 01/04 18:52