R Studio (data analysis) R code is already done and I just need someone to write a report and add screen shot for me.

library(tidyverse) library(pROC) library(aod) pain<-read.csv("pain.csv") # 1 summdata<-pain%>%group_by(Treatment)%>%summarise( SeverityAtBaseline=sum(Severe),SuccessAfterTrial=sum(Pain)) # 2 p1<-35/100 p2<-17/100 p<-(35+17)/(100+100) a<-(p1-p2)/sqrt(p*(1-p)*(1/100+1/100)) prop.test(c(35,17),c(100,100),correct=F) # 3 pain$Treatment <- ifelse(pain$Treatment == "New Treatment", 0, 1) m <- glm(pain$Pain ~ pain$Treatment, family=binomial) summary(m) exp(cbind(OR = coef(m), confint.default(m,level=0.99))) #c-statistic pain$x<-predict(m,type=c("response")) cs<-roc(pain$Pain~pain$x) # 4 m1 <- glm(pain$Pain ~ pain$Treatment+pain$Severe+pain$Age, family=binomial) summary(m1) exp(cbind(OR = coef(m1), confint.default(m1))) #global test wald.test(b=coef(m1), Sigma=vcov(m1), Terms = 2:4) #odds ratio for a 10 year increase exp(-0.19576*10) #c-statistic pain$x1<-predict(m1,type=c("response")) cm<-roc(pain$Pain~pain$x1) # 5 plot(1-cm$specificities, cm$sensitivities, type="l", xlab="1-specificity", ylab="Sensitivity", main="Multiple Logistic") abline(a=0, b=1) grid()