# use calculate_summary script to generate stats
# load stats and dictionary
#set path for file output
#path <- '/home/ahopkins/Documents/emj/ImmunoseqResults/adjuvant_study/plots/'
dir.create(path)
pdf(paste0(path,'number-vs-richness.pdf'),height=6,width=6, title='Number vs Richness')
plot(dict[-c(1,2),'cells'],stats$Richness,
xlab='Input Cell Number',
ylab='Richness')
dev.off()
pdf(paste0(path, 'number-vs-clonality.pdf'),height=6,width=6, title='Number vs Clonality')
plot(dict[-c(1,2),'cells'],stats$Clonality,
xlab='Input Cell Number',
ylab='Clonality')
dev.off()
#pdf(paste0(path, 'age-vs-richness.pdf'),height=6,width=6, title='Age vs Richness')
#plot(dict[-c(1,2),'age'],stats$Richness,
# xlab='Age at Diagnosis',
# ylab='Richness')
#dev.off()
#pdf(paste0(path, 'age-vs-clonality.pdf'),height=6,width=6, title='Age vs Clonality')
#plot(dict[-c(1,2),'age'],stats$Clonality,
# xlab='Age at Diagnosis',
# ylab='Clonality')
#dev.off()
#pdf(paste0(path, 'age-vs-overlap.pdf'),height=6,width=6, title='Age vs Overlap')
#plot(dict[-c(1,2),'age'],stats$overlap,
# xlab='Age at Diagnosis',
# ylab='Overlap')
#dev.off()
#pdf(paste0(path , 'number-vs-overlap.pdf'),height=6,width=6, title='Number vs Overlap')
#plot(dict[-c(1,2),'cells'],stats$overlap,
# xlab='Input Cell Number',
# ylab='Overlap')
#dev.off()
#pdf(paste0(path, 'clonality-vs-overlap.pdf'),height=6,width=6, title='Clonality vs Overlap')
#plot(stats$Clonality,stats$overlap,
# xlab='Clonality',
# ylab='Overlap')
#dev.off()
#pdf(paste0(path, 'richness-vs-overlap.pdf'),height=6,width=6, title='Richness vs Overlap')
#plot(stats$Richness,stats$overlap,
# xlab='Richness',
# ylab='Overlap')
#dev.off()
##################################Confounders######################################
#t <- t.test(dict$age[dict$response=='NR' & dict$type=='PRE'],dict$age[dict$response=='R' & dict$type=='PRE'])
#pdf(paste0(path, 'age-by-response.pdf'),height=8,width=4.5, title='Age by Response')
#stripchart(dict$age[-c(1,2)] ~ factor(dict$response[-c(1,2)]),#
# pch=19,
# main=paste('p=',round(t$p.value,4)),
# vertical=TRUE,
# at=c(1.25,1.75),
# xlim=c(1,2),
# ylab='Age',
# par(mai=c(1,1,1,1))
# )
#dev.off()
tt <- t.test(dict$cells[dict$response=='NR' & (dict$type=='PRETREAT')],
dict$cells[dict$response=='R' & (dict$type=='PRETREAT')])
pdf(paste0(path, 'number-by-response.pdf'),height=8,width=4.5, title='Number by Respone')
stripchart(dict$cells[-c(1,2)] ~ factor(dict$response[-c(1,2)]),
pch=19,
main=paste('p=',round(tt$p.value,4)),
vertical=TRUE,
at=c(1.25,1.75),
xlim=c(1,2),
ylab='Input Cell Number',
par(mai=c(1,1,1,1))
)
dev.off()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.