|
r*****t 发帖数: 286 | 2 ☆─────────────────────────────────────☆
whistleruk (whistleruk) 于 (Fri Feb 9 15:02:53 2007) 提到:
1。给出uniform()在每次调用gasdev()过程中被调用次数的估计。
double uniform(long *seed);
/* uniform() returns a uniform random variable in (0,1) */
double gasdev(long *seed)
/* Returns a normal random variable */
{
static int iset = 0;
static double gset;
double fac, rsq, v1, v2;
if (iset == 0) {
do {
v1 = 2.0 * uniform(seed) - 1.0;
v2 = 2.0 * uniform(seed) - 1.0;
rsq = v1 * v1 + v2 * v2;
} while (rsq >= 1.0 || rsq = |
|
G******2 发帖数: 579 | 3 他们没有report p value 和t statistics .我还在想是我不懂的原因,一般来说我先
看p value. 不过我可以要求report
R square 看起来是model 的。
给我的sad code 打不丶开。我家电脑没有sas.
说是用的Bayesian schawz criteria selection. 跟我平常用的stepwise 和backwards
selection 不一样。不懂呀。
最小的SBC 和最大的RSQ的variable 被淘汰了,说intuitive judgement 不make
sense 是个dummy variable
用了两个variable 说 SBC 最小值,然后rsq 也最大。然后multivaraite 用了这两个
variable 一起,beta 都不大,没有report p value。f statistic 也没有给我。
我真希望当初读书时多念点Stats 。这都是phd 写的report。我看得稀里糊涂的。可是
平常看文献这些p value是一定要report的。
我是不是应该问这些phd要求数据支持他们选择... 阅读全帖 |
|
G******2 发帖数: 579 | 4 【 以下文字转载自 Quant 讨论区 】
发信人: Gigi2012 (Gigi), 信区: Quant
标 题: Re: 请问multi variate linear regression 选择risk factor 问题 (
发信站: BBS 未名空间站 (Fri Apr 10 01:32:27 2015, 美东)
他们没有report p value 和t statistics .我还在想是我不懂的原因,一般来说我先
看p value. 不过我可以要求report
R square 看起来是model 的。
给我的sad code 打不丶开。我家电脑没有sas.
说是用的Bayesian schawz criteria selection. 跟我平常用的stepwise 和backwards
selection 不一样。不懂呀。
最小的SBC 和最大的RSQ的variable 被淘汰了,说intuitive judgement 不make
sense 是个dummy variable
用了两个variable 说 SBC 最小值,然后rsq 也最大。然后multivara... 阅读全帖 |
|
f*******i 发帖数: 8492 | 5 我用的是SAS 9.2,参照网上的例子,使用同样的命令行来运行不同的数据,但是输出结果的项目名称
和数量却完全不一样。
人家的输出结果,在“Centroid Hierarchical Cluster Analysis”下有如下几个项目
Variable,Mean,Std,Dev,Skewness,Kurtosis,Bimodality。
而我的输出结果,却只有如下四个项目Eigenvalue,Difference,Proportion,Cumulative。
在“Cluster History”, 人家有NCL,Clusters Joined,FREQ,STD,SPRSQ,RSQ和Dist。
而我的输出结果却只有NCL, cluster joined,FREQ, Norm cent dist,tie。
请问,这是怎么回事啊? 命令明明是一样的。
|
|
n*****t 发帖数: 41 | 6 谢谢,文件如下。我只要读取最后一个数,就是那个likelihood的数值。
MERLIN 1.1.2 - (c) 2000-2007 Goncalo Abecasis
References for this version of Merlin:
Abecasis et al (2002) Nat Gen 30:97-101 [original citation]
Fingerlin et al (2004) AJHG 74:432-43 [case selection for
association studies]
Abecasis and Wigginton (2005) AJHG 77:754-67 [ld modeling, parametric
analyses]
Fingerlin et al (2006) Gen Epidemiol 30:384-96 [sex-specific maps]
Chen and Abecasis (2007) AJHG 81:913-26 [qtl assoc... 阅读全帖 |
|
h***i 发帖数: 3844 | 7 这个Rsq其实不靠谱啊。有些点weight太高了。 |
|