由买买提看人间百态

topics

全部话题 - 话题: trt
1 2 3 下页 末页 (共3页)
a***e
发帖数: 1627
1
题目是:
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.... 阅读全帖
w*******e
发帖数: 2651
2
来自主题: Running版 - 奔TRT
TRT - Tahoe Rim Trail.
I didn't run the trail, but hiked on the trail.
Hiking on 9000 feet is quite a refreshing experience.
c*****n
发帖数: 195
3
来自主题: Chemistry版 - Fmoc-His(Trt)-OH 工业制法
最近要从同位素标记的Histidine合成Fmoc-His(Trt)-OH, 产物买不到,查了一下
Scifinder居然没有很高效的相关合成方法,有过类似经验的能否提供一下工业上大量
制成的有效方法?原料很贵,失误不起。谢过了
c*****n
发帖数: 195
4
来自主题: Chemistry版 - Fmoc-His(Trt)-OH 工业制法
谢谢,可能是我描述的有些误差,我是买到13C labeled Histidine (no protection
group), 现在想加上Fmoc 和Trt两个保护基,如果哪位有相关经验?
m***c
发帖数: 118
5
来自主题: Statistics版 - 求问一个分布
在做一个clinical trial项目,目的是研究一种treatment对于某种病发病(event)次数
的影响,数据大概如下(non-trt一共有100人,trt一共有 200人):
event trt non-trt
1 15 49
2 12 35
3 9 22
4 7 19
5 10 21
6 6 17
7 4 10
8 1 4
9 1 7
10 2 3
11 1 2
... ... ...
___________________________________
100 200
问题:
... 阅读全帖
r***o
发帖数: 127
6
两块儿地:sub1 & sub2.
在上面做实验,每年的1月1日做,从91年做到00年,一共10年(year1-year10).
在每块儿地上做test vs control,得出trt1 @ trt2.
91年做10次相同试验,n1-n10
92年做20次相同试验,n1-n20,
.....
00年做100次,n1-n100.
需要我写一个3-way mixed model的code
我的code是:
proc mixed;
class trt sub year;
model y=sub year sub*year;
random trt trt*sub trt*year trt*sub*year;
repeated ?????/ sub=sub type=cs r un;
这里的?????应该写什么啊?
在每个最小的cell里面observation是repeated.好像我以前从来没有遇到过这种问题啊.
另外,我这个的code没有错吧,人家好像很满意的样子
h****s
发帖数: 16779
7
来自主题: Statistics版 - SAS问题请教
这样应该行了:
proc iml;
use data;
read all var {ID} into ID;
read all var {AA} into AA;
read all var {Trt} into Trt;
close data;
ID_Trt1_AA3 = ID[loc(AA = 3 & Trt = 1)];
Num_ID_Trt1_AA3_unique = nrow(unique(ID_Trt1_AA3)) * ncol(unique(ID_Trt1_AA3
));
Num_Trt1 = nrow(Trt = 1) * ncol(Trt = 1);
print Num_ID_Trt1_AA3_unique Num_Trt1;
quit;
手头没SAS,不过大约就是这样。想要你的输出格式需要再控制输出流。
e****t
发帖数: 766
8
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... 阅读全帖
a******p
发帖数: 414
9
来自主题: Statistics版 - 有关sas macro 的问题
我有个问题想请教一下这里的高手:
我在用proc lifetest 画survival plot时,
我想把strata &trt; 后面加个option: strata &trt / order=internal;
但&trt是character,所以我需要用format 和 informat.
问题是&trt是macro variable, 我不知道如何用format 和 informat.
多谢!
a******p
发帖数: 414
10
来自主题: Statistics版 - 请教一个call symput的问题
下面是我的code:
data new;
set old end=last;
retain total;
if _n_=1 then all=0;
call symput(trt,compress(count));
all=all+count;
if last then call symput('all', compress(all));
run;
这里trt 是treatment,有四个值。
我想请问,call symput(trt,compress(count)); 不能用在macro中么?只能用在data
step?
下面是我想要的。
%macro count(data=,var=);
data new;
set old end=last;
retain total;
if _n_=1 then all=0;
call symput(trt,compress(count));
all=all+count;
if last then call symput('all', compress(all));
run;... 阅读全帖
a******p
发帖数: 414
11
来自主题: Statistics版 - 请教一个call symput的问题
下面是我的code:
data new;
set old end=last;
retain total;
if _n_=1 then all=0;
call symput(trt,compress(count));
all=all+count;
if last then call symput('all', compress(all));
run;
这里trt 是treatment,有四个值。
我想请问,call symput(trt,compress(count)); 不能用在macro中么?只能用在data
step?
下面是我想要的。
%macro count(data=,var=);
data new;
set old end=last;
retain total;
if _n_=1 then all=0;
call symput(trt,compress(count));
all=all+count;
if last then call symput('all', compress(all));
run;... 阅读全帖
a*******g
发帖数: 80
12
来自主题: Statistics版 - proc glm questions
I have a group of patients, 2 trt arms and two centers. the outcome is a
continuous variable. I am interested the effect of trt but like to control
the effect of centers. Here is the code I used:
proc glm data=test;
class recon treatment_arm;
model scores=treatment_arm|recon/ss3;
estimate "trt diff"
treatment_arm 1 -1
treatment_arm*recon 0.5 0.5 -0.5 -0.5;
run;
But the log returned me with
"NOTE: trt diff is not estimable.". (no error message)
When I use the proc mixed, it wor... 阅读全帖
T**u
发帖数: 54
13
来自主题: Statistics版 - SAS问题请教
我有个DATASET按ID,TRT排序后如下:
ID AA TRT
1 3 1
1 3 1
1 8 1
2 3 1
2 5 1
3 1 2
4 2 2
4 3 2
5 5 3
6 5 3
6 8 3
7 3 4
7 3 4
8 2 4
...
请问我怎么能统计在每个TRT(1,2,3,4)下, AA=3的ID(如果一个ID有1个以上AA=3,只计
算一次)的总数和百分比.
谢谢大家!
T**u
发帖数: 54
14
来自主题: Statistics版 - SAS问题请教
重排版:
我有个DATASET按ID,TRT排序后如下:
ID AA TRT
1 3 1
1 3 1
1 8 1
2 3 1
2 5 1
3 1 2
4 2 2
4 3 2
5 5 3
6 5 3
6 8 3
7 3 4
7 3 4
8 2 4
...
请问我怎么能统计在每个TRT(1,2,3,4,共有4个处理)下, 出现AA=3的ID(如果一个ID有1
个以上AA=3,只计算一次)的总数N和占那个处理的百分比.
最后结果格式为:
TRT1 TRT2 TRT3 TRT4 TOTAL
N( %) N( %) N( %) N( %) N( %)
谢谢大家!
A*********u
发帖数: 8976
15
来自主题: Statistics版 - SAS问题请教
典型的clinical trial
我给你个macro把

