SAS中文论坛

 找回密码
 立即注册

扫一扫,访问微社区

查看: 628|回复: 2
打印 上一主题 下一主题

求助:如何对预测值与实验值进行拟合度的检验、相关分析?

[复制链接]

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
楼主
 楼主| 发表于 2009-5-16 18:20:27 | 只看该作者

求助:如何对预测值与实验值进行拟合度的检验、相关分析?

各位大侠:
   小女子通过对不同实验条件进行优化设计后,根据实验测定值进行响应面(RSREG程序)优化分析后得到一个预测方程。现在需对各实验条件下的预测值、实验值进行分析以检测模型的可信度,并以相关系数来表示。
   那么,请问:如何对不同条件下的预测值、实验值进行相关分析、拟合度检验呢?因为本人不是统计学专业,而只是平时喜欢、也习惯了用SAS的固定程序粘贴过来后进行统计分析,所以本人只能笨笨的把实验数据列出来,然后请求大侠列出具体的程序语言、或命令语式及其用法格式。

  若能得到大侠的帮助,本人将不胜感激!同时也感谢论坛为此提供了这样一个好的平台!

  1-32号实验条件下的预测值、实测值数据分别列如下:

4.64         4.35
6.66         7.02
7.03         7.62
7.13         7.34
5.33         5.77
6.97         7.60
5.46         5.80
5.94         6.05
5.28         6.34
6.92         7.20
7.29         6.75
7.77         8.69
5.60         5.60
7.62         7.89
6.11         6.68
6.21         6.27
6.71         6.43
7.62         7.21
2.43         2.28
1.99         1.50
9.13         8.47
10.11         9.70
4.14         3.78
6.27         5.57
8.05         8.43
8.43         6.99
7.16         7.87
7.16         8.24
7.16         7.24
7.16         6.62
7.16         7.15
7.16         6.92
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
沙发
 楼主| 发表于 2009-5-19 21:16:06 | 只看该作者

Re: 求助:如何对预测值与实验值进行拟合度的检验、相关分析?

1.相关分析:很简单,用proc corr过程。
2.拟合有度检验:
方法一:用proc npar1way过程,最好对EDF进行检验,因为这样才能够看出其分布特征是否一致。当然,你也可以检验中位数。
这两个过程在SAS中都是非常简单的,需要提及一点,你的数据格式可不是你列出来的,需要增加一列flag.如原始值标记为1,预测值标记为0。 然后把你提问时所列的两列数据合成一列如newcol.
程序如下:
proc npar1way edf  data=nkwilling;
      class flag;
      var newcol;
run;

方法二:用配对检验,用proc univariate过程,此时数据格式如你所写,但是也要增加一列deviance=原始值-预测值。程序如下:
proc univariate data=nkwilling loccount;
var deviance;
run;
建议采用第二种方法。因为你的需求来看,似乎配对检验更合适。注意如果使用方法二,在结果中一定要看符号秩检验(signed rank)所对应的p-value.
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
板凳
 楼主| 发表于 2009-5-20 09:07:00 | 只看该作者

Re: 求助:如何对预测值与实验值进行拟合度的检验、相关分析?

结果都在proc rereg的output输出结果里面,从你的描述看来,你是看不懂结果,而不是需要再行进行检验分析。
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|SAS中文论坛  

GMT+8, 2026-2-5 17:39 , Processed in 0.073025 second(s), 20 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表