全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Gauss专版
3747 7
2012-12-03
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怎么做?
不胜感激!
二维码

扫码加我 拉你入群

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

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

全部回复
2012-12-3 20:25:46
帮定,等高人!
二维码

扫码加我 拉你入群

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

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

2012-12-3 21:20:06
坐等高人啊
二维码

扫码加我 拉你入群

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

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

2012-12-5 12:46:33
再次求助,急急急
二维码

扫码加我 拉你入群

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

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

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

二维码

扫码加我 拉你入群

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

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

2012-12-16 13:38:43
@ 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;
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群