重排版:
我有个DATASET按ID,TRT排序后如下:
ID AA TRT
1 3 1
1 3 1
1 8 1
2 3 1
2 5 1
3 1 2
4 2 2
4 3 2
5 5 3
6 5 3
6 8 3
7 3 4
7 3 4
8 2 4
...
请问我怎么能统计在每个TRT(1,2,3,4,共有4个处理)下, 出现AA=3的ID(如果一个ID有1
个以上AA=3,只计算一次)的总数N和占那个处理的百分比.
最后结果格式为:
TRT1 TRT2 TRT3 TRT4 TOTAL
N( %) N( %) N( %) N( %) N( %)
谢谢大家!
A*********u
发帖数: 8976
16
来自主题: Statistics版 - SAS问题请教
对了,才意识到你只要AA=3,那不一定要用proc freq
用proc freq会比较容易的的把所有的AA出现过的值都统计给你
如果只要3
可以这样(trt比较多的话可以用array)
proc sort data=xxx out=xxx2 nodupkey;
by id aa;
run;
data xxx3;
set xxx2 end=eof;
by id aa;
if aa=3 then do;
if trt=1 then n1+1;
if trt=2 then n2+1;
.
.
nt+1;
end;
if first.id then n+1;
if eof;
pct1=n1/n*100;
pct2=n2/n*100;
.
.
pctt=nt/n*100;
col1=put(n1,3.)||"("||put(pct1,3.)"%)";
..
..
run;

