rstool命令的用处??
上周学习了线性回归,今天上机。牙膏销量与价格、其他公司价格、广告费的关系。
X1=[-0.05 0.25 0.60 0 0.25 0.20 0.15 0.05 -0.15 0.15 0.20 0.10 0.40 0.45 0.35 0.30 0.50 0.50 0.40 -0.05 -0.05 -0.10 0.20 0.10 0.50 0.60 -0.05 0 0.05 0.55]
X2=[5.50 6.75 7.25 5.50 7.00 6.50 6.75 5.25 5.25 6.00 6.50 6.25 7.00 6.90 6.80 6.80 7.10 7.00 6.80 6.50 6.25 6.00 6.50 7.00 6.80 6.80 6.50 5.75 5.80 6.80]
X3=X2'.*X2'
X4=X1'.*X2'
X=[ones(30,1) X1' X2' X3 X4];
Y=[7.38 8.51 9.52 7.50 9.33 8.28 8.75 7.87 7.10 8.00 7.89 8.15 9.10 8.86 8.90 8.87 9.26 9.00 8.75 7.95 7.65 7.27 8.00 8.50 8.75 9.21 8.27 7.67 7.93 9.26]
[b,bint,r,rint,stats]=regress(Y',X);
b,bint,r,rint,stats;
pause (1)
rcoplot(r,rint)
X5=[X1',X2']
rstool(X5,Y,'interaction')
最后一个rstool 是次响应曲面模型函数。'interaction':表示包括常数项、一次项和交叉项。
但是这个函数到底是干什么用?那个图改变X1,X2是可以预测么??