c*******n 发帖数: 892 | 1 y是binary,positive/negative
x是一个categorical variable,有20个level,每个level有30-50个obs不等
现在想比较x小于40个obs的那组negative rate是不是比大于40个obs的那组大,或者是
说x每个level的negative rate会不会随着obs的增多而减小?
请教一下,能不能把x每个level的rate和frequency拿出来作为一个新的data,根据
frequency大于40小于40分组做ttest 或者直接rate和frequency做regression?
如果不能,那应该怎么做呢? | w*******9 发帖数: 1433 | 2 坐等大牛来解释。
个人觉得直接拿rate做regression不妥,首先每个rate基于的obs个数不一样,所以每
个rate的
reliability (variance) 不一样,更重要的是你的obs的个数怎么来的肯定能起到决定
性作用。obs不是Independent variable!
如果你能自己编程的话,可以考虑以下方法:
1) Given ni, the rate pi follows G(n_i). (e.g. Beta regression for
fractional data)
2) Given pi, Yij follows binomial B(pi)
3) 1)and 2) determines P(Yij>0|n_i)=int_0^1 pG(p,f(n_i))dp
4) Write down the conditional likelihood of Yij given ni, i=1,2,.., j=1, 2,
..ni.
5) Maximize the likelihood wrt to parameters contained in G(). Assess the
correlation between pi and ni through G().
【在 c*******n 的大作中提到】 : y是binary,positive/negative : x是一个categorical variable,有20个level,每个level有30-50个obs不等 : 现在想比较x小于40个obs的那组negative rate是不是比大于40个obs的那组大,或者是 : 说x每个level的negative rate会不会随着obs的增多而减小? : 请教一下,能不能把x每个level的rate和frequency拿出来作为一个新的data,根据 : frequency大于40小于40分组做ttest 或者直接rate和frequency做regression? : 如果不能,那应该怎么做呢?
| c*******n 发帖数: 892 | 3 多谢楼上的回答!
可是我不会编程,只会用SAS做些简单的分析。。。有没有其他的方法呢?望大牛们指
点! |
|