z********n
发帖数: 710 | 1 数据是2-levels的survey design sample. 只有level-1 也就是individual level的
weight, 这种情况下面，如何在写code的时候把weighting考虑进去呢？是不是直接用
proc glimmix data = X;
weight level1weight;
model DV = ..../solution dist=binomial covb;
random intercept /subject = level2id;
run;
但是发现这样运行出来的结果都是P<.0001,显然是有问题的。但是我不知道问题在哪里
。感觉是glimmix 没有考虑sample design来计算variance。如果真的是这样话，要用
什么SAS procedure呢？ |
|
s**********7
发帖数: 24 | 2 最近在做一个project，有一个数据，是一群人的考试数据，有 ID， binary
response 0/1， 0 是错误回答，1是正确回答, 以及age, gender, years of
education 这样的variable.
Psychology/measurement 里有IRT model， 所以要求用这一堆数据和Rasch
Model (one parameter logistic Model)来估计item difficulty parameter. 这个比
较简单，用HLM，item 是level1, 人的能力是 level 2.item nested within 每个人，
PROC GLIMMIX 就得出来了。
还有一个要求是用 SAS PROC GLIMMIX，用这个binary response 数据和age
, gender, yearsofeducation variable 来model 人的能力. 我有一点没想通，人的能
力是latent trait，在HLM (Multilevel mixed-m... 阅读全帖 |
|
n**m
发帖数: 156 | 3 我现在有一个后续问题，LD要求proc glimmix出来了以后画图。
dependent variable has 3 levels. I need to plot the predicted probability,
with an independent variable on x-axis. Is there any easy way within proc
glimmix to do it?. Thanks. |
|
R******o
发帖数: 83 | 4 大概是这样的数据：
ID Time Group Outcome
1 1 A Y
2 1 A N
3 1 A N
4 1 B Y
5 1 B N
6 1 B N
7 2 A Y
8 2 A N
9 2 A Y
10 2 B N
11 2 B N
....
....
Proc glimmix data=...;
Class group?
model Outcome=Group ？
Random ID ? or Random intercept/ subject=ID ??
Random time/subject=ID Rside?
Run;
我实在很糊涂，不明白怎样同时分析group effect 和 time effect，应不应该把time
放到Model 里面。
应该是Random ID 还是 Subject=ID，两种得到的结果不一样。
而且好像 用 Rside，就必须把time 放到 Class 后面。 |
|
S******y
发帖数: 1123 | 5 What is PROC GLIMMIX (Generalized Linear Mixed Models) s equivalent in R?
Thanks. |
|
c**********5
发帖数: 653 | 6 HI,你们好，
这个是我的临床实验设计。
医院----randomization-----A or B---在每个subgroup的病人，再进行
randomization，C OR D。在接下来的12 个月内，每个病人有3 次repeated
measurements for the outcome response。
Q1。 想知道intervention effect on outcome.
OUTCOME1(BINARY)
proc glimmix ;
class Hosp_id time trt_grp patient_id;
model outcome=trt_grp|time/link=logit dist=binomial solution;
random intercept/subject=Hosp_id type=un solution;
random time/subject=patient_id type=cs residual S;
run;
Fist it i... 阅读全帖 |
|
S*****z
发帖数: 109 | 7 我可以在local run,但是提交到server (windows)后,
就出现ERROR: Procedure GLIMMIX not found.)
SAS版本是9.13.
请多指教!谢了 |
|
|
d*******o
发帖数: 493 | 9 glimmix一定要9.2版本。你用mixed不行吗 |
|
S********a
发帖数: 359 | 10 正在做generalized linear mixed model with Poisson distribution, 目前理解是用proc glimmix procedure in SAS or lmer in R, 谁能给推荐一下这方面的书，最好有例子的，谢谢 |
|
z**********i
发帖数: 12276 | 11 SAS for Mixed Models?
code and examples.
我是用NLMIXED做的。
用proc glimmix procedure in SAS or lmer in R, 谁能给推荐一下这方面的书，最好
有例子的，谢谢 |
|
n**m
发帖数: 156 | 12 有没有人知道在proc glimmix里怎么跳出 likelihood ratio or AIC, 谢谢 |
|
n**m
发帖数: 156 | 13 I'm using proc glimmix to run cumulative logit model. Boss asks me to plot
the predicted probability to make the result easier to understand. But can't
find an easy approach to do it. Any one has idea? Thanks. |
|
f*********3
发帖数: 313 | 14 用Proc Glimmix建立3-level random effect logistic regression后需要 predicted
probability，用output 命令只能给出每个样本对应的概率，但我需要自变量不同取值
的概率。请问大家有什么好的建议么？谢谢了。 |
|
q********i
发帖数: 795 | 15 加几列fake data进去，把你要的自变量值包括进去, run glimmix, output predicted
probabilities. fake data doesnt contribute to your model estimation, but
the program will calculate those probabilities in the end. |
|
S******6
发帖数: 752 | 16 用SAS run proc glimmix multilevel model (3 levels). 数据大, 13 hours
processing without result.
请教有什么好办法替代。
xiexie |
|
|
m**********a
发帖数: 370 | 18 想对一个数据建模，预测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 ;
... 阅读全帖 |
|
d**********l
发帖数: 183 | 19 尝试过用glimmix，但是因为glimmix没有repeated measure这个选项，multinomial
response必须至少在random statement中加上intercept，这样model的g matrix is
negative，因为有很多id都没有repeated measure，但是有些id有10次左右的repeated
，所以random intercept takes a lot of df.
各位高手能给点提示吗？我想model一个multinomial response，with repeated
measure，link function是genlogit，能得出baseline odds ratio的？
谢谢！！！！ |
|
o****o
发帖数: 8077 | 20 GLIMMIX里面的等效repeated statement是通过random statement来实现的，读读手册
另外你可以用NLMIX/GLIMMIX的macro，那个支持repeated statement
SAS |
|
d******3
发帖数: 93 | 21 那可以用compound symmetric的covariance matrix
就是假设所有窝里所有老鼠之间的correlation是一样的
SAS里还是用subject=来指定cluster，每个cluster里的每个obs都会被认为是一个不同
的老鼠
outcome如果不是连续的，dichotomous的话还是可以用genmod的，不是01的我没做过，
肯定可以用glimmix
不过glimmix有时候convergence有问题 |
|
s*r
发帖数: 2757 | 22 gee似乎主要还是在categorical outcome上吧，proc genmod/logistic就可以实现
但在clinical trial里，repeated measure categorical outcome似乎还是用glimmix
多，glimmix就是反复调用proc mixed |
|
w***o
发帖数: 1411 | 23 5个病人，左耳CONTROL, 右耳朵TREATMENT.FOLLOW UP 时间1，3，5 ,10月测试听力。
OBJECTIVE: 比较每个时间点TREATMENT EFFECT ON 听力。 比较每一组不同时间的听
力。
glimmix model 比较每个时间点TREATMENT EFFECT ON 听力
proc glimmix data=..;
CLASS id group month;
MODEL score= group /SOLUTION;
RANDOM INTERCEPT/SUBJECT = id;
estimate 'contral vs treatment' group 1 -1;
where mth=1;
RUN;
这个行不行？大家有什么建议？ |
|
w***o
发帖数: 1411 | 24 对，REPEATED IN MIXED IS THE SAME AS RANDOM_RESIDUAL_ IN GLIMMIX.
我的问题是，需不需要 REPEATED STATEMENT IN THE mixed MODEL or random_
residual_ in glimmix. |
|
k****i
发帖数: 347 | 25 我一般不用glimmix，如果不是normal的outcome我用GEE
你是只要比较1月和5月的区别呢，还是要比较整体时间上的趋势？
只比较control group内1月和5月的区别，用proc mixed：
proc mixed;
class month id;
model y=month /s;
random int /sub=id;
where month in (1 5);
run;
glimmix code应该类似 |
|
z**********i
发帖数: 12276 | 26 I was told "NLMIXED allows you to program your own likelihood whereas
GLIMMIX works with common distributions. I don't think GLIMMIX has the beta-
binomial distribution as one of its options and the last time I looked it
also wasn't capable to doing the zero-inflation." |
|
n*****3
发帖数: 1584 | 27 problem is , most of the company does not need
覆盖glimmix全部功能的软件包. maybe just a very
small part of it.
BTW if you use several nonlinear and mixed model
package of R, you will get most of glimmix functions. |
|
t*****u
发帖数: 26 | 28 由于爱人在Atlanta附近找到新职位,本人求ATLANTA附近的工作推荐.有５年左右的
Healthcare data analyst工作经验.SAS Programming Skill熟练.
SUMMARY OF QUALIFICATIONS
• 7 years of statistical data analysis experience.
• 7 years of SAS programming experience.
• 5 years of experience working in the Healthcare Industry.
• 5 years of experience managing and analyzing large and complex
datasets.
• 2 years of experience in ArcGIS map creation.
• National Healthcare Safety Network training a... 阅读全帖 |
|
|
s***e
发帖数: 2139 | 30 没用过,我来看一下
genmod不收殓,用glimmix就收敛了？ |
|
j*****e
发帖数: 182 | 31 Because the outcome is binary, you should use generalized linear mixed model
. The term "mixed" means you have random terms in your model. This is caused
by your design(repeated measurement). Do it either in genmod or glimmix.
The parameter estimation method is GEE.GEE is not a model. It is an
estimation method, different from the traditional likelihood method. |
|
|
|
s*r
发帖数: 2757 | 34 你有没有看过proc glimmix的manual |
|
s*********t
发帖数: 3 | 35 看了proc mixed的manual, 没有相关输出项。proc glimmix倒没看。
我目前的做法是在R中写出loglikelihood function，然后用nlm,可以得到Hessian
matrix，再求逆；
多谢大家的回复～ |
|
h****s
发帖数: 16779 | 36 是 PROC GLIMMIX. typo啊，呵呵
with |
|
h*****5
发帖数: 14 | 37 谢谢！
貌似GLIMMIX只能做MIXED EFFECT MODEL USING RANDOM STATEMENT，没有REPEATED
STATEMENT.
想找一个GEE macro，然后改写CODE， 插入一个TRUNCATED POISSON 分布。能否用SAS
实现？ |
|
h****s
发帖数: 16779 | 38 repeated / subject = ...
在Glimmix里用
random _resid_ / subject = ...
SAS |
|
h****s
发帖数: 16779 | 39 自己写一个吧
我自己写了一个proportional odds model的： GLIMMIX 在proportional odds model
的时候很多keyword不能用.
SAS |
|
h*****5
发帖数: 14 | 40 谢谢大家。刚看了下，nlinmix macro只能用于正态分布。正在看proc GLIMMIX 和
macro，希望能work。
to hypnos： 太菜了，自己写不出gee code. 方便的话，能否借鉴一下？非常感谢。 |
|
m****r
发帖数: 237 | 41 目前为止没觉得自己用过啥SAS。。。当然，还是会一些的。。不过，似乎做科研用SAS
的不多。。
我一个同学用得比较多，不过基本上也主要是在用proc glimmix和proc iml。。。至于
数据载入这些，似乎主要是以后工作要是SAS programmer的要用吧。。。
btw，统计学博士在读。。主要用R和C++来编程 |
|
|
|
S******y
发帖数: 1123 | 44 Thanks dreamer and fanta for quick reply!
That was a great help!
Will look into these in R! |
|
p******s
发帖数: 229 | 45 病人，医生，医院，3级资料。
outcome是catagorical。想控制clustering effect at 医生和医院。
以前2-级的资料，我用nlmixed, glimmix, genmod。
这3级别的，有什么SAS proc可用吗？STATA有合适的命令吗？
谢谢！ |
|
s*********e
发帖数: 1051 | 46 1) although i don't know much about "glmer", 2-level random effect is common
and very doable in sas glimmix procedure. Do take a look at the "random"
statement in this procedure.
2) GLMM "assume the data points are NOT independently distributed" is not a
correct statement. Instead, it ought to be "GLMM does NOT have to assume the
data points are independently distributed". GLMM can take care of this
dependency by various co-variance structure in either G-side matrix or R-
side matrix. |
|
w****m
发帖数: 1245 | 47 y numeric with lots of 0,
generate two variables y1 binary, y2 numeric
Binary:
if y=0, y1=0
if y^=0, y1=1
Numeric:
if y=0, y2=.
if y^=0, y2=y.
log(y2) is normal.
有牛人指教一下吗？ |
|
k****i
发帖数: 347 | 48 不用吧，他的objective已经说的很明确了
OBJECTIVE: 比较每个时间点TREATMENT EFFECT ON 听力
model我觉得没什么问题。我的问题是为什么用glimmix而不是mixed
你的听力score不是normal distribution么？ |
|
h****s
发帖数: 16779 | 49 glimmix 里面没有repeated statement,是用random statement 代替的。 |
|
w***o
发帖数: 1411 | 50 objective 2: 比方说，control group 内，听力在各个不同时间点（比方说1 月和5月
的）的比较，变化显著否？
请赐教如果用 glimmix 该如何建model. |
|
