s******s
发帖数: 145 | 1 有两组数，两个mean,只知道下面四个值
N1=1000, mean1=12
N2=1010, mean2=10.5
用哪个test可以来看看mean1和mean2是不是significant different? 还是说条件不足
， 需要知道StdDev什么的。
晕了
第一个给出正确意见的发5个包子 |
|
y***n
发帖数: 51 | 2 我来说说吧。这应该算longitudinal data了，最简单的就是比较他们的mean. 相信你
数据有一列是时间吧。 你可以用proc means 分别算mobility 和 memory 的均值，
where time=week0 (baseline), time=month1 （只有第一个月的） or month12. 然
后对两个treatment比较。
你要是一起算呢， 可以先把数据对treat time排序，然后
proc means data=yourdata;
var mobility; *memory;
class treat time;
output out=mean1;
run;
然后用mean1里面的数据可以画出每个treatment的均值随时间变化的曲线。
如果两个treatment的均值在baseline相差很大，你可以考虑算mobility （time）-mobility （
baseline）在每个时间点的均值，同样对memory。
如果要模型，你可以考虑用GEE， Genmod。mobility和memory的correlation可以在模
型里面... 阅读全帖 |
|
z****e
发帖数: 54598 | 3 计算一个数组的mean和variance
由于这是两个normal相加，所以mean和variance都是两个独立的mean和variance的关系
翻统计书，我不太记得了，不过应该是如下关系
mean = (mean1 + mean2)/2
variance = variance1 + variance2 + 2Cov12
从题目看，这个Cov12应该是0
然后就是解二元二次方程组
java写sqrt什么要用Math.sqrt来做 |
|
y***j
发帖数: 11235 | 4 基本概率问题，我这种土人都能回答,p=0.05时候，１００随机表达的基因，中间平均
会有５个显著。
ｐｖａｌｕｅ和ｆｏｌｄ ｃｈａｎｇｅ完全两个不同的量度。
ｐｖａｌｕｅ只说明这两个东西是不是来自一个群体。如果ｐ很小说明这两个东西很大
概率上是不同的，具体有多大的不同，是不能确定的。例如你方差很小的时候mean 1和
mean 1.1会有区别，但是有可能1和1.1的区别在统计上显著但是生物学上没有意义。
但是方差大了的时候mean1 和mean10可能统计上不显著。这里说明这两个东西可能来自
一个群体，如果你不相信，就要算一下power，如果power很小的话很可能是你的实验重
复不够，加些实验次数解决了。
fold change上的区别对与生物体来说可能更重要，但是这个fold change是真正的差别
还是随机误差呢？p value是回答的这个问题。
另外仔细想想fold change也可能是无意义的，尤其是在基因组水平的研究。例如，一
个基因可能只需要2-4个特异的转录因子，group A 平均表达100个，group B平均表达
1000,其实对生物的影响也是一样的。
老夫经常... 阅读全帖 |
|
d*****n
发帖数: 65 | 5 proc summary data=one;
var x y z;
output out=temp1(drop=_type_ _freq_) mean= ;
proc transpose data=temp1 out=temp2 prefix=mean;
proc summary data=temp2;
var mean1;
output out=data(drop=_type_ _freq_) mean= ;
run; |
|
w******n
发帖数: 309 | 6 I have 100 files(from the simulation results), say dat1 ... dat100, for each
file, I would like to some analysis, for example, get the mean of each file
. Then I would like to save the analysis results into different files: say
mean1.txt .... mean100.txt. I know how to do for each file, however, I don't
know how to do the loop from file 1 to fil1 100.
thanks a lot1 |
|
o****o
发帖数: 8077 | 7 你这个用BY statement在SAS里面也不难，分条件运行不同回归也很容易。比如
data original;
array _x{*} var1-var10;
do i=1 to 10000;
do _j=1 to dim(_x); _x[_j]=rannor(7655)+sin(i + _j); end;
output;
drop i _j;
end;
run;
ods select none;
proc surveyselect data=original out=samp rep=10
sampsize=1000 method=srs;
run;
proc means data=samp;
by replicate;
var var1;
output out=_mean mean(var1)=mean1;
run;
data samp |
|
b*w
发帖数: 2062 | 8 多谢，那么如何估计两个（不完全独立，先当独立对待）truncated normal random
variable的和的mean/std呢？是不是也是mean = mean1+mean2, std^2=std1^2+std2^2？ |
|
x****n
发帖数: 43 | 9 sample1 <- rnorm(n1, mean=mean1, sd=sd1)
sample2 <- rnorm(n2, mean=mean2, sd=sd2)
t.test(sample1, sample2) |
|
k*******a
发帖数: 772 | 10 data a;
input value freq_1 freq_2;
datalines;
3 12 11
4 4 3
5 6 15
6 7 13
7 9 11
8 11 4
;
run;
proc sql;
select sum(value*freq_1)/sum(freq_1) as mean1, sum(value*freq_2)/sum(freq_2)
as mean2
from a;
quit; |
|
s***r
发帖数: 1121 | 11 How can I estimate t value using the heteroskedasticity-consistent standard
errors in a PROC MEAN (see below)? (NOT in a PROC REG)
proc means data=a1 noprint; var var1; by byvariable;
output out=a2 mean=mean1 t=tret;
run;
3 baozi will be given for the first 3 replies. thanks. |
|
