library(survival)
anemia <- read.delim('https://raw.githubusercontent.com/IowaBiostat/data-sets/main/anemia/anemia.txt')
S <- with(anemia, Surv(Time,Status!=0))
fit <- survfit(S~1)

table1<- cbind(fit$time[1:5], fit$n.risk[1:5], fit$n.event[1:5], round((fit$n.risk[1:5]-fit$n.event[1:5])/fit$n.risk[1:5],4), round(fit$surv[1:5],4))

colnames(table1)<- c("time", "n(t)", "d(t)", "[n(t)-d(t)]/n(t)", "cumproduct")

table1

plot(fit, ylab = "Probability", xlab = "Time")

fit2 <- with(anemia, survfit(S~Trt))
plot(fit2, ylab = "Overall Survival", xlab = "Time", col =
c("red","blue"))
legend("bottomleft", c("MTX","MTX + CSP"), text.col = c("red","blue"),
bty = "n")

survdiff(Surv(anemia$Time, anemia$Status) ~ anemia$Trt)

boringTable<-as.data.frame(matrix(data=c("A","B","C","D"),nrow=2,byrow=T),
                           row.names=(c("Reject","Fail to Reject")))
colnames(boringTable)<-c("H0 True","H0 False")
print(boringTable)

table<- rbind(c(358, 2), c(10, 886))

colnames(table)<- c("HIV", "No HIV")

rownames(table)<- c("Test +", "Test -")

print(table)

(z <- (100-90)/8) 
(p <- pnorm(z,lower.tail=FALSE))

1-dbinom(0,5,p)

(z <- (100-90)/(8/sqrt(5)))
(p <- pnorm(z,lower.tail=FALSE))