library("ggpubr")ggscatter(data,x="CE10",y="OS",add="reg.line",conf.int=TRUE,cor.coef=TRUE,cor.method="pearson",xlab="CE10",ylab="OS")cor.test(formula=~CE9+OS,data=data,subset=MFP=="D")t.test(formula=CE10~MFP,data=data)#CE1CE2CE3CE4BA-51520.153685660.2767248670.194387700.044953297BA-51530.093426640.0468381600.023344770...[
阅读全文]
KmeansR#runningclusterofgenestostratifygrouplibrary(factoextra)library(cluster)dat<-read.csv("C:/Users/?/Desktop/dat.csv")df<-na.omit(dat)df<-scale(df)#DETERMINEHOWMANYCLUSTERSISOPTIMAL(numberatthecurve)#plotnumberofclustersvs.totalwithinsumofsquaresfviz_nbclust(df,kmeans,method="wss")#calculategapstatisticbasedonnumberofclustersgap_stat<-clusGap(df,FUN=kmeans,nstar...[
阅读全文]
https://bioinformatics-core-shared-training.github.io/RNAseq-R/align-and-count.nb.htmlhttps://bioinf.wehi.edu.au/Rsubread/annot/#ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_34/GRCh38.primary_assembly.genome.fa.gzlibrary(Rsubread)library(edgeR)library(limma)library(annotables)fastq.files<-list.files(path="~/Desktop/RNAseq",pattern=".fastq.gz$",full.names=T...[
阅读全文]
扁平疣特效方香附15g
木贼15g山豆根15g板蓝根15g加水烧开,凉到不会烫伤,用纱布外洗7-10天,10分/次,1-2次/天.同时口服VitB1.[
阅读全文]
Watchzaharaenglishatyoutube.Thisisthebestoneformetolearnenglish.[
阅读全文]
D<-apply(as.matrix(dat),2,as.numeric)d<-dat%>%group_by(Group)%>%summarise_if(is.numeric,mean,na.rm=TRUE)HNSC_TT<-HNSC[,grepl("01A",names(HNSC))]packageRsubreadforfastqtobam/countpackageDESeq2fordifferentlygeneexpressionpackageDOSEforGSEApackageyarrrlibrary(Hmisc)res<-rcorr(as.matrix(data))[
阅读全文]