| 5 | 1/1 | 返回列表 |
| 查看: 3293 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[求助]
50金求解:關(guān)于lack-of-fit test的意義及應(yīng)用
|
||||
|
我得到的標(biāo)準(zhǔn)曲線是用origin進(jìn)行線性擬合的,但審稿人說加個lack-of-fit test,因為我的r值不太好。有位大神對我的數(shù)據(jù)用matlab進(jìn)行了test,,請見http://www.gaoyang168.com/bbs/viewthread.php?tid=4020335。但問題是,我看不懂…… 我想請問:標(biāo)準(zhǔn)曲線是用origin進(jìn)行線性擬合的還是用lack-of-fit test得到的呢,帖子里的P值能說明什么,應(yīng)該如何應(yīng)用,我看到直接給了個P值,但也有人用F值,真的不明白,想請教下,對下那位大神的數(shù)據(jù),我應(yīng)該如何應(yīng)用。 求解……急 我特意取了一些文獻(xiàn)的表達(dá),請見圖,我應(yīng)該用哪個好。 PS:我是做分析方法的。 |
榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |

榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |
|
見帖 http://en.wikipedia.org/wiki/Lack-of-fit_sum_of_squares http://www.gaoyang168.com/bbs/viewthread.php?tid=4020335&pid=7&page=1#pid7 LOFRTEST Lack-of-fit test for regression model with independent replicate values. LOFRTEST(D,alpha) is a statistical test that gives information on the form of the model under consideration. A significant lack-of-fit suggest that there may be some systematic variation unaccounted for in the hypothesized model (chosen model does not well describe the data). It arises when there are exact replicate values of the independent variable in the model that provide an estimate of pure error. Pure error is in essence the amount of error that cannot be accounted for by any model. Then allows a test on whether there is error present aside from pure error. For the construction of the lack-of-fit test we need to examine three common types of linear models: - single mean (one parameter) - slope and intercept or common regression model (two parameters) - separate means for each x-value or one-way ANOVA (many parameters). So, the pure error is the error of the separate means on ANOVA and the total error in the residual resulting in the regression analysis: the lack-of-fit results to be the difference between this two sources of error, SS(LOF) = SSR(Model) - SSE(ANOVA). Syntax: lofrtest(D,alpha) Inputs: D - matrix data (=[X Y]) (last column must be the Y-dependent variable). (X-independent variable entry can be for a simple [X], multiple [X1,X2,X3,...Xp] or polynomial [X,X^2,X^3,...,X^p] regression model). alpha - significance level (default = 0.05). Outputs: A complete summary (table) of analysis of variance partitioning sources of variation for testing lack-of-fit. Example from the data on height and weight of 19 students in Psy 202. Assigment 3 of Psych 3030 from the Department of Psychology of the York University. Available on Internet at the URL address http://www.psych.yorku.ca/lab/psy3030/assign/assign3.htm We are interested to test with a significance-value = 0.05 if there is a lack-of-fit on the regression model due to the height replicate values. ------------------- ------------------- Height Weight Height Weight ------------------- ------------------- 60 90 68 140 60 100 68 135 62 110 70 160 62 116 70 145 62 120 70 148 66 140 71 143 66 170 71 135 68 130 74 195 68 117 74 164 68 155 ------------------- ------------------- Data matrix must be: D=[60 90;60 100;62 110;62 116;62 120;66 140;66 170;68 130;68 117;68 155;68 140;68 135; 70 160;70 145;70 148;71 143;71 135;74 195;74 164]; Calling on Matlab the function: lofrtest(D) Answer is: Lack-of-fit test for regression model with independent replicate values. -------------------------------------------------------------------------- SOV SS df MS F P -------------------------------------------------------------------------- Model 7658.359 1 7658.359 31.720 0.0000 Residual 4104.378 17 241.434 -------------------------------------------------------------------------- Lack-of-fit 2142.011 5 428.402 2.620 0.0797 Pure error 1962.367 12 163.531 -------------------------------------------------------------------------- Total 11762.737 18 -------------------------------------------------------------------------- If the associated P-value for any F test is equal or larger than 0.05 The corresponding null hypothesis is met. Otherwise it is not met. Created by A. Trujillo-Ortiz, R. Hernandez-Walls, A. Castro-Perez and F.J. Marquez-Rocha Facultad de Ciencias Marinas Universidad Autonoma de Baja California Apdo. Postal 453 Ensenada, Baja California Mexico. atrujo@uabc.mx Copyright (C) March 4, 2005. To cite this file, this would be an appropriate format: Trujillo-Ortiz, A., R. Hernandez-Walls, A. Castro-Perez and F.J Marquez-Rocha. (2005). lofrtest:Lack-of-fit test for regression model with independent replicate values. A MATLAB file. [WWW document]. URL http://www.mathworks.com/matlabcentral/ fileexchange/loadFile.do?objectId=7074 References: Department of Psychology of the York University. Available on Internet at the URL address http://www.psych.yorku.ca/lab/psy3030/assign/assign3.htm Zar, J. H. (1999), Biostatistical Analysis (2nd ed.). NJ: Prentice-Hall, Englewood Cliffs. p. 345-350. |

新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +7 | 谷雨上岸 2026-03-23 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|