24小时热门版块排行榜    

查看: 2829  |  回复: 4
【悬赏金币】回答本帖问题,作者09yschen将赠送您 40 个金币

09yschen

金虫 (正式写手)

[求助] 50金求解:关于lack-of-fit test的意义及应用

我得到的标准曲线是用origin进行线性拟合的,但审稿人说加个lack-of-fit test,因为我的r值不太好。有位大神对我的数据用matlab进行了test,,请见http://muchong.com/bbs/viewthread.php?tid=4020335。但问题是,我看不懂……
我想请问:标准曲线是用origin进行线性拟合的还是用lack-of-fit test得到的呢,帖子里的P值能说明什么,应该如何应用,我看到直接给了个P值,但也有人用F值,真的不明白,想请教下,对下那位大神的数据,我应该如何应用。
求解……急
我特意取了一些文献的表达,请见图,我应该用哪个好。
PS:我是做分析方法的。






回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
09yschen(金币+10): 2012-01-08 20:51:17
见帖
http://en.wikipedia.org/wiki/Lack-of-fit_sum_of_squares
http://muchong.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.
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2012-01-06 22:37:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

lack of fit,失拟检验,评估模型的拟合度。如果 p 值小于您选择的 a 水平,则证明模型未能与数据准确拟合。您可能需要添加项,或者变换数据,以便更准确地为数据建模。
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
3楼2012-01-06 23:15:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

09yschen

金虫 (正式写手)

引用回帖:
3楼: Originally posted by dbb627 at 2012-01-06 23:15:49:
lack of fit,失拟检验,评估模型的拟合度。如果 p 值小于您选择的 a 水平,则证明模型未能与数据准确拟合。您可能需要添加项,或者变换数据,以便更准确地为数据建模。

能不能帮我写一个程序,谢谢
4楼2012-01-07 23:58:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaopanpan

新虫 (小有名气)

楼主 你好 我现在也遇到了这样的问题  审稿人意见:Linearity should be tested either by ANOVA or a lack of fit test.请问您现在会了吗  能不能教一下
5楼2014-11-22 10:32:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 09yschen 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 37/1850 2024-06-17 02:17 by xmclyn
[基金申请] 最新消息:2024国自然下载文件名变了 +9 dream200606 2024-06-16 12/600 2024-06-17 00:34 by kobe0107
[找工作] 江西双非一本和四川双一流高校如何选择? 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[硕博家园] 博士毕业高校和就业的相关问题 +5 SCITOPPP 2024-06-14 8/400 2024-06-16 21:33 by kingkongking
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +9 zju2000 2024-06-12 15/750 2024-06-16 20:45 by 北溟鱼1318
[考博] 34岁读博士晚吗 +27 emitdne 2024-06-13 27/1350 2024-06-16 18:57 by yugege2009
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +7 恐龙爸爸 2024-06-14 7/350 2024-06-16 14:48 by redfish105
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 8/400 2024-06-16 10:14 by 午睡未进行
[基金申请] 博后面上今天有bug可以看到是否资助? +20 lyfbangong 2024-06-12 31/1550 2024-06-15 21:18 by since—2010
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[基金申请] 博后特助这周出结果吗?往年都是啥时候啊? +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
信息提示
请填处理意见