Summer School STATS108 & 208 官方问答贴(恭喜得到A+的同学们)

z-score

我帮你，你还bs我。我没事撑得？

^曼陀罗^

我错了~~呵呵。。。。。指条明路给我吧~~~谢谢~

^曼陀罗^

你QQ都不上了~~看来我以后最少一天要来顶一次贴了~~~

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:01 发表
哈哈。。。先在这里小小BS下某人。。。。

好嘛。。重复下问题。。。

我做第三题的时候.....用> resid.plot(body.fit)来画....那个...p value=0.0567...那是不是不需要再TAKE LOG 或者QUADRATIC...

得到这个结果的确难决定要不要用quadratic或者log transformation，因为P-value是some to weak evidence。你可以先fit一个quadratic然后看一下summary(body.fit1)的quadratic term P-value是不是significant. 如果你要take log的话，从resid.plot(body.fit)的图是看不出来的。唯一用log transformation的时候是在你的residuals有non-constant scatter.

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:05 发表
你QQ都不上了~~看来我以后最少一天要来顶一次贴了~~~

不好意思，我这两天是没怎么上。你觉得QQ上问题方便是吗？

在这里问不是很好是吗？

z-score

咳，别提QQ了，就是MSN我都潜水很久了.

^曼陀罗^

第三题。。。。。。。12个VARIABLES。。。一个个试过去。。偶头都会白了呢~呵呵~~~

在我看来。。。是CONSTANT SCATTER咯。。。- -！

QQ的好处就是可以截图啊。。。这里。。发不了图给你看~~~

^曼陀罗^

去STATS 20X的FORUM看了下~~~~

Stats!  Feb 5 2008, 11:02 AM Post #1  
For question three i got weak evidence over necessity for a quadratic term, however it is made so difficult by the fact there are nearly ten variables. Should I attempt to fit a quadratic?  
   

mike  Feb 5 2008, 11:23 AM Post #2  
No.  
   

0000  Feb 5 2008, 03:55 PM Post #3  
would it be any harm to put a quadrtic term on the model if it improves it? I have tried to put a quadrtic term to the model and made residual plot of the fitted values look better than it did.
(after the model building process, the p-value of no quadratic term needed is significant at the 5% level, so I thoght it was a good idea to fit a quadrtic term on the model.)  
   

mike  Yesterday, 09:44 AM Post #4  
The residual plot shows no real need for a quadratic term.  
   
guest  Yesterday, 07:23 PM Post #5  
but its way better if you log(height) ?

QUOTE(mike @ Feb 7 2008, 09:44 AM)
The residual plot shows no real need for a quadratic term.

mike  Today, 10:07 AM Post #6  
residual plot from the initial model shows constant scatter, so why would you log the response????????

Remember, a good model is a simple model.  
这么看的话。。就是两个都不需要做咯？？

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:13 发表
第三题。。。。。。。12个VARIABLES。。。一个个试过去。。偶头都会白了呢~呵呵~~~

在我看来。。。是CONSTANT SCATTER咯。。。- -！

QQ的好处就是可以截图啊。。。这里。。发不了图给你看 ...

12个variables不用一个一个fit, 那不得累死。你看最开始的pairs(body.df) pairs plot. 你看第一行，first row，哪个explanatory variable与response有curvi-linear relationship就先fit 他的quadratic就好。

我发帖子想帮助大家读stats的朋友们。你发到这里可以让其他的人都看到，如果有图片你也可以插图是不？

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:16 发表
去STATS 20X的FORUM看了下~~~~

Stats!  Feb 5 2008, 11:02 AM Post #1  
For question three i got weak evidence over necessity for a quadratic term, however it is made so difficult by the fact the ...

对，我和Mike的看法一样。你不妨试一下。

^曼陀罗^

这个老师也真好。。这不是告诉我们不需要LOG也不需要QUADRATIC么。。。

我最开始写的> pairs.20x(data.frame(height,chest,waist,thigh,knee,calf,ankle, bicep,wrist,age,weight,gender))是把全部的VARIABLES都写进去。。。。。这个是PAST ASSIGNMENT的格式。

但是CB上面好像只有用到quantitative VARIABLE。。。

怎么办。。两个做法会得到一样的答案么

^曼陀罗^

可以插图的吗。。。不会。。。。

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:27 发表
这个老师也真好。。这不是告诉我们不需要LOG也不需要QUADRATIC么。。。

我最开始写的> pairs.20x(data.frame(height,chest,waist,thigh,knee,calf,ankle, bicep,wrist,age,weight,gender))是把全部的VA ...

不用那个data.frame什么的，直接pairs.20x(body.df)就好了。

pairs plot不能显示qualitative variables.

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:28 发表
可以插图的吗。。。不会。。。。

这个帖子教你怎样插图。
http://www.skykiwi.com/bbs/frame.php?frameon=yes&referer=http%3A//www.skykiwi.com/bbs/viewthread.php%3Ftid%3D648505%26extra%3Dpage%253D1

^曼陀罗^

pairs.20x(body.df)和我  pairs.20x(data.frame(height,chest,waist,thigh,knee,calf,ankle, bicep,wrist,age,weight,gender))得到的圖一樣..- -!老師爲什麽不教這種簡單的......

還有....然後看summary(body.fit)...其中有幾個很大的p value...我們是不是選最大的那個..然後去掉它..從新弄一個fit1..然後再看哪個p-value最大..再去掉....一直到
Coefficients:
             Estimate Std. Error t value Pr(>|t|)   
