带置信区间的画图问题 [推广有奖]

xiaomuyoucun

发表于 2013-3-15 12:06:28 |显示全部楼层 |坛友微信交流群

epoh 发表于 2011-4-1 21:48
糊涂了,封掉a没发觉.
这问题,我再仔细看一下function rqss().

epoh 老师，我运行上边程序出现错误：

> xx1<- cbind(as(w,"matrix.diag.csr") %x% X,w %x% Z)
错误于cbind(deparse.level, ...) : args have differing numbers of rows
此外: 警告信息：
In as.matrix.csr(x, nrow = nrow, ncol = ncol, eps = eps) :
ncol*nrow indivisable by length(x)

使用道具举报

ywh19860616 发表于 2013-3-15 12:18:45 |显示全部楼层 |坛友微信交流群

xiaomuyoucun 发表于 2013-3-15 12:06
epoh 老师，我运行上边程序出现错误：

> xx1

哈哈，epoh老师太忙了
你可以直接利用rqpd 包了。

一份耕耘，一份收获。

使用道具举报

xiaomuyoucun

发表于 2013-3-15 12:45:03 |显示全部楼层 |坛友微信交流群

ywh19860616 发表于 2013-3-15 12:18
哈哈，epoh老师太忙了
你可以直接利用rqpd 包了。

这个程序之前你也用过，那请问怎么回事？是哪里出问题了，你能修改过来吗？

使用道具举报

xiaomuyoucun

发表于 2013-3-15 15:11:29 |显示全部楼层 |坛友微信交流群

ywh19860616 发表于 2011-4-1 13:54
我试加了，程序也是可以运行的，只是结果不一样
> library(SparseM)
> library(quantreg)

你好，我发现这一段程序总是运行不过去，你之前运行过没有错误吗？为什么我的总出现以下错误？
急于求解：
> xx1<- cbind(as(w,"matrix.diag.csr") %x% X,w %x% Z)
错误于cbind(deparse.level, ...) : args have differing numbers of rows
此外: 警告信息：
In as.matrix.csr(x, nrow = nrow, ncol = ncol, eps = eps) :
ncol*nrow indivisable by length(x)

使用道具举报

ywh19860616 发表于 2013-3-15 15:15:11 |显示全部楼层 |坛友微信交流群

xiaomuyoucun 发表于 2013-3-15 15:11
你好，我发现这一段程序总是运行不过去，你之前运行过没有错误吗？为什么我的总出现以下错误？
急于求解 ...

你的最后目的是什么？运行面板分位回归模型吗？
如果答案是，那你直接利用rqpd包就可以了，作者
已经把这些程序打包为package了，可以直接用

一份耕耘，一份收获。

使用道具举报

ywh19860616 发表于 2013-3-15 15:15:53 |显示全部楼层 |坛友微信交流群

xiaomuyoucun 发表于 2013-3-15 15:11
你好，我发现这一段程序总是运行不过去，你之前运行过没有错误吗？为什么我的总出现以下错误？
急于求解 ...

你是为了运行面板分位回归吗？
如果是，那用rqpd包就可以了

一份耕耘，一份收获。

使用道具举报

xiaomuyoucun

发表于 2013-3-15 15:18:53 |显示全部楼层 |坛友微信交流群

ywh19860616 发表于 2013-3-15 15:15
你是为了运行面板分位回归吗？
如果是，那用rqpd包就可以了

我不是很清楚这个包rqpd怎么用。。主要是我还实现惩罚项和工具变量

使用道具举报

xiaomuyoucun

发表于 2013-3-15 15:20:00 |显示全部楼层 |坛友微信交流群

ywh19860616 发表于 2013-3-15 15:15
你是为了运行面板分位回归吗？
如果是，那用rqpd包就可以了

还有我不清楚为什么我运行那一段程序时总汇报那个错误，你当时没遇到吗？怎么解决的？

使用道具举报

ywh19860616 发表于 2013-3-15 15:30:36 |显示全部楼层 |坛友微信交流群

xiaomuyoucun 发表于 2013-3-15 15:18
我不是很清楚这个包rqpd怎么用。。主要是我还实现惩罚项和工具变量

library(SparseM)
library(quantreg)
library(Ecdat)
data("Produc")
x1=log(Produc$pcap)
x2=log(Produc$pc)
x3=log(Produc$emp)
x4=Produc$unemp
X=cbind(x1,x2,x3,x4)
y=log(Produc$gsp)
N<-length(y)
s= rep(1:48,rep(17,48))
X <- as.matrix(X)
p <- ncol(X)
n <- length(levels(as.factor(s)))
Z <- as.matrix.csr(model.matrix(~as.factor(s)-1))
lambda=1
w=c(0.1,0.25,0.5,0.25,0.1)
taus=c(0.1,0.25,0.5,0.75,0.9)
K <- length(w)
xx1<- cbind(as(w,"matrix.diag.csr") %x% X,as.matrix(w) %x% Z)
Penalty <- cbind(as.matrix.csr(0,n,K*p),lambda*as(n,"matrix.diag.csr"))
xx<-rbind(xx1,Penalty)
y <- c(w %x% y,rep(0,n))
a <- c((w*(1-taus)) %x% (t(X)%*%rep(1,N)),sum(w*(1-taus)) * (t(Z) %*% rep(1,N)) + lambda * rep(1/2,n))
rq.fit.sfn(xx,y,a)