矩阵运算

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> n<-100
> t<-runif(n,0,1)
> x1<-rnorm(n,5,1)
> x2<-rnorm(n,1,1)
> u1<-rnorm(n,0,1)
> u2<-rnorm(n,0,1)
> e<-rnorm(n,0,1)
> X1<-x1+u1
> X2<-x2+u2
> y1<-0
> for(i in 1:n)
+ {
+   y1[i]<-1+t[i]*t[i]+(X1[i]-u1[i])*sin(t[i])+(X2[i]-x2[i])*exp(t[i])+e[i];
+       }
> X<-X2
> Y<-t(y1)
> I<-diag(n)
> S<-t(X)%*%X
> P<-X%*%solve(S)%*%t(X)
> dn<-t(I)%*%P%*%I+n
> a<-solve(dn)%*%(t(I)%*%(I-P)%*%Y)
Error in solve.default(dn) :
  system is computationally singular: reciprocal condition number = 2.42432e-20
> b<-solve(S)%*%t(X)%*%Y-solve(S)%*%t(X)%*%I%*%a
Error in solve(S) %*% t(X) %*% Y : non-conformable arguments
>
a,b一直报错，求大神指点!!

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：矩阵运算 For system Error

你的dn有问题！dn无法求逆！
> dn<-t(I)%*%P%*%I+n
> solve(dn)
Error in solve.default(dn) :
  system is computationally singular: reciprocal condition number = 4.49055e-20
>

king64 发表于 2015-11-18 12:49
你的dn有问题！dn无法求逆！
> dn solve(dn)
Error in solve.default(dn) :

是的，dn的行列式值为0，不能求逆，可是变系数不就是这样求的？