label.true<-ifelse(testset[,6]=="good",1,0) for(i in 1:length(p)){ p0<-p[order(p)[i]] label.pred<-1*(p>=p0) Recall[i]<-sum(label.pred*label.true)/sum(label.true) Precision[i]<-ifelse(is.nan(sum(label.pred*label.true)/sum(label.pred)),1,sum(label.pred*label.true...
length(mtcars$mpg) ## [1] 32 sam<-sample(1:length(mtcars$mpg),size = length(mtcars$mpg)/2) mpg1<-mtcars$mpg[sam] mpg2<-mtcars$mpg[-sam] wilcox.test(mpg1,mpg2,correct = FALSE,paired = T) ## Warning in wilcox.test.default(mpg1, mpg2, correct = FALSE, paired = T): cannot...
for(i in 1:length(unique(dfsp$site))){ # 循环的次数就是i in 1:length(unique(dfsp$site)) dir[[i]]<-paste("p",unique(dfsp$site)[i],".pdf",sep = "") # dir的目的是构建循环对应图顺序的导出格式,这里是以pdf为例 # 因为我要用的导出函数为ggsave,该函数还可以导出其他格式,根据自己需...
plot(hc,labels=DM[,1],cex=.4,main=("Method:ward.D2"),xlab="") ## 作图后,可使用identify函数,直接在图中选择分层位置 Id <- identify(hc) ## 返回的Id列表,记录了分层情况,可以输出分层后的结果 for(i in 1:length(Id)){ print(DM[Id[[i]],1]) } 2. K均值聚类 数据2: setwd("/Use...
向量增加元素可以直接通过“vector[n+1]<-0”方式增加 a<-c()——向量初始化 vector <- numeric(<int>)——创建初始向量<int>个数,并赋初值为0 length(vector)<- leg——修改对象长度为leg names(vector) <- c("A","B","C")——给向量起名称 ...
方法1: 采用 for + if 语句实现 1: 自定义函数 # digital was translated into month's englishnameMonth_name_for_if<-function(month){Month_name<-c()for(i in1:length(month)){if(month[i]==1)Month_name[i]<-"Jan"if(month[i]==2)Month_name[i]<-"Feb"if(month[i]==3)Month_name[i...
ACAP1 215.31 321.62 38.617 1055.613 1790.097## ACHE 15.12 16.65 22.124 13.428 27.188## ACP5 605.90 1935.20 1120.105 306.313 744.657# 同上篇《生信数据库2: TCGA(RNA-Seq)》一致,但此处提取fpkm数据matrix=data.frame(matrix(nrow=60660,ncol=0))for(i in1:length(count_file)){path=paste0('counts/'...
2、seq,有两种用法:①seq(起点,终点,步长); ②seq(length=9, from=1, to=5) seq还有一种简写:seq(x) #相当于1:length(x),但当length(x)为0时,返回integer(0) 3、c(1,2,8) 4、使用scan(),可以等待键盘输入。输入过程中,可以使用空格分隔每个元素,也可以一行输入一个元素。输入完毕只需键入一个...
1. R对象都有两个基本属性:moede()类型属性和length()长度属性,类似的向量的类型有logical(逻辑性)、numeric(数值型)、complex(复数型)、character(字符型)。 Eg: Z<-0:9 Z<-as.character(Z) 把数值型Z转化为字符型Z 2. 缩短向量的长度: Eg: x<-c(2,4,6,8,10), ...
for (i in 1:length(les)){ vaus <- strsplit(le, " ")[[1]] wd<- vaus[[1]] fsiemgndx[[word]] = as.double(vaes[-1]) etxPressar(pb, i) } # 创建我们的嵌入矩阵 faikimbddngim = 300 fawkiebiix = array(0, c(mx_ords, faii_mdig_m)) ...