SAS中文论坛

标题: t值计算 [打印本页]

作者: shiyiming 时间: 2010-10-27 20:01
标题: t值计算
data a: code Time day r _n
600001 20060524 20060531 -0.008281157 -5
600001 20060525 20060531 -0.026505991 -4
600001 20060526 20060531 0.065703898 -3
600001 20060529 20060531 0.086088239 -2
600001 20060530 20060531 0.03979086 -1
600001 20060531 20060531 0.107274286 0
600001 20060601 20060531 0.059430877 1
600001 20060602 20060531 -0.038942064 2
600001 20060605 20060531 -0.051635697 3
600001 20060606 20060531 0.038702317 4
600001 20060607 20060531 -0.062411551 5
600001 20060608 20060531 -0.017669674 6
600001 20060609 20060531 0.013122828 7
600001 20060612 20060531 0.016912312 8
600001 20060613 20060531 -0.051505329 9
600001 20060614 20060531 -0.038340433 10
600003 20070226 20070305 0.015148111 -5
600003 20070227 20070305 -0.056897604 -4
600003 20070228 20070305 0.034666183 -3
600003 20070301 20070305 0.112772922 -2
600003 20070302 20070305 -0.01225781 -1
600003 20070305 20070305 0.086679702 0
600003 20070306 20070305 -0.051920234 1
600003 20070307 20070305 -0.010350856 2
600003 20070308 20070305 -0.015317982 3
600003 20070309 20070305 -0.019223534 4
600003 20070312 20070305 -0.002773618 5
600003 20070313 20070305 0.014388151 6
600003 20070314 20070305 -0.008704821 7
600003 20070315 20070305 -0.002119721 8
600003 20070316 20070305 -0.029184168 9
600003 20070319 20070305 -0.022995407 10
600007 20060425 20060509 -5
600007 20060426 20060509 -4
600007 20060427 20060509 -3
600007 20060428 20060509 0.114413478 -2
600007 20060508 20060509 0.042151868 -1
600007 20060509 20060509 0.069974391 0
600007 20060510 20060509 0.018565172 1
600007 20060511 20060509 -0.008876887 2
600007 20060512 20060509 -0.008373075 3
600007 20060515 20060509 -0.050114431 4
600007 20060516 20060509 0.017978201 5
600007 20060517 20060509 0.026904472 6
600007 20060518 20060509 -0.015287463 7
600007 20060519 20060509 -0.030309676 8
600007 20060522 20060509 0.07694451 9
600007 20060523 20060509 -0.02235848 10
要求分别按照code 计算n= -5， -4，-3,..., 10天的均值和对应的t值，
另外需要把-2——0天的值加起来作为一个新的值求均值和t值，如何计算。
在利用如下程序时，计算出来的-2——0天的均值是分别计算-2，-1，0天的均值的1/3，当然手算也可以，但是不知道t值如何计算，请教大虾，谢谢
proc format;
value _n1Fmt -5=1 -4=2 -3=3 -2-0=4 1=5 2-3=6 4-10=9 11-20=10;
run;

作者: shiyiming 时间: 2010-10-28 10:07
标题: Re: t值计算
proc means里面输出的时候除了mean，再加上t就行了

作者: shiyiming 时间: 2010-10-28 10:19
标题: Re: t值计算
我知道，但是两种算法不一样，不知道怎么回事

作者: shiyiming 时间: 2010-10-28 14:47
标题: Re: t值计算
to vcvi
不明白你的问题是什么，你已经用proc format对n分了组，那么分组求mean和t就好了啊。

作者: shiyiming 时间: 2010-10-28 17:45
标题: Re: t值计算
proc format对n分了组，那么分组后 -2——0作为一个新的量要求 mean和t ，
我用单独的_n= -2算的均值例如是0.2， _n=-1的均值是0.3， _n=0的mean =0.2, 那么利用的proc format 分组后-2——0的均值是0.233，但我想要的是-2——0作为一个整体的均值是0.7，所以相应的t我想是有问题的. 程序如下
data A: code time day r _n
proc format;
value _n1Fmt -5=1 -4=2 -3=3 -2-0=4 ;
run;
proc means data = a noprint;
class code day _n; ways 3;
var r;
format _n _n1Fmt.;
output out = MeansOut mean = mean;
run;
proc means data = MeansOut noprint;
class _n; ways 1;
var mean;
output out = MeansOut1(drop = _type_ _freq_) n = n mean = mean t = t;
run;

作者: shiyiming 时间: 2010-10-28 17:51
标题: Re: t值计算
那就对每一code每一组别先求sum啊，然后在用proc mean求每一组的均值和t。

作者: shiyiming 时间: 2010-10-28 17:59
标题: Re: t值计算
to dengzi:
不会，在原来的程序上怎么改

作者: shiyiming 时间: 2010-10-29 02:56
标题: Re: t值计算
这个程序应该是这样的：
1.如果code day _n 三个变量一起决定了一个独特不重复的观测，
2.在求t的时候，每个code的-2-0观测值（3个）只用到平均值（一个值）。
如果你要的不是这个，那么这个程序是错的。
另外你给的数据是错误的。其中有三个缺失。
[quote:2h9cviae]600007 20060425 20060509 -5
600007 20060426 20060509 -4
600007 20060427 20060509 -3[/quote:2h9cviae]

作者: shiyiming 时间: 2010-10-29 03:06
标题: Re: t值计算
[quote:27mnzr25]我用单独的_n= -2算的均值例如是0.2， _n=-1的均值是0.3， _n=0的mean =0.2, 那么利用的proc format 分组后-2——0的均值是0.233，但我想要的是-2——0作为一个整体的均值是0.7，所以相应的t我想是有问题的.[/quote:27mnzr25]
明明是要sum，非要说是均值。
如果只这样改动如下

proc means data = a noprint;
class code day _n; ways 3;
var r;
format _n _n1Fmt.;
output out = MeansOut sum = mean;
run;
原因是如果只有一个值，sum =mean

作者: shiyiming 时间: 2010-10-29 21:28
标题: Re: t值计算
to jingju, thank you , 是我没表达清楚，如下就可以了？ way 3 ,way 1什么意思，可以解释下吗，谢谢。
proc means data = a noprint;
class code day _n; ways 3;
var r;
format _n _n1Fmt.;
output out = MeansOut [color=#FF0000:2m1htjp0]sum[/color:2m1htjp0] = mean;
run;
proc means data = MeansOut noprint;
class _n; ways 1;
var mean;
output out = MeansOut1(drop = _type_ _freq_) n = n mean = mean t = t;
run;

作者: shiyiming 时间: 2010-10-29 22:54
标题: Re: t值计算
[code:2i36zwiz]proc format;
value _n1Fmt -5='1' -4='2' -3='3' -2-0='4' 1='5' 2-3='6' 4-10='9' 11-20='10';
run;

proc means data =a nway noprint;
class code day _n;
var r;
format _n _n1Fmt2.;
output out = MeansOut sum = _v_;
run;
proc means data =MeansOut nway noprint;
class _n;
var _v_;
output out = MeansOut1(drop = _type_ _freq_) n =n mean =mean t =t;
run; [/code:2i36zwiz]
我也说不清楚，你查一下help。另外可以用nway代替如上。

作者: shiyiming 时间: 2010-11-1 13:17
标题: Re: t值计算
to jingju11
你这里的proc means 后面的nway是什么意思啊？

作者: shiyiming 时间: 2010-11-2 21:12
标题: Re: t值计算
nway就是输出最细分的组合们（highest level）。n也就是ways 的最高数字。三个class变量就等同于ways 3；
4个class 变量就是ways 4；等等

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)