YJJ18636810884的博客

文章目录R语言绘图1难点一：修改图像参数的两种方法难点二：如何将多条曲线画在一张图上难点三：给点添加文本标注难点四：如何将图片导出为不同格式R语言绘图1难点一：修改图像参数的两种方法#方法一：用par参数进行修改#par()将生成一个含有当前图像参数设置的列表#par(no.readonly = T)可以生成一个可以修改的当前参数列表#输入运行数据dose = c(20, 30, 4...

x <- c(4,6,2)1、 sort(x)是对向量x进行排序，返回值排序后的数值向量, 即(2,4,6)。2、 rank()是求秩的函数，它的返回值是这个向量中对应元素的“排名” (2,3,1)。3、 order()的返回值是对应“排名”的元素所在向量中的位置，即(3,1,2) 当然可以设定decreasing 来设置是从大到小还是从小到大，default是从小到大的顺序。...

如果列表中间的某一元素没有输入，则将该元素自动保存为NULLa <- NULLa[[1]] <- matrix(1:12, ncol = 3)a[[3]] <- c("B","C","F")结果输出如下：a[[1]][,1] [,2] [,3][1,] 1 5 9[2,] 2 6 10[3,] 3 7 ...

R语言中的集合运算A <- 1:10B <- 1:5C <- seq(2,8,2)#求A和B的交集union(A,B)#求A和B的并集intersect(A,B)#A-Bsetdiff(A,B)#B-Asetdiff(B,A)#判断A是否属于Ball(A%in%B)#判断B是否属于Aall(B%in%A)#判断集合A和B是否相同setequal(...

用PCA的方法确定最佳聚类数M 为计算出共识矩阵Fn = ecdf(M[lower.tri(M)]) 提取出共识矩阵下三角的数据，然后将用ecdf 方法生成拟合曲线计算0.1到0.9之间的面积面积最小值对应K为最佳KKvec = 2:maxKx1 = 0.1; x2 = 0.9 # threshold defining the intermediate sub-intervalPAC...

文章目录1. 用NMF示例方法产生默认的判断rank数的图片2. 提取cophenetic 数据3. 实现自动判断最佳rank值1. 用NMF示例方法产生默认的判断rank数的图片library(NMF)data("esGolub")esGolub <- esGolub[1:200,]esGolub$Sample <- NULLestim.r <- nmf(esGol...

第一种方法，直接用list() 函数生成列表#method onea <- "my first list"b <- 1:3c <- matrix(0, ncol=3,nrow=4)mylist <- list(name=a, other=b,c)第二种方法 先生成一个空文件，然后用[[i]] 生成列表，第二种方法在循环生成列表中用的很多a <-...

问题阐述：同样的数据，同样的命令，重复运行命令时，发现每个cluster里面的个数都在变化问题原因：k-means算法中起点是随机或者认为给定的，如果初始类中心发生改变，可能会导致结果改变解决问题的方法：1、把初始随机数进行固定2、自己设定初始随机种子R语言中代码实现R里 kmeans(x, centers, iter.max = 10)可以通过centers 来进行调节R语言初始随...

1.将可能会报错的命令行放到try()函数中，并将其赋值到一个变量名fit（可以是任意变量名）中2.如果有错误，将fit将称为一个类型为“try-error”的变量3.可以用R 'try-error' %in% class(fit)放置在if()函数中进行判断4.try() 函数一般与next 和 break 进行联用，next跳过本次循环，break跳出循环 library(Conse...

文章目录聚类树的切割方法在R语言中的实现静态切割比较简单，现在主要谈动态切割##dynamicTreeCut 函数包的使用这个包主要是探测层次聚类中的clusters2.1 执行动态切割比较重要的包（“cutreeDynamic”）input ：clutering treedissimilarity matrix that was used to produce the ...

R语言对生存分析的处理用R语言，对生存分析数据进行整理，按照相应格式进行输出题目要求：

文章目录mcl $ operator is invalid for atomic vectorsmcl $ operator is invalid for atomic vectors问题阐释：在展示mcl运行结果mcl_result$Cluster时出现如下问题：$ operator is invalid for atomic vectors问题解决：先看看mcl函数mcl...

问题来源：探寻：NCR2 表达对GBM肿瘤芯片数据中cell-cycle 和 cytokines 通路基因的相关关系通过KEGG.db 来下载KEGG相关通路的基因if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager")BiocManager::install("KEG...