(Intercept) 213.32496    6.45815  33.032  < 2e-16 ***
chest        -0.15178    0.06884  -2.205   0.0280 *  
waist        -0.66706    0.05725 -11.652  < 2e-16 ***
thigh        -0.58953    0.09730  -6.059 3.07e-09 ***
calf         -0.62432    0.13233  -4.718 3.26e-06 ***
bicep        -1.02794    0.14782  -6.954 1.39e-11 ***
wrist         0.63009    0.36441   1.729   0.0845 .  
weight        1.42233    0.06371  22.324  < 2e-16 ***
gender        6.22182    0.95690   6.502 2.28e-10 ***

wrist就不需要再去掉..然後重新弄個fit4了吧?

BuyQueen

阿童木的帖子，顶一个

z-score

原帖由 ^曼陀罗^ 于 2008-2-8 22:53 发表
pairs.20x(body.df)和我  pairs.20x(data.frame(height,chest,waist,thigh,knee,calf,ankle, bicep,wrist,age,weight,gender))得到的圖一樣..- -!老師爲什麽不教這種簡單的......

還有....然後看summar ...

恩，阿童木知道很多老师不告诉你的窍门。

恩，对，在drop掉wrist。

^曼陀罗^

为什么连WRIST都要DROP掉呢？那不是有WEAK EVIDENCE AGAINST THE HYPOTHESIS THAT GENDER IS NOT NEDDED IN THE MODEL吗？因为P-VALUE=0.0845..

CB上的CASE STUDY33：TECHNITRON是这样的呀。。。。

t-multiplier

原帖由 ^曼陀罗^ 于 2008-2-9 09:20 发表
为什么连WRIST都要DROP掉呢？那不是有WEAK EVIDENCE AGAINST THE HYPOTHESIS THAT GENDER IS NOT NEDDED IN THE MODEL吗？因为P-VALUE=0.0845..

CB上的CASE STUDY33：TECHNITRON是这样的呀。。。。 ...

The famous Professor Alan Lee once said that "there is no correct model in regression model building, but some models will be more useful  than others". If you want to copy case studies than you can do so.

t-multiplier

If I was building the model I will drop wrist as I believe keeping the model simple is more important.

^曼陀罗^

沒有正確答案..- -!MARKING的时候不是有个标准答案的么?和老师的不一样也没事?

我drop掉試試看...怎麽樣決定是留着好..还是删掉好呢?

t-multiplier

原帖由 ^曼陀罗^ 于 2008-2-9 09:57 发表
沒有正確答案..- -!MARKING的时候不是有个标准答案的么?和老师的不一样也没事?

我drop掉試試看...怎麽樣決定是留着好..还是删掉好呢?

If people only want the model answers from me, I guess there is not much significance for my existence then.

Please refer to your assistance room tutor from now on.

晨曦.com

sammie????????

t-multiplier

原帖由 晨曦.com 于 2008-2-9 12:22 发表
sammie????????

Sammie is a lab demonstrator and the official 301 tutor.

^曼陀罗^

我不是要答案啦..我的意思是...你說.there is no correct model in regression model building..但是..改作業的時候不是都有model answer的麽... 你想歪了....

如果一個是有drop一個是沒有drop wrist的..那麽p-value還有 r square的值不是就會是不一樣的麽...
那樣有影響的嗎...

t-multiplier

原帖由 ^曼陀罗^ 于 2008-2-9 12:29 发表
我不是要答案啦..我的意思是...你說.there is no correct model in regression model building..但是..改作業的時候不是都有model answer的麽... 你想歪了....

如果一個是有drop一個是沒有dro ...

刚才我看了一下cast study Technitron, 他的那个gender有weak evidence没有drop是因为它是qualitative variable, qualitative variable需要看anova output P-value不是regression output P-value. 所有我认为wrist一定要drop掉，你看pairs plot wrist和bicep多么的correlate.

^曼陀罗^

不好意思...在問一下....."pairs plot wrist和bicep多么的correlate."這個...和drop wrist有什麽直接關係嗎...我不懂得怎麽看圖來選擇要不要drop~~

t-multiplier

原帖由 ^曼陀罗^ 于 2008-2-9 12:46 发表
不好意思...在問一下....."pairs plot wrist和bicep多么的correlate."這個...和drop wrist有什麽直接關係嗎...我不懂得怎麽看圖來選擇要不要drop~~

Pairs plot里面wrist和bicep的correlation是那个数字，好像是0.85吧，就是matrix main diagonal的对称数字。如果两个explanatory variables的correlation很高，这就代表这两个variables有multi-collinearity. Drop掉一个和保存两个没有什么区别，所以就drop因为model会更simple.

晨曦.com

請教一下,怎么做 time series的multi-regression?

^曼陀罗^

又学到新的知识了~哈哈...
再问一下...Q2..我画FITTED LINES ON THE ORIGINAL SCATTER PLOT...
> plot(Sales~Advertising,type="n",main="Fitted Model")
> text(Sales,Advertising,Promotions)
> abline(sales.fit$coef[1],sales.fit$coef[2])
> abline(sales.fit$coef[1],sales.fit$coef[3],sales.fit$coef[2])
> abline(sales.fit$coef[1],sales.fit$coef[4],sales.fit$coef[2])
我是哪里出错了....为什么画出来的只有两条线...而且不是平行的~~~~