摘要:
001、 协方差(Covariance): 可以通俗的理解为:两个变量在变化过程中是同方向变化?还是反方向变化?同向或反向程度如何? 计算公式: 公式简单翻译一下是:如果有X,Y两个变量,每个时刻的“X值与其均值之差”乘以“Y值与其均值之差”得到一个乘积,再对这每时刻的乘积求和并求出均值。 a <- 阅读全文
摘要:
001、 相关系数是最早由统计学家卡尔·皮尔逊设计的统计指标,是研究变量之间线性相关程度的量,一般用字母 r 表示。由于研究对象的不同,相关系数有多种定义方式,较为常用的是皮尔逊相关系数。 由于研究对象的不同,相关系数有多种定义方式,较为常用的是皮尔逊相关系数。 研究变量之间线性相关程度的量。 a 阅读全文
摘要:
001、 .Primitive 与.Internal的函数通常都是 C语言写成的,所以在R的开发环境中中是无法显示的。到网上下载R的源代码,然后在 src/main/*.c (通常位于这个位置)下查找。 参考:https://blog.csdn.net/weixin_36054993/article 阅读全文
摘要:
001、plink软件生成raw文件 root@PC1:/home/test# ls outcome.map outcome.ped ## maf过滤,并将基因型数据转换为0、1、2的形式 root@PC1:/home/test# plink --file outcome --sheep --maf 阅读全文
摘要:
001、 a <- 1:5 b <- 11:15 dat <- rbind(a, b) colnames(dat) <- paste0("c", 1:5) dat ## 测试数据框 dimnames(dat) ## 同时获取行名和列名,并生成列表 002、可以利用该函数给数据框或者矩阵命令行名和列名 阅读全文
摘要:
missing函数作用就是如果y没有具体的值,那么输出为真,如果y有具体的值,那么输出为假,这么说可能有点抽象,举个例子: 001、 myplot <- function(x,y) { if(missing(y)) { y <- x x <- 1:length(y) } plot(x,y) } my 阅读全文
摘要:
概述: 公式: 公式里的U就是下面运行结果的a$u,Σ就是下面的a$d变换成的对角矩阵(diag(a$d)),V就是下面运行结果的a$v 001、 > x <- matrix(1:20,nrow = 4,ncol = 5) ## 生成测试矩阵 > x [,1] [,2] [,3] [,4] [,5] 阅读全文
摘要:
001、 seq_len(3) ## 用于生成指定长度的1到指定长度的序列 seq_len(9) seq_len(0) 阅读全文
摘要:
相当于poyhton中assert函数, 用于嵌入函数中, 在检测某一确定的变量或条件是否为真, 如果判断为假,则为终止程序。 001、 stopifnot(3 > 0) 002、 stopifnot(3 < 0) 阅读全文
摘要:
001、 is.finite(Inf) ## 判断是否为有穷大的数 is.infinite(Inf) ## 判断是否为无穷大的数 is.nan(NaN) ## 判断是否为非数 002、 1.无穷大 R中的无穷大用Inf表示(即Infinity,无穷大),负无穷表示为-Inf。要检查一个数是否为无穷, 阅读全文
摘要:
001、 a <- c(8, 2, 3, 7, 6, 1, 2, 4) b <- c(2, 3, 6, 2, 6, 9, 3, 7) c <- c(9, 2, 7, 1, 4, 2, 8, 4) dat <- data.frame(a, b, c) ## 测试数据框 dat dat2 <- scal 阅读全文
摘要:
001、 使用methods函数 methods("prcomp") ## moethods函数返回的带*号的函数不能通过直接输入函数名获取函数的源代码,需要使用getAnywhere函数 002、使用getAnywhere函数 getAnywhere(prcomp.default) 参考:http 阅读全文
摘要:
001、 root@PC1:/home/test# ls outcome.map outcome.ped root@PC1:/home/test# cat outcome.map ## 测试数据,均为24号染色体数据 24 snp1 0 260996 24 snp2 0 363750 24 snp3 阅读全文