随笔分类 - R语言
摘要：001、无放回抽样 > a <- 1:5 ## a为数值向量 > a [1] 1 2 3 4 5 > sample(a, 3) ## 从a中任意抽取3个 [1] 3 1 2 > sample(a, 3) ## 从a中任意抽取3个 [1] 4 1 2 默认为无放回抽样，一个元素不能抽取两次， 除非本身
阅读全文
摘要：列表是R的数据结构中最为复杂的一种。 一般来说，列表就是一些对象的有序集合。 对象可以是向量、矩阵、数据框、列表等，可以使用list函数来创建列表。 > g <- "xxxxx" ## 字符串 > h <- c(2, 4, 3) ## 数值向量 > j <- matrix(1:15, nrow =
阅读全文
摘要：数组： 不同于矩阵和数据框，维度大于2。 R中最简单的数组：3维。（行， 列， 面），如下两行、三列、四面的数组。 > dim1 <- c("a1", "a2") > dim2 <- c("b1", "b2", "b3") > dim3 <- c("c1", "c2", "c3", "c4") >
阅读全文
摘要：1、 > for (i in seq(1, 10)){ + var <- paste0('var', i) ## 变量名 + data <- i + 10 ## 变量值 + assign(var, data) ## 把变量值分配给变量名 + }
阅读全文
摘要：1、 install.packages("devtools") library("devtools") install_github('royfrancis/pophelper') library(pophelper)
阅读全文
摘要：1、问题 unable to start device PNG 2、排查，进入R，输入：capabilities() 3、测试 > png("a.png") Error in .External2(C_X11, paste0("png::", filename), g$width, g$height
阅读全文
摘要：1、 > a <- c(3, 8, 2, 4, 7) > b <- c(2, 9, 4, 7, 3) > a [1] 3 8 2 4 7 > b [1] 2 9 4 7 3 > pmax(a, b) ## 返回两组向量中的最大值 [1] 3 9 4 7 7 > pmin(a, b) ## 返回两组向
阅读全文
摘要：1、测试1 > a <- 1:10 > a [1] 1 2 3 4 5 6 7 8 9 10 > rev(a) ## 逆向输出a [1] 10 9 8 7 6 5 4 3 2 1 > a <- c("aaa", "bbb", "ccc", "ddd") > a [1] "aaa" "bbb" "cc
阅读全文
摘要：1、R语言中cumsum函数用于求累积和 > a <- 1:5 > a [1] 1 2 3 4 5 > cumsum(a) ## 求累积和 [1] 1 3 6 10 15
阅读全文
摘要：1、选种需要注释的代码， ctrl + shift + c进行注释， 再次选中，执行ctrl + shift + c则取消注释 （1）选中： （2）、执行ctrl + shift + c进行注释 （3）选中 4、执行ctrl + shift + c，取消注释
阅读全文
摘要：1、测试 > par(mfrow=c(2,2)) > plot(1:10, main = "001") > plot(1:10, main = "002", xaxt = "n") ## 去除x轴 > plot(1:10, main = "003", yaxt = "n") ## 去除y轴 > pl
阅读全文
摘要：1、 #ECDF指的是Emperical Cumulative Density Function,即经验累积概率密度函数 > test <- rnorm(1000,mean=172,sd=12) ## 生成符合随机正态分布的1000个数， 平均值为172， 标准差为12 > length(test)
阅读全文
摘要：1、直接测试 > layout(matrix(c(0,1,2,3),2,2,byrow=TRUE)) ## matrix划分画布，按行排列，两行两列的矩阵， 0表示空出来，然后按照数字顺序进行绘图 > plot(1:10,main = "001") > plot(1:10,main = "002")
阅读全文
摘要：1、空白 > x <- 1:10 > y <- x^2 > plot(x,y) 2、expression用于上标 > x <- 1:10 > y <- x^2 > plot(x,y, ylab = expression(x^y)) ## 用于上标 3、用于下标 > plot(1:10, ylab =
阅读全文
摘要：1、测试1 当前目录下文件 > getwd() [1] "C:/Users/75377/Desktop/r_test2" > dir() ## 查看当前目录下文件 [1] "a.txt" > file.exists("a.txt") ## 判断当前目录下是否存在a.txt [1] TRUE > fi
阅读全文
摘要：1、quantile生成分位数 > data1 <- c(0, 10) > quantile(data1, c(0.2, 0.5)) ## 取0-10 20%和50%点的数值 20% 50% 2 5 > data2 <- c(1, 23) > quantile(data2, c(0.2, 0.5))
阅读全文
摘要：1、测试1 直接生成文件 > dir() character(0) > cat("aaa", "bbb", file = "test.txt") ## 直接生成文件 > dir() [1] "test.txt" 2、在文件上追加 > cat("1111", "222", file = "test.t
阅读全文
摘要：1、R语言中REmap安装遇到如下问题 2、解决方法 install.packages("devtools") library(devtools) install_github("Lchiffon/REmap")
阅读全文
摘要：1、 > dat1 <- c(3, 1, 4, 7, 6) ## 测试数据1， 未排序 > all(dat1 == sort(dat1)) ## dat1未排序， 判断原始向量和排序后的向量是否一致 [1] FALSE > dat2 <- c(1, 3, 4, 6, 7) ## 测试数据2， 已排序
阅读全文
摘要：1、生成矩阵，使用matrix函数 > set.seed(111) ## 设定随机数种子 > vect <- sample(1:10,16, replace = T) > vect [1] 4 3 9 5 3 8 10 1 10 4 8 10 9 8 1 7 > dat <- matrix(vect
阅读全文
浙公网安备 33010602011771号