求用Gauss做非参数回归的方法，不胜感激

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

x=[0.1875 0.3125 0.4375 0.625 0.875 1.125 1.375 1.625 1.875 2.25 2.75 3.25];
y=[24124 499 2221 648 97 31 7 2.5 1.7  0.388  0.113 0.026];
对x,y做线性拟合时用的公式为y=a*(x^b)*exp(-c*x);
如果我想进行非参数化拟合x,y。请问用Gauss怎么做？
不胜感激！

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：非参数回归 GAUSS 不胜感激 非参数 USS 不胜感激

帮定，等高人！

坐等高人啊

再次求助，急急急

简单啊 两边取ln 原来的等式变为：
y=ln(a) + b*ln(x) - c*x;
然后做线性回归就好了

@ Lession 7.6: Nadaraya-Watson Nonparametric Regression Estimation @
new;cls;library pgraph;
n=100;
x=seqa(0,(2*pi)/(n-1),n);
mx=sin(x);
y=mx+0.5*rndn(n,1);
othband=1.06*stdc(x)*(n^(-1/5)); @最优光滑参数@
xy(x,mx~y);
print " 显示sin(x)和y图形 press any key to continue...";
waitc;cls;
{px,py}=npreg(y,x,1,othband);
print " 显示sin(x)、y和NW核估计py图形 ";
xy(x,mx~py~y);

proc(2)=npreg(y,x,krn,h); @ Kernel Regression @
local strt,endd,kern,pts,px,py,i,t,yt;
    strt=minc(x); endd=maxc(x); pts=100;
px=seqa(strt,(endd-strt)/(pts-1),pts);
py=px;
    i=1; do while i<=pts;
      if krn==1; t=nkernel(py[i,1],x,h);
      else; t=ekernel(py[i,1],x,h);
      endif;
      yt=t.*y; py[i,1]=sumc(yt)/sumc(t);
    i=i+1; endo;
retp(px,py);
endp;

proc nkernel(x,xi,h); @ Gaussian kernel @
local u,kern;
   u=(x-xi)./h;kern=(1/sqrt(2*pi))*exp(-(1/2).*(u^2));
retp(kern);
endp;

proc ekernel(x,xi,h); @ Epanechnikov kernel @
local u,kern;
  u=(x-xi)./h;kern=(0.75*(1-u^2)).*(abs(u).<=1);
retp(kern);
endp;

xuelida 发表于 2012-12-16 13:38
@ Lession 7.6: Nadaraya-Watson Nonparametric Regression Estimation @
new;cls;library pgraph;
n=100 ...

您好，请问能请您来我的帖子下回答一个非参数相关的问题吗？感激不尽！

同求同求