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

8494

收藏 2015-12-04

请教各位大神，针对非参和半参模型，如何通过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.

先多谢啦!

扫码加我拉你入群

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

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

全部回复

victorchan0633

2015-12-4 16:27:20

第一种情况：
> 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
要具体问题具体分析，非常复杂，也是先随机模拟出一些数据，然后半参数非参数窗宽设置核函数等不同设置，方法应该会随着数据不同而不同。

扫码加我拉你入群

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

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

vapor__9

2015-12-7 07:54:36

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)}

好的，我试试看，谢谢哇

扫码加我拉你入群

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

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

vapor__9

2015-12-7 07:55:23

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)}

好的，我试试看，谢谢哇

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群