求教：非参和半参模型如何用R得到MSE (mean squared error)?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请教各位大神， 针对非参和半参模型，如何通过R得到mean squared error?
希望能写出code （不用现成的package）。

举个例子，假设我的模型是: y=sin(u)*x1+4*x2+e, 其中u 是uniform distribution, x1,x2 是变量.现在需要求两种情况下的MSE：1，假设我们知道这个具体的模型的形式来求出MSE。2， 假设我们不知道这个半参的具体形式，即通过y=a(u)*x+e （这里的X其实分别是x1 x2, 所以a（u）也有两个，分别是x1 和x2的“系数”）也就是说不知道哪个系数是常数，哪个是function， 的非具体形式求出MSE.

先多谢啦!

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Squared Square Error mean red package error 模型 如何

第一种情况：
> rm(list = ls())
> a<-seq(0, 1, length.out = 11)
> b<-seq(0, 1, length.out = 11)/11
> d<-seq(2,3,length.out = 11)
> y<-c()
> for(i in 1 :11){y[i]=sin(b[i])*a[i]+4*d[i]+rnorm(1)}
> mse<-c()
> for(i in 1:11){mse[i]<-(y[i]-(sin(b[i])*a[i]+4*d[i]))}
> tem<-(sum(mse))^2+mse%*%mse/11
要具体问题具体分析，非常复杂，也是先随机模拟出一些数据，然后半参数非参数窗宽设置核函数等不同设置，方法应该会随着数据不同而不同。

victorchan0633 发表于 2015-12-4 16:27
第一种情况：
> rm(list = ls())
> a b d y for(i in 1 :11){y=sin(b)*a+4*d+rnorm(1)}

好的，我试试看，谢谢哇

victorchan0633 发表于 2015-12-4 16:27
第一种情况：
> rm(list = ls())
> a b d y for(i in 1 :11){y=sin(b)*a+4*d+rnorm(1)}

好的，我试试看，谢谢哇