l**********8 发帖数: 10 | 1 个人意见:先做一些transformation(取log、开方之类的),把分布整的像正态分布
一点,然后One-way ANOVA,如果显著的话再加pairwise comparison(Tukey test什么
的。。。) |
|
d********m 发帖数: 3662 | 2 用什么软件?主要几个就是tukey, bofferoni,就是加个命令行的问题。 |
|
j*****u 发帖数: 285 | 3 用的是graph pad,one-way anova 有让选择 post test, tukey, bofferoni 之类的,
但two-way anova就看不到有这个选项了 |
|
S*****h 发帖数: 21 | 4 不如这样说, 你一个图中只有2组数据的时候,用t test,3组或以上用one way anova.
其结果是在有 Corrections的情况下,告诉你组间数据有无明显差异。如果有差异,
此时你不知
道哪些组间有差异,所以你要辅以post hoc 分析。post hoc 又有很多种,最常见用
Bonferroni, tukey's之类。 |
|
u*h 发帖数: 397 | 5 在你的实验中, 有四组老鼠:
Group 1: WT-Saline
Group 2: WT-Drug
Group 3: TR-Saline
Group 4: TR-Drug
H0: u1=u2=u3=u4
H1: at least one inequality.
One-way ANOVA result: F value = 7.53, Pr = 0.0008
Reject H0, accept H1: at least the mean value of one group is different from
the others.
Use Tukey method to analysis the mean of the four groups:
u1 != u2: DRUG对WT老鼠有作用。
u1 != u4: TR-Drug 组的老鼠, 和WT-Saline的老鼠不同。
因为u1 != u4, 一个进一步的问题就是, u3和u1, u4 的关系是怎样的:
显然, 因为u1!=4, 那么 u3=u1, u3=u4不可能同时成立。 对此, 有三个假设:
1. u1 = u3, u3!... 阅读全帖 |
|
T***y 发帖数: 43 | 6 【 以下文字转载自 Linux 讨论区 】
发信人: Tukey (突击), 信区: Linux
标 题: 请教作科学计算哪种Linux比较好?
发信站: BBS 未名空间站 (Sun Sep 9 11:07:14 2007), 转信
主要是自己写C或者Fortran的程序
还有Matlab等科学计算软件
谢谢! |
|
l******s 发帖数: 82 | 7 在数学界,几个数学家虚构一个人作为他们共同的笔名并不是
Bourbaki的专利,1938年,princeton的数学家Boas,Ralph
Tukey,John和Smithies,Frank以Petard,Hector的笔名在
American mathematical monthly上发表了一篇题为“A con-
tribution to the mathematical theory of big-game hu-
nting”的搞笑文章,并用Petard为作者写了一份讲义,1939年
Boas在Cambridge大学遇见了weil,他们居然产生了要Petard
和Bourbaki的女儿Betti Bourbaki(当然也是虚构的)结婚的
想法,并用法语写了婚礼请柬发给当地的数学家
1950年,Boas应邀为《大不列颠百科全书》写了篇文章,其中
提到Bourbaki写的几卷书,并说Bourbaki其实是一群人的一个
笔名,Boas的这种出卖“岳父”的行为激怒了Bourbaki,不久
Boas收到一封匿名信,开头就说“你这个不知天高地厚的家伙,
你竟敢说Bourbaki不存在!”最后 |
|
n*********0 发帖数: 22 | 8 请问你说的multiplicity adjusted p-value是simultaneously testing, adjusted
family p-value吗?就是同时进行多项比较,控制family p-value小于一定值,而不是
控制单项p-value。
如果是的话,可以用sas软件计算。有三种procedures, Tukey, Scheffe,和
bonferroni。 具体的用法和code网上有,可以搜索到。 |
|
e*******n 发帖数: 4912 | 9 101. 1917年,日本数学家挂谷宗一(Kakeya,S(1886-1947))提出
一个问题:一位武士上厕所时遭到袭击,他只有一根短棒,为了
挡住射击,短棒应旋转360度(支点可以变化),但厕所很小,
问短棒最少要扫过多大面积?
这个问题引起当时很多人的兴趣,如1925年Birkhoff在他
写的The origin, nature, and influence of relativity
一书中提到“近几年日本数学家挂谷宗一提出的问题,是同样令
人感兴趣的问题”
1928年,苏联数学家Besicovitch,Abram Samoilovitch
(1891-1970)解决了这个问题,答案是可以任意小,1960s
Besicovitch在美国数学会就挂谷问题专门做了期科普电影,
有意思的是中间他打了个喷嚏,Besicovitch觉得很不雅,
坚持要求在录像中把这个镜头剪掉,于是现在人们看到的是
镜头突然转向一边,然后是一声闷响....
102. William James seemed to have what seemed vitally important
ideas in d... 阅读全帖 |
|
f***a 发帖数: 329 | 10 So I think you are looking for a statistical test employed with two
independent samples. Depends what you want to test, there are several tests:
for testing if they have same mean, use t Test;
for testing if they have same median, use Mann-Whitney U Test;
for testing if they have same distribution, use Kolmogorov-Smirnov Test;
for testing equal variability, use Siegel-Tukey Test or Moses Test;
Hope this helps, cheers! |
|
g********r 发帖数: 8017 | 11 兄弟胸无大志,是为了不想做千老才学统计的。从出发点来说,就已经等而下之了。转
统计之前,也在别的行业作了几年研究,最后虽然混了个统计学位,还是没有统计的视
角。做的是应用,也一直不敢跟人说自己是搞统计的。一定要强调是做应用的。这几天
看见大伙讨论统计的实质,收获很多,也忍不住说几句外行话。
从外面看统计,它是个桥梁学科,就是联系理论和应用的学科。所以统计里面偏理论的
有,偏应用的也有。打个比方,人要盖房子。出工修房子的就是搞工程的,搞医学的;
研究怎么冶铁,怎么混合粘土才能烧出砖来的就是搞数学物理的。统计呢?就是拿了别
人炼好的钢,研究做出什么工具好使的,铁锹,铁镐,铁桶什么的等等工具在修房子的
时候都用得上。工具匠人,有的更关心选择钢的材质,有的更关心工具的尺码形状.意见
不统一总是难免.
既然是制造工具的,就必须要看这个工具的目的。面对现实问题的时候,就分化出两类
统计学家来。一种是坚守理论阵地的:“你的这个问题呢,我只能根据现有理论方法,
给你解决一个相似的问题。因为我们的理论延伸不到你的问题那里。”什么时候能延伸
到呢?十年二十年都有可能。Tukey同学有个著名的句子被引用多次 |
|
D******n 发帖数: 2836 | 12 make sure u do multiple comparison adjustment
such as bonferoni, sheffe, and tukey |
|
h******e 发帖数: 1791 | 13 【 以下文字转载自 Biology 讨论区 】
发信人: hehehehe (入了大圣教,心中充满喜乐), 信区: Biology
标 题: multiple comparison方法的选择。
发信站: BBS 未名空间站 (Sun Feb 21 22:41:55 2010, 美东)
有什么根据吗?在什么情况下选择tukey,什么时候用bonferroni,什么时候用scheffe
?谢谢。 |
|
s*r 发帖数: 2757 | 14 i guess tukey and sheffe are for pairwise comparisons between multiple
levels of two factors
scheffe |
|
d*******1 发帖数: 854 | 15 google tukey trend test. i recalled ppl wrote a SAS macro for this on SUGI..
... |
|
f*******r 发帖数: 383 | 16 Fisher, Pearson, Neyman
Tukey, Rao, Wald, Le Cam, Stein
Efron, Huber, Donoho, ... |
|
a****m 发帖数: 693 | 17 如果我想做两个group, control group(16 patients), treatment group (14
different patient given the treatment), 而且对每个patient做五个不同时间的测
量,如
果对于任何时间点的对照组和treatment,看有没有差异,用什么teset?
一般的t test with bonferroni correct是不是太严格了?
还有Fisher's LSD and Tukey's HSD? 好像很多人perfer Turkey's HSD.
还有什么repeated ANOVA,不知道谁有这方面经验?
3x |
|
y*********e 发帖数: 298 | 18 我刚刚接触统计这方面,有没有入门级介绍
False discovery rate
one way ANOVA
two way ANOVA
tukey's test
Bonferroni correction
之类的书,中英文的都可以。中文更好,入门嘛,便于理解。
哪位能给个提示。
谢谢 谢谢 |
|
s*****r 发帖数: 790 | 19 有个笑话说: 没有正确的模型,但是我们也从不犯错, 因为我们从不说某个东西是错的
或者正确的。
所以, statistically we never deny dangers from anything,we just say how
likely it is.
你说“Tukey的脑海里早已有了这个概念”,“这个概念”是什么?我不认为turkey的
概念和你说的有一点点关系(what I mean is statistically unrelated)。 |
|
a***g 发帖数: 2761 | 20 多余理那帮看见首页就进来喷一嘴的人
在他们看来这和甲板的小三坑基本是一个性质
的言辞对我没有任何意义,因为我自认为我的肉体早已死过more than 10 times,我在
这里是用我的灵魂在与所有人对话。你们在任何地方看到的我的画面不过是一个灵魂幻
化成的图片而已。
信一个国内医学院毕业的卫生统计学硕士以独自一人的能力为统计学筑起了一道崭新的
地平线!毫无疑问,这是科学史上的一个奇迹。
理终将获胜,因为随机系统不可被假定;我们唯一能够假定的是,它是非确定的,因而
任何确定性假设下的方法论都是不可接受的。
Danger of Optimization’,遗憾的是,由于数学系统中长期缺乏对‘随机对应’这一
重要概念的准确定义,因而,Tukey无法从理论上论证最优化的错误。这个工作最终由
我在2009年以不是100%地严谨的方式初步完成了。
否定足以引发整个学术界的地震,因为在我看来,最优化是数学家们在随机系统中犯下
的另一个严重错误或悲剧。 |
|
|
A*******s 发帖数: 3942 | 22 why? if my memory serves me right Tukey should be more conservative than
Bonferroni if the number of tests is small. |
|
k*******a 发帖数: 772 | 23 不是很了解楼主的情况,假设这4组是独立的, 那么一共4组,所以可能的pair数是6组
. 如果同时比较这6组,那么Tukey最好,不过楼主只比较其中的4组,所以得具体计算
,假设error的df=20
B=qt(1-.05/(2*4),20)=2.74
T=qt(.95,4,20)/sqrt(2)=2.80
所以Bonferroni的确好点:)
不过如果 df比较小,比如6,那么T会比B小 |
|
A*******s 发帖数: 3942 | 24 what is ur purpose?
using Bonferroni means u test a family of hypotheses simultaneously. Say if
you would like to see which pair of groups are more/less different, then do
pairwise comparison one by one and rank the P-values. If you would like to
test a hypothesis that all/some pairs are equal, then do simultaneous test
with Bonf/Tukey/Scheffe.
From what u said before, i think u just want to find out the significantly
different pairs. u don't need Bonferroni.
Just my 2 cents. |
|
a****g 发帖数: 8131 | 25 SAS中好像没有像tukey,bon这样的可以对nonparamtric work的东西
直接降低critical value又too conservative
大家是怎么解决的? |
|
z**********i 发帖数: 12276 | 26 proc glm data=XX;
class time;
model OUTCOME = time;
lsmeans time / pdiff cl adjust=tukey;
repeated time/ type=un sub=SUBJECT;
run;
这样行吗?
多谢回复! |
|
g*****l 发帖数: 424 | 27 Formal name: Terry Speed --- The Hero of the World of Statistics
正式名: 泰瑞·速度 --- 统计世界的英雄
------------------------------------------------------------------
Another version with pictures and music is avialable at:
http://user.qzone.qq.com/176497662/blog/1332220225
题记: 谨以此文,向我们眼中那个儒雅可亲而又坚强勇敢的老人致以最高的敬意!
此文为非正式文,虽然属实,但仅供娱乐。勿拍砖。
--- by 灵致
Note: 请跳过下面两行虚晃一枪的引言直接往下阅读。
*************************** Quote: ***************************************
***
In my opinion, stat... 阅读全帖 |
|
c*********r 发帖数: 1802 | 28 没有实际遇到过,课上是这么说的。
通常interaction项可以加到error里,
如果实在要用,可以先用tukey test检验是否显著。
如果显著,就从error term里扣除。 |
|
z**********i 发帖数: 88 | 29 I will appreciate any input you may give.
I am runing glm procedure, but I get one more test significant in contrast
than in pairwise multiple comparison. I don't know why. Below is the code
and output.
proc glm data=out.final;
class dx;
format dx dx4x.;
model &v =dx;
lsmeans dx;
means dx/tukey;
contrast "demented vs normal" dx 1 0 -1;
contrast "demented vs mci" dx 1 -1 0;
contrast "mci vs normal" dx 0 1 -1;
run;
1. Least Squares Means;
Least Squares Means
Q1_PARTAPP
dx ... 阅读全帖 |
|
a***e 发帖数: 1627 | 30 题目是:
Perform the analysis of the data and identify the best treatment(一个6个,
已经量化为1-6,就是program中的trt).
SAS program :
data latin;
input row col trt yield;
cards;
1 1 3 3.10
1 2 6 5.95
1 3 1 1.75
1 4 5 6.40
1 5 2 3.85
1 6 4 5.30
2 1 2 4.80
2 2 1 2.70
2 3 3 3.... 阅读全帖 |
|
|
|
T***y 发帖数: 43 | 33 【 以下文字转载自 Mathematics 讨论区 】
发信人: Tukey (突击), 信区: Mathematics
标 题: 请推荐基于测度论的随机书
发信站: BBS 未名空间站 (Tue Aug 14 21:22:27 2012, 美东)
最近做的一个题目要用到类似Levy Measure, Ito Semi-martingale之类的东西,文献
看不懂-_-。我只学过基于测度论的概率论(具体是大数定律,中心极限定理和鞅),
想请大家介绍几本好的基于测度论的Stochastic Processes, Stochastic Calculus教
科书。最好不是干巴巴的定义一二三,定理ABC之类的,能够讲一些为什么需要定义某
个概念,某个命题在其实是说些什么之类的讲解的书。
多谢了! |
|
s*r 发帖数: 2757 | 34 http://en.wikipedia.org/wiki/Familywise_error_rate
Tukey first coined the term experimentwise error rate and "per-experiment"
error rate for the error rate that the researcher should use as a control
level in a multiple hypothesis experiment.
Since not all tests done in an experiment should constitute a single family
(for example: in a multiple-stage experiment, a separate family might be
used for each stage), the terminology was changed (by Miller) to "family-
wise error-rate" |
|
l*****o 发帖数: 61 | 35 有一组数据。Dependent variable yield, independent variables: rate, timing,
form, crop, tillage, year, site. 所有independent variables都是categorical
variable.
我用proc mixed 看这些independent variables 怎么影响dependent variable的:
proc mixed data=st0810;
class year site timing rate form Crop Tillage;
model yield=rate timing form Crop Tillage year;
random site;
lsmeans Nrate nform prevcrop tillage/adjust=t... 阅读全帖 |
|
w*******d 发帖数: 396 | 36 请教一个关于student's t test的基本问题。
实验目的:
测试维生素D是否影响钙的吸收
实验方法:
20个接收测试的人分为两组。
第一天, 第一组的人食用含有同位素钙42且不含维生素D的食物, 第二组的人食
用含有同位素钙44且含维生素D的食物。
第二天, 第二天,进行食物互换,即,第一组的人食用含有同位素钙44且含维生素
D的食物,第二组的人食用含有同位素钙42且不含维生素D的食物。
两周之后,测量血液中钙42和钙44的含量。
实验结果:
Subject 42Ca 44Ca w Vitamin D
1 a1 b1
2 a2 b2
3 a3 b3
4 a4 b4
5 a5 b5
6 a6 b6
7 a7 b7
8 a8 b8
9 a9 b9
10 a10 b10
11 a11 b11
12 a12 b12
13 a13 b13
14 a14 b14
15 a15 b15
1... 阅读全帖 |
|
w*******d 发帖数: 396 | 37 实验方法也可以这样表达。
食物A:含有同位素钙42,不含维生素D。
食物B:含有同位素钙44,含维生素D。
食物A和B 中含有同位素钙的数量相同,并且除维生素D之外,其他营养成分都相同。
1-10人:第一天吃食物A,第二天吃食物B。
11-20人:第一天吃食物B,第二天吃食物A。
两周之后,测量所有人血液中钙42和钙44的含量。
简单说就是,钙44和维生素D同时进食,钙42没有和维生素D同时进食。
体内留有的钙44和钙42量可以反映维生素D对钙吸收的影响。
我的问题是:
(1) 如果使用Student’s t test分析的话,样本是independent samples 还是 paired
samples?
(2) 是否可以用ANOVA and Tukey’s test? |
|
s****b 发帖数: 2039 | 38 比如tukey,dunnett,等。
如果已经用了two-way/multi-way ANOVA,还需要multimple comparison吗?
multiple comparison 和 multi-way ANOVA 的区别是什么呢 |
|
d********t 发帖数: 837 | 39 Do you know what mutiplicity is ?
比如tukey,dunnett,等。如果已经用了two-way/multi-way ANOVA,还需要multimple
comparison吗?multiple compar........ |
|
g*****o 发帖数: 812 | 40 你眼睛是瞎的么,还是文盲?
陈老前面说统计要加强应用,并不是说要抛弃数学。
而后面提到tukey的说法,是说他看到的很多研究既没有实用性也没有数学性。
你脸真是比地壳还厚,不知道是本来就那么厚,还是自己打脸肿成这么厚的 |
|
w*****r 发帖数: 165 | 41 求助这里的大牛:
我的数据有3个group,5个时间点。用 R 里面的 lme 来做repeated measure ANOVA:
imgdata.lme <- lme(measdat ~ Group * Time, random = ~1 | Subject,
data=imgdata, method="REML")
我现在想比较每个时间点的group间(A,B,C) 的pair-wise的区别 (A vs. B, B vs. C
, A vs. C)。不需要跨时间点比较。所以会是:
T1: A-B
T1: A-C
T1: B-C
T2: A-B
T2: A-C
T2: B-C
T3: A-B
T3: A-C
T3: B-C
T4: A-B
T4: A-C
T4: B-C
T5: A-B
T5: A-C
T5: B-C
但需要correct for multiple comparisons,比如Tukey方法。我试了用glht,但对于
该怎么设contrast matrix感觉很晕。求指点。
谢谢大家。 |
|
m**********a 发帖数: 370 | 42 想对一个数据建模,预测target跟 year county Site_ID Salinity water_temp 这5个
变量的关系。
response variable是target,另5个是exloantary variables。
现在确定 year是fixed variable,而 county Site_ID这2个地点变量是random
variable。
还有2个变量Salinity,和 water_temp。这2个变量是continous variable。
现在的问题是,这Salinity,和 water_temp 应该归为 fixed 还是random的变量?
proc glimmix data=XXX method=laplace plots=residualpanel(conditional
marginal);
class Year county Site_ID ;
model target = year Salinity water_temp /dist=lognormal link
=identity ;
... 阅读全帖 |
|