a****g
发帖数: 8131 | 1 model中有4个factors
对于这个factor, global test结果显示有差异,但是lsmeans找不到差别,means可以
对于这种情况,都应该用lsmeans吗?
请赐教 |
|
a********a
发帖数: 346 | 2 我觉得lsmean 的结果可靠，如果不balance 的话，lsmean和mean的结果就不一样。 |
|
v*****c
发帖数: 44 | 3 How to call for overall adjusted means and standard error in Proc GLM? Is
there any short forms other than the 'estimate' statement below (I have 9
levels of analsite)?
proc glm;
class treat analsite;
model pchg = analsite treat/ ss3;
lsmeans treat /pdiff cl stderr e;
estimate 'overall' intercept 1 analsite 0.11111111 0.11111111 0.
11111111 0.11111111 0.11111111 0.11111111 0.11111111 0.11111111 0.11111111
treat 0.25 0.25 0.25 0.25 ;
run; |
|
m*********n
发帖数: 413 | 4 你说的global test是？
lsmeans和means是两种不同的contrast，看你要用哪一种。 |
|
x***I
发帖数: 91 | 5 9.2的help说proc mixed的lsmeans有个at的option. 可是我试了很多次都不好用，有人
成功用过这个么，能给个例子么，谢谢亚。 |
|
d****z
发帖数: 301 | 6 一直折腾不出来.好几个月了.
牛细胞生长情况, 有两种细胞,3个时间点,每个有两个重复.
想找出差异表达的基因...
结果一直疯狂报错:
warning: Stopped because of infinite likelihood.
跪求建议意见!!!!
proc mixed data=long;
class rep genotype time;
model y=rep genotype time genotype*time /outp=check;
random rep*time;
lsmeans genotype*time;
ods output tests3=results;
ods output lsmeans=lsmeans;
by gene;
run; |
|
w*****n
发帖数: 375 | 7 Your random statement might be wrong, and I suspect that you don't
have enough degrees of freedom to estimate so many parameters.What's
your sample size?
Generate variable "id", which has two reps, and try the following
simple model first:
proc mixed data=long;
class id genotype time;
model y = genotype time genotype*time /outp=check;
random intercept/sub = id;
lsmeans genotype*time;
ods output tests3=results;
ods output lsmeans=lsmeans;
by gene;
run; |
|
w******4
发帖数: 488 | 8 I think this is a split-plot design, where genotype is the whole-plot factor
, time is the sub-plot factor. You can try this:
proc mixed data=long;
class rep genotype time;
model y=genotype time genotype*time /outp=check;
random rep(genotype);
lsmeans genotype*time;
ods output tests3=results;
ods output lsmeans=lsmeans;
by gene;
run; |
|
a****y
发帖数: 1035 | 9 我想把下面几行SAS code 在R里实现，请教牛人指点一二。
请问SAS code里ddfm=kr 和 lsmeans 在R 里有类似的option吗？
proc glm data=sample;
class a b c;
model y = a b c*a;
random c*a /test;
lsmeans a /stderr pdiff;
run;
proc mixed data=sample;
class a b c;
model y = a b /ddfm = kr;
random c*a;
lsmeans a /pdiff;
run;
谢谢先！！！ |
|
c*******o
发帖数: 8869 | 10 这么说吧：
for blood pressure"
proc mixed;
model bloodpressure=treatment baseline;
lsmeans treatment/pdiff;
run;
for change of blood pressure from baseline:
proc mixed;
model change =treatment baseline;
lsmeans treatment/pdiff;
run;
I used baseline in both model. should I expect same p value for contrast?
ANCOVA |
|
a****y
发帖数: 1035 | 11 补充一下，也就是说，已有log（y）的lsmeans 了，如何得到y 的lsmeans？
欢迎指教！
谢谢！ |
|
z**********i
发帖数: 88 | 12 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 ... 阅读全帖 |
|
e****t
发帖数: 766 | 13 :( ls mean diff and SE is from statement:
lsmeans trt * time / pdiff;
for example, for time 4, trt 1 - trt 2, use this statement, we obtain
lsmeans diff and se.... |
|
m**********a
发帖数: 370 | 14 想对一个数据建模，预测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 ;
... 阅读全帖 |
|
w*****1
发帖数: 473 | 15 【 以下文字转载自 Statistics 讨论区 】
发信人: wz99331 (dotti), 信区: Statistics
标 题: 请问大家有关mixed model
发信站: BBS 未名空间站 (Wed Jul 19 12:57:18 2017, 美东)
大家好，我需要用mixed model 对三个group的人群（教育程度1，2，3）进行分析，看
看这些人的认知能力在不同时间的变化有没有显著性的不同。
proc mixed data=outban;
class timepoint educ;
model rbantotal= timepoint | educ/ solution;
random int /subject=id ;
lsmeans educ*timepoint;
run;
因为我是3个group，而不是两个group,如果我想看相互作用项educ*timepoint的p-值，
是不是需要用f-test?如何在mixed model中用f-test呢？我看到现在的结果里就已经有
f-test了：
... 阅读全帖 |
|
w*****1
发帖数: 473 | 16 【 以下文字转载自 Statistics 讨论区 】
发信人: wz99331 (dotti), 信区: Statistics
标 题: 请问大家有关mixed model
发信站: BBS 未名空间站 (Wed Jul 19 12:57:18 2017, 美东)
大家好，我需要用mixed model 对三个group的人群（教育程度1，2，3）进行分析，看
看这些人的认知能力在不同时间的变化有没有显著性的不同。
proc mixed data=outban;
class timepoint educ;
model rbantotal= timepoint | educ/ solution;
random int /subject=id ;
lsmeans educ*timepoint;
run;
因为我是3个group，而不是两个group,如果我想看相互作用项educ*timepoint的p-值，
是不是需要用f-test?如何在mixed model中用f-test呢？我看到现在的结果里就已经有
f-test了：
... 阅读全帖 |
|
|
s*r
发帖数: 2757 | 18 remove the 'class x'
or use lsmeans |
|
D******n
发帖数: 2836 | 19 蛮做
trt block y
1 1 1.1
2 1 2.1
1 2 3.3
2 2 4.5
3 3 4.5
3 4 6.7
proc mixed data=a1;
model y=trt;random block;
lsmeans trt/adjust=bon;run; |
|
h******e
发帖数: 1791 | 20 使用glm，mixed这些procedure中的estimate statement，算出的结果是什么呀？该如
何解释？于lsmean statement算出的结果有什么不同？谢谢。 |
|
q****k
发帖数: 1023 | 21 lsmeans treatment*time | slice=time?
I like the SLICE statement.
follow- |
|
c**d
发帖数: 104 | 22 when you look for mean change over time, you first have to check interaction between month and treatment. if significant interaction,
treat month as a class variable (you can create a new time variable ) to find mean change over time.
Then you can get what you want by using lsmeans trt*monthclass/
By the way, model can be complicated to check linear and quadratic trend with random effects
follow- |
|
z**********i
发帖数: 12276 | 23 proc glm data=XX;
class time;
model OUTCOME = time;
lsmeans time / pdiff cl adjust=tukey;
repeated time/ type=un sub=SUBJECT;
run;
这样行吗?
多谢回复! |
|
o****o
发帖数: 8077 | 24 分析分好多层次，从最底层的hardcore programming到最表层的大忽悠，各个层次需要
相适应的工具和技术，能赚钱的分析也从小到大。
有些项目需要大规模的计算，使用超大的数据，超高的数据更新率，这种分析就需要你
提到的那些高度技巧的能力；有些就是很简单的相关分析，要的是简单明了，能很容易
的跟管理层交流，这个时候像SAS和R这种对编程要求较低的工具就体现作用了。
其实对于工业界的分析师，我个人觉得追逐编程技巧是舍本逐末，因为随着年龄增加，
这种能力终究是比不过年轻人的；而真正应该关注的是本行业的业务知识，管理知识，
对于大部分行业这些都是年纪越大越吃香。编程嘛，招几个年轻人来就好，一天工作10
个小时不嫌多的。
另外，SAS从出生那天就不能算正规的编程语言；SAS本来就是提供一个方便的界面好用
计算机语言来描述统计分析，那些LSMEAN,CONSTRAST statement就是干这活的。所以
SAS的本意就是给你一个工具让统计师能摆脱编程的痛苦，集中精力分析数据和业务。
至于说SAS的年费，对大公司真的是不值一毛，主要看老板吹效果，而且这种软件一旦
用上了，价格弹性很小。另... 阅读全帖 |
|
a****y
发帖数: 1035 | 25 对log 的 values， 其LSmeans如何做 transform back？ |
|
l******o
发帖数: 3764 | 26 现在有一个dependent variable (dep1), 几个predictors （var1 var2 var3 var4)
var1是个categorical variable，求dep1在var1的几个level上的adjusted means
我以前都这么做（只会这一个）
proc glm data=xxx;
class var1;
model dep1=var1 var2 var3 var4;
lsmeans var1;
run;
但是这次我的var1是ordinary的（其实var1是我把一个continuous variable分成的
quintile)，我这样用是不是就把var1的几个categories当成无序的处理了，会lose
power吧？
我记得老师讲过大概应该是先得出var1=0时dep1的值，然后var1=1时dep1的值就是加上
1 x beta of var1, 以此类推，是这样吗？
可是俺不会sas code 请好心的高手指教
多谢多谢 |
|
e****t
发帖数: 766 | 27 suppose, each patient will be measured at 1,2,3,4 time, but there is missing
value at different timepoint.
we fit model as
proc mixed ;
class sub time trt ;
model sub = time trt time * trt;
repeated time / sub = sub type = un ;
lsmeans trt * time / pdiff ;
run;
when calculating effect size,
we use ls Means diff / (SE / sqrt (1/n1 + 1/n2))
here n1 and n2 should be patient number in different group.
or should be patient number in different group by visit (due to missing,
observations at each visit... 阅读全帖 |
|
i**z
发帖数: 194 | 28 unbalanced design or missing data will make a difference in the two
procedures. |
|
s*r
发帖数: 2757 | 29 proc glm 里面的means 多半是把mean和0比吧 |
|
i****e
发帖数: 151 | 30 两组人做REPEATED MEASUREMENTS(2visits)，每次做两只眼。
目的一：比较两组OUTCOME在每次VISIT的不同；
目的二：比较每组两次VISIT的变化
需要考虑到每个人两只眼的CORRELATION
PROC MIXED DATA=A;
CLASS GROUP EYE VISIT;
MODEL Y=GROUP VISIT GROUP*VISIT ;
REPEATED VISIT /SUB=ID;
RANDOM ID;
LSMEANS GROUP*VISIT;
这个CODE对吗？尤其是RANDOM？GROUP*VISIT那个自动设了一个作为REFERENCE GROUP，
我想设每组把第一个VISIT作为自己的REFERENCE怎么弄（比如程序把第一组第一个
VISIT设为REFERENCE,那么第二组的第二个VISIT跟这个比的P就没有意义） |
|
j******k
发帖数: 91 | 31 应该是你LSMEAN后面那些数字写得不对....
我现在就是面临这个问题.... |
|
|
s*r
发帖数: 2757 | 33 如果是从lsmeans statement 出来的(不是你自己写estimate statement出来的), 我所
见的原因是class variable involved in an interaction effect,然后有一个level里
面的人比较少 |
|
l*****o
发帖数: 61 | 34 有一组数据。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... 阅读全帖 |
|
a********e
发帖数: 256 | 35 比如two factor experiment的话，没有missing cell, 只是少数cell里
面的reps比别的多，Ｒ用lm, ＳＡＳ用proc mixed, 这种情况下R和SAS
的结果会差很多吗？
我的data 算出来，lsmeans是一样的，但是每个treatment level的t-value和p-value
不一样。 |
|
a********e
发帖数: 256 | 36 我的data, lsmeans算出来一样，但是t-value 和p-value不一样。 |
|
R******d
发帖数: 1436 | 37 试了一下 lsmeans time*group / slice=time ;
感觉和glm出来的结果相差很大啊。 |
|
z*****w
发帖数: 118 | 38 我的model:
proc glm data=f;
class n year;
model tc=year n n*year/solution;
lsmeans n n*year/pdiff;
estimate "compare m1-m2" n 1 -1 0;
estimate "compare m1-m3" n 1 0 -1;
estimate "compare m2-m3" n 0 1 -1;
run;
quit;
n有３组，year有6年，想test每２组之间slope的差异，不知道这样的sas code对不对？
还有estimate statement这样写，好象没有考虑interaction n*year,对吗？ |
|
w*****1
发帖数: 473 | 39 大家好，我需要用mixed model 对三个group的人群（教育程度1，2，3）进行分析，看
看这些人的认知能力在不同时间的变化有没有显著性的不同。
proc mixed data=outban;
class timepoint educ;
model rbantotal= timepoint | educ/ solution;
random int /subject=id ;
lsmeans educ*timepoint;
run;
因为我是3个group，而不是两个group,如果我想看相互作用项educ*timepoint的p-值，
是不是需要用f-test?如何在mixed model中用f-test呢？我看到现在的结果里就已经有
f-test了：
Type 3 Tests of Fixed Effects
Num Den
Effect ... 阅读全帖 |
|
w*****1
发帖数: 473 | 40 【 以下文字转载自 Statistics 讨论区 】
发信人: wz99331 (dotti), 信区: Statistics
标 题: 请问大家有关mixed model
发信站: BBS 未名空间站 (Wed Jul 19 12:57:18 2017, 美东)
大家好，我需要用mixed model 对三个group的人群（教育程度1，2，3）进行分析，看
看这些人的认知能力在不同时间的变化有没有显著性的不同。
proc mixed data=outban;
class timepoint educ;
model rbantotal= timepoint | educ/ solution;
random int /subject=id ;
lsmeans educ*timepoint;
run;
因为我是3个group，而不是两个group,如果我想看相互作用项educ*timepoint的p-值，
是不是需要用f-test?如何在mixed model中用f-test呢？我看到现在的结果里就已经有
f-test了：
... 阅读全帖 |
|
