| 5 | 1/1 | 返回列表 |
| 查看: 3294 | 回復(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:我是做分析方法的。 |
榮譽版主 (著名寫手)
![]() |
專家經(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. |

榮譽版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |

新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +7 | c297914 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|