| 共 224 篇文章 |
|
在R软件pwr包中用pwr.chisq.test()函数可以评估卡方检验的功效、效应值和所需的样本大小。格式为：pwr.chisq.test(w = NULL, N = NULL, df = NULL, sig.level = 0.05, power = NULL)，其中w是效应值，N是总样本大小，df是自由度。pwr.chisq.test(w=ES.w1(P0,P1),N =, df = 1, sig.level = 0.05, power = 0.90)。代码如下，pwr.chisq.test (w=ES... 阅1 转自田瑞t59i6... 公众公开 22-11-13 15:47 |
R语言系列第四期：③R语言表格数据率的比较。在这个部分我们会介绍一系列用于分析表格数据的函数，我们会着重看prop.test()，binom.test() , chisq.test() 以及 fisher.test()函数。> caesar.shoe.yes<-caesar.shoe[‘Yes’,]> prop.test(caesar.shoe.yes,caesar.shoe.total)为了分析行列都多于两个分类的表格数据，可以使用函数chisq... 阅1 转自新用户864... 公众公开 22-11-13 15:45 |
Call:coxph(formula = Surv(days, censor) ~ group) n= 25, number of events= 20 coef exp(coef) se(coef) z Pr(>|z|)group 0.3652 1.4408 0.4592 0.795 0.426 exp(coef) exp(-coef) lower .95 upper .95group 1.441 0.694 0.5858 3.544Concordance= 0.536 (se = 0.064 )Likelihood ratio test= 0.64 on 1 df, p=0.4Wald test = 0.63 on 1 ... 阅78 转1 评0 公众公开 22-11-10 16:53 |
> #使用table()函数生成简单的频数统计表示例如下> mytable <- table(Arthritis$Improved)> mytable None Some Marked 42 14 28 >#使用prop.table()将这些频数转化为比例值> #注意一维列联表不用指定margins(行或者列)> prop.table(mytable) None Some Marked 0.5000000 0.1666667 0.3333333 >#使用prop.table()*100转... 阅189 转1 评0 公众公开 22-11-10 16:20 |
R之四格表与列联表数据的统计处理四格表与列联表数据的统计处理。Example13_13 <- read.table ("example13_13.csv", header=TRUE, sep=",")attach(Example13_13)library(gmodels)CrossTable(a, b)mytable <- xtabs(~a + b)chisq.test(mytable)fit <- lm(a~b)summary(fit)coefficients(fit)#回归系数confint(fit)#... 阅186 转1 评0 公众公开 22-11-10 16:17 |
> #使用table()函数生成简单的频数统计表示例如下> mytable <- table(Arthritis$Improved)> mytable None Some Marked 42 14 28 >#使用prop.table()将这些频数转化为比例值> #注意一维列联表不用指定margins(行或者列)> prop.table(mytable) None Some Marked 0.5000000 0.1666667 0.3333333 >#使用prop.table()*100转... 阅76 转0 评0 公众公开 22-11-10 15:45 |
通过上述χ2检验，明确了三种食物的真菌检出率并不相同，此时我们还需要进一步考虑三种食物真菌检出率到底谁与谁之间的差异存在统计学意义，这里就需要用到“卡方分割”，通俗讲就是把R*C列联表拆分成若干个四个表分别进行χ2检验，进而判断不同组两两比较差异是否用统计学意义，但是，因为多组比较可能会增加犯I类错误概率，所以还需要对χ2检... 阅177 转2 评0 公众公开 22-11-03 11:29 |
# 运行结果：> data77 <- matrix(c(42,30,48,72,21,36),nrow=2,dimnames=list(c("DN组","无DN组"),c("DD","ID","II")))> chi.result77 <- chisq.test(data77)> chi.result77$expected DD ID IIDN组 32.09639 53.49398 25.40964无DN组 39.90361 66... 阅1 转自勤悦轩 公众公开 22-11-02 18:39 |
按照既定的研究方案开展临床试验、收集数据，具体说来包括以下几个内容：3. 根据研究方案的规定进行期中分析（如果有），按照研究方案的规定随访患者，收集数据，填写CRF表格。2. 按照研究方案的规定进行最终的统计分析，所有亚组分析、敏感性分析都应该按照研究方案和统计分析计划的规定进行，不建议进行无计划的亚组分析。3. 按照研究方案中... 阅2333 转7 评0 公众公开 22-05-01 10:50 |
一张图搞懂所有临床研究分类。(2). 观察性研究，传统意义上的观察性研究主要用于推断病因，可分为四类：描述性观察研究、横断面研究、病例-对照研究与队列研究。真实世界研究是相对于理想世界研究而言的，一言以蔽之，可以把传统的随机对照试验意外的干预性与观察性研究均理解为真实世界研究。真实世界研究中也可以有随机对照试验，只是这种随... 阅346 转2 评0 公众公开 22-05-01 10:48 |