用proc freq怎么能让一个有一个以上AA=3的ID,不重复计次呢?
谢谢!
D******n
发帖数: 2836
17
来自主题: Statistics版 - 有这么三组数。
蛮做
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;
c*********n
发帖数: 87
18
来自主题: Statistics版 - 问一个SAS lifereg的问题!
我觉得我犯了一个大错,你在1楼提到使用lifereg对吧,但是我把lifereg和lifetest
混用了,因为lifetest是nonparametric的,常用的方法包括KM和life table (2个
group),log rank test对多组,但是lifereg是parametric的,所以可能要改成:
proc lifereg ;
class trt;
model time*status(0) = trt / dist=???;
run;
你得specify一个survial time的distribution,常用的有exponential,gamma等等,
你看看这个网页:http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/statug_lifereg_sect011.htm
run完后的output里有Analysis of Parameter Estimates,你可以看到对trt effect,
看他们有没有作用,但是貌似没有pairwise的test。如果是我的话我就用l
E****2
发帖数: 247
19
请教大家一个问题:
At first, I used,
proc phreg data = pfsdays;
class trt;
model ept*eptcen(1) = trt / rl;
assess var = (trt) ph;
run;
I got : ERROR: Cannot assess the functional form of ATRSORT since it is not
a continuous explanatory variable.
Then I created a numeric var, trtpn, 0, 1 correspoding to two treatment
groups. No error. And I got
two plots for the assess statemetns, Cumulative Martingale function form,
and Standardized Score Process Plot
ods graphics on;
proc ... 阅读全帖
k****i
发帖数: 347
20
来自主题: Statistics版 - 请教个简单的统计问题(包子)
原来d是这个意思,我还以为是day...
这样的话,d1里的x1,x2应该当做repeated measures
6 out of 10个人有反应,应该按照binomial distribution的events/trials形式来分析
具体来说,应该在SAS里用proc genmod
data大概要做成这个样子
D TRT X N
1 0 n1 10
2 1 n2 10
...
SAS code大概是
proc genmod data=data;
class d;
model x/n= trt / dist=bin;
repeated subject=d / type=cs;
run;
如果想估计不同D的效果,就直接adjust for d
proc genmod data=data;
class d;
model x/n= d trt / dist=bin;
run;
e****t
发帖数: 766
21
:( 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....
s*****a
发帖数: 2735
22
来自主题: Statistics版 - 有关sas macro 的问题
试试把lifetest放到macro里?
%macro test1;
proc lifetest XXX;
strata &trt / order=internal;
....;
run;
%mend test1;
%test1;
看看这样行不行?
--------------------
不知道你的&trt是怎么来的,能把&trt产生那一步也给出吗?
a******p
发帖数: 414
23
来自主题: Statistics版 - 有关sas macro 的问题
是的,我写的是survival的macro.
&trt是macro parameter
请看这个link:
http://support.sas.com/documentation/cdl/en/statug/65328/HTML/d
我的&trt是treatment, 可能是one two three 或更多arms
我在画survival plot后,在图的最下面treatment label: one three two.
我希望是one two three.
如果我想用option, order=internal,我需要把我的&trt先format ,然后在informat,
但我不知怎么把一个macro variable format 和 informat.
多谢!
s******8
发帖数: 102
24
来自主题: Statistics版 - 有关sas macro 的问题
不清楚究竟是啥问题。“问题是&trt是macro variable, 我不知道如何用format 和
informat.”
是否可以照猫画虎?
%macro test1;
data xxx;
format &trt yourmatnamehere.;
set xxx;
run;
proc lifetest XXX;
strata &trt / order=internal;
....;
run;
%mend test1;
%test1;
a****g
发帖数: 8131
25
来自主题: Statistics版 - 用2-sided test做superiority test
这个sample size如何计算呢?
比如trt是20% (by hypothesis), control 是30%
it's easy to calculate 1-sided
when it comes to 2-sided, should i calculate the sample size?
can i do it in two calculations
trt 20% vs control 30%
and then trt 40% vs control 30%
and then pick the larger one?
thanks a lot

size
practice,
a******p
发帖数: 414
26
来自主题: Statistics版 - 请教一个call symput的问题
谢谢你的回复,但我没有写错。
如果我这样写,其中 trt 是变量,取值例如 placebo, active1, active 2,
call symput(trt, compress(count))
那么我是可以生成三个宏变量的: &placebo, &active1,&active2
如果加上引号call symput('trt', compress(count));
那么只能生产一个宏变量

发帖数: 1
27
土耳其总统埃尔多安周日表示,他23日将与所领导执政正义发展党的国会议员开会,届
时将会针对沙特阿拉伯知名记者卡舒吉(Khashoggi)丧命案,宣布土耳其方面的调查
结果。
据美联社消息,埃尔多安在伊斯坦布尔发表演说时表示,他会宣布土耳其方面调查此案
的详情,伸张正义。
据报道,埃尔多安的演说发表于沙特阿拉伯星期六早些时候的一份声明之后,该声明最
终承认59岁的卡舒吉已经在领事馆去世,尽管其解释说他在一次“斗殴”中被杀,但遭
到了国际上的怀疑和对掩护的指控。沙特阿拉伯方面还表示,涉案的18名沙特人被捕,
几名高级情报官员被解雇。
土耳其的亲政府媒体则对此作了不同的披露,报道称一支15人的沙特袭击小组前往土耳
其杀害华盛顿邮报的这位专栏作家,在几小时后乘私人飞机离开土耳其。
“这15个人为什么来这里?为什么有18人被捕?所有这些都需要在其所有细节中加以解
释,“埃尔多安质问道。稍早前,正义发展党(AK Party)副主席柯土尔穆斯(Numan
Kurtulmus)也表示,卡舒吉在沙国驻伊斯坦布尔领事馆遇害,安卡拉绝不允许任何包
庇。他还警告,本案的后果可能相当“严重”。柯土尔穆斯接受有线... 阅读全帖
j******l
发帖数: 2790
28
当别人的导师可不容易,你做的研究课题有意思有影响力你自己活的意气风发潇洒如云
,就算是烂校也会感染有抱负有激情的学生;反之如果老板形容猥琐科研上mediocre一
副混日子得过且过的样子,我真想不出来有什么可以忽悠住有理想的青年
[在 TRT (TRT) 的大作中提到:]
:这个话题可能会引起很大的争议。我指明中国学生并不是说明其他国家的学生比中国
学生好,而是作为华人教授,一般情况下,中国学生比较多。

:...........
l****a
发帖数: 23
29
【 以下文字转载自 Faculty 讨论区 】
发信人: TRT (TRT), 信区: Faculty
标 题: 中国老板, 中国学生,访学, 找工作,评tenure.
关键字: 中国学生,访学, 找工作,评tenu,中国老板
发信站: BBS 未名空间站 (Sun Apr 2 16:33:04 2017, 美东)
在这个版上看了几年。发现有几个话题讨论比较多:中国老板和中国学生,访学, 找
工作,评tenure. 很少讨论怎么写好proposal, 哪里找机会申请funding。
中国老板和中国学生的问题其实很简单。中国老板尽量不找中国学生,中国学生尽
量不要找中国老板。中国学生对中国老板有比老外更高的要求。当然中国老板对中国学
生也有不一样的要求, 觉得应该更勤奋。我现在就直接和找我的中国学生说:我知道
你们的第一选择肯定是老外,你们现在找到我是因为你们没有找到老外。大部分直接拒
绝,让他们去找别人。个别不错的,就让他们直接和组里的学生聊。如果还愿意来,就
让他们试试。合适的 就好好培养,给他们足够空间。试了不合适就直接和他们说,让他
们找别人。 如果中国学生自己试了觉得... 阅读全帖
c***d
发帖数: 142
30
来自主题: Biology版 - 请教long PCR引物设计
设计了一对用于long PCR的引物
F引物:
5'- CTC GGC AAA CAY AAA CTC CGC CTG TTT ACC AAA AAC -3'
GC content (%): 44.44 to 47.22; Tm (degrees C): 73.39 to 74.96
R引物:
5'- C AAC GAT AGT TTT TCA AGT TRT TRG TCT TGG -3'
GC content (%): 31.25 to 37.50; Tm (degrees C): 65.40 to 69.16
两个引物的GC含量相差太大,Tm值也差一些。看网上说GC含量相差太大的话,扩增不容
易启动。F引物已和别的试过,很好。R引物前后几个碱基都不适合加进去来改进引物,
别的位置也没法设计引物。所以我在5'端随机加了几个G和C以提高GC含量和Tm温度,如
下:
R引物(前面加了几个G和C,不和模版match):
5'- gcg gcC AAC GAT AGT TTT TCA AGT TRT TRG TCT TGG -3'
GC content (%): 41.67 t... 阅读全帖
s*****g
发帖数: 2666
31
来自主题: Chemistry版 - 求教一个固相多肽合成的问题
是不是cys的保护基trt很难除掉,质普里面老是有比较强的trt峰
a******r
发帖数: 786
32
ENT don't usually treat tinnitus without an underlining medical condition
such as acoustic neuroma, because there is nothing they can do. Some may try
to prescribe anti-anxiety meds but these meds themselves may have tinnitus
as a side effect. ENT usually check for any underlining medical conditions,
if it is sudden or rapid onset, they will give oral steroids for 3-6 days.
If nothing, they will sent to audiologists to do things like TRT or other
masking devices/hearing aids.
The effectiveness o... 阅读全帖
A*******r
发帖数: 194
33
来自主题: Statistics版 - what is the problem?
从treatment的第一个字符开始取一个字符。
treatment='Active'
符给trt的就是‘A’
treatment='Placebo'
符给trt的就是‘P’
q********i
发帖数: 795
34
来自主题: Statistics版 - 急问一个analysis of covariance的问题
looks to me you are doing two different tests. The first is a typical ANCOVA
(mean post score among different trt/veh,while controlling basline), but
the second is an ANOVA, coz you are simply testing the mean difference
across different category of trt/veh.
g**r
发帖数: 425
35
呵呵,离题了耶。我们公司的好像都是这么做的,我是跟着做:
用PRIMARY ENDPOINT来看,如果每个TRT少于两个SUBJECT的,就POOL。如果要做很多
POOL,按地域。只对PRIMARY ENDPOINT做POOL,其他的不管。
对于TRT×CENTER,我见过放在MODEL里面的,但很少。有这个就麻烦了,每个ENDPOINT
你都要顾及。
在PLANNING的时候,要反复跟CRP和PM强调,每个CENTER能够MEET你的最小BLOCK SIZE
很重要很重要很重要,不能MEET的时候要跟你来商量下一步怎么办。当然这个还是不能
避免小CENTER(人家就那么多人你也没办法不是?),但让他们在选CENTER的时候起码
对我们这块儿有个考虑。
i*********5
发帖数: 62
36
来自主题: Statistics版 - help rep
请教:
input individual treatment replication y1 y2;
(y1, Y2 are response variables, individual means patient, replication means
replicate each treatment 3 times, thus get 3 times y1 data )
I want to use the next model:
if no replication
proc mixed;
class treatment individual;
model y1 =treatment;
random individual;
run;
If add replication
proc mixed;
class trt individual replication;
model y1 =trt;
random individual replication, replicastion*treatment;
run;
Is this right?
我不太清楚 replication 如何处理? ... 阅读全帖
z**********i
发帖数: 12276
37
来自主题: Statistics版 - sample size calculation
算PROPORTION需要EFFECT SIZE吗?不是很确定.
但知道,他们的病人很少,应用新方法的可实施例很少.所以,他们想计算每年的有效率.
然后跟常规方法的有效率比较.
多谢回复!

then. What is the ratio of subj # in new trt vs old trt? What is the
porportion in the controlled arm?
c****y
发帖数: 94
38
把sampling 看成 block. 这样应该可以比较吧。 Sampling 1-5 就想象成 block 1-5
, 你 就用mixed model.
proc mixed;
class treatment
model response=treatment;
random sampling;
estimate "trt a vs. trt b" treatment 1 -1 0 0 0;
run;
当然你要多问问给你数据的人,这样是否恰当。 我估计 他说 “但最终又要把
sampling 1-5看成一个整体” 应该是这样吧。
s*****a
发帖数: 2735
39
来自主题: Statistics版 - 有关sas macro 的问题
有两个问题:
1. 如果不是要放在loop里,或者要生成很多同样类型的图,为什么要生成&trt。生成&
trt的目的是什么?
2. 如果只是单纯的想在最下面legend里头改变treatment label,可以把survival
curve从proc lifetest output到一个dataset,然后用proc gplot 画出来。在gplot里
头用legend statement:
legend1 across=1 frame
label=(h=2 position=bottom j=c 'Treatment label:')
value=( h=2 'One' 'Three' 'Two');
不好意思,理解不太清楚。如果你能够把你最终想要的描述一下,或者能帮你更多。
a******p
发帖数: 414
40
来自主题: Statistics版 - 有关sas macro 的问题
谢谢你的回答。我也是赶鸭子上架,表述不清楚。
下面的code是我从网上找到的。我的&trt 相当于下面code的里的group,
我的问题是,&trt是macro variable, 我不知怎么向下面一样写format
proc format;
invalue $bmtifmt 'ALL' = 1 'AML-Low Risk' = 2 'AML-High Risk' = 3;
value bmtfmt 1 = 'ALL' 2 = 'AML-Low Risk' 3 = 'AML-High Risk';
run;
data Bmt2;
set sashelp.BMT(rename=(Group=G));
Group = input(input(G, $bmtifmt.), 1.);
label Group = 'Disease Group';
format Group bmtfmt.;
run;
proc LIFETEST data=Bmt2 plots=s(atrisk(outside maxlen=13)=0 to 2500 by 500);... 阅读全帖
f*******n
发帖数: 2665
41
来自主题: Statistics版 - 请教一个call symput的问题
不需要放入宏中。
call symput( 'trt'||put(trt,1.),.......
P****D
发帖数: 11146
42
来自主题: Statistics版 - 请教一个call symput的问题
你写的逻辑不对罢了。一个宏变量只能有一个值。例如你要count的第一个值成为宏变
量trt,那就 if _n_=1 then do; call symput('trt', compress(count)); end;
f*******n
发帖数: 2665
43
来自主题: Statistics版 - 请教一个call symput的问题
不需要放入宏中。
call symput( 'trt'||put(trt,1.),.......
P****D
发帖数: 11146
44
来自主题: Statistics版 - 请教一个call symput的问题
你写的逻辑不对罢了。一个宏变量只能有一个值。例如你要count的第一个值成为宏变
量trt,那就 if _n_=1 then do; call symput('trt', compress(count)); end;
g******7
发帖数: 1433
45
来自主题: Statistics版 - 问一个R的data manipulation问题
有如下data set:
id gender TRT smoke0 smoke1 smoke2 smoke3
1 0 1 0 2 3 2
2 1 2 2 0 3 4
3 1 1 1 0 1 5
想用R create一个variable smoke 和visit成如下:
id gender TRT smoke visit
1 0 1 0 1
1 0 1 2 2
1 0 1 3 3
1 0 1 2 4
2 1 2 2 1
2 1 2 0 2
2 1 2 3 3
2 1 2 4 4 。。。。。。
请问该如何操作,谢谢!
c********h
发帖数: 330
46
来自主题: Statistics版 - a quesiton about random effect
你是说computationally?
try lme4 package in R
glmer(y ~ trt + (trt|study), family = "binomial")
Y****a
发帖数: 243
47
来自主题: Statistics版 - statistical test 问题请教
central limit theory 就是这时候用的啊,前提是你有足够多的observations
无论是diff.trt还是diff.com,看起来都不是normal distributed。但是当有足够多的
observations时,可以认为他们的mean是normal distributed。所以,可以用t-test来
检验他们的mean是否相同,即 mean(diff.trt) - mean(diff.com) = 0
e****t
发帖数: 766
48
Snatner Snell 1980 proposed the exact unconditional CI for two independent
proportion difference.
if one group is 0/n... is the CI avaible ?
it is strange for me
as i use proc freq; tables trt*y; exact riskdiff;
(i could not generate exact CI if one group is 0/n )
but if i changed to
proc freq; tables trt * y ; exact or riskdiff; run;
the exact CI for rate difference would be generated in outputs.
is this CI from Santner Snell 1980 ?
D*V
发帖数: 3096
49
来自主题: ChinaNews版 - 耳其v5:拦截载有俄公民客机
中新网10月11日电 综合报道,一家从莫斯科飞往大马士革的叙利亚客机10日晚在土
境内遭土耳其战机拦截,因为土当局怀疑该机运有非法物品。虽然客机在被截留数小时
后重新起飞,但在土叙紧张关系不断升温之际,拦截客机事件引发多方关注。俄罗斯外
交部11日发表声明,坚持要求土耳其就拦机事件作出解释。
叙客机遭拦截 被怀疑运武器
土耳其10日晚出动两架战机拦截了一架从莫斯科飞往大马士革的叙利亚客机,将其
迫降在安卡拉。机上有35名乘客和2名机组人员。客机于数小时后被土耳其放行,继续
飞往大马士革。
土耳其安纳托利亚通讯社援引达武特奥卢的话称,“根据民航条例,该客机上存在
应当汇报的非法货物,该物品被认为可疑。”他还指出,土耳其将暂时保管被没收的物
品,以进行进一步检查。但土方没有说明查获了何种物品。
据悉,两架土耳其战机当天迫使这架叙利亚航班在过境土耳其领空时在安卡拉降落
,客机上被怀疑搭载武器和弹药。
土耳其私营电视台(NTV)称,被没收的物品据信是导弹零件。而土官方电视台(TRT)
推断货物可能是运往大马士革的军用通信设备。
事件发生后,土耳其当局提醒土耳其航空公司避免使用叙利亚领空,以免发... 阅读全帖
a***s
发帖数: 12296
50
据土耳其媒体报道,土耳其最大城市伊斯坦布尔阿塔图尔克国际机场28日晚发生两起爆炸,机场内还传出枪声。
据美联社报道,目前已知近50人死亡。
当地媒体记者称,现场有手雷扔出,警方向袭击者开火。土耳其官员表示两起爆炸中至少一起为自杀式袭击。
另据央视新闻报道,两名自杀式袭击者遭到警察开枪射击后,引爆身上的炸弹。大批警察封锁机场出入口,航班全部暂停。据中国驻伊斯坦布尔总领馆透露,目前暂未接到中国公民伤亡报告。
目前,该机场所有进出航班皆已暂停,目前尚无任何组织或人员宣称策划这起袭击。
据悉,事发机场是土耳其最大国际机场。
相关阅读
土耳其总理:愿为击落俄战机支付赔偿金
据俄罗斯“卫星”新闻网6月27日消息,土耳其总理比纳利·耶尔德勒姆表示,土耳其愿意为被击落的苏-24战机向俄罗斯支付赔偿金。
俄罗斯总统新闻秘书德米特里·佩斯科夫27日向记者表示,土耳其总统埃尔多安在致俄罗斯总统普京的信函中为被击落的苏-24战机道歉,对飞行员的牺牲表示同情与慰问。埃尔多安还表示,土耳其不想破坏与俄罗斯的关系,俄罗斯是土耳其的战略伙伴。
耶尔德勒姆在土耳其TRT电视台节目中说:“信函内容很清楚。我们表达了我们的... 阅读全帖
1 2 3 下页 末页 (共3页)