借助rmvnorm()生成一个多维正态随机矩阵
安装mvtnorm包
> install.packages("mvtnorm")> library(mvtnorm)
生成随机矩阵x=rmvnorm(n,mean=rep(0,nrow(sigma)),sigma=diag(length(mean)…)...
借助rmvnorm()生成一个多维正态随机矩阵
安装mvtnorm包
> install.packages("mvtnorm")
> library(mvtnorm)
生成随机矩阵
x=rmvnorm(n,mean=rep(0,nrow(sigma)),sigma=diag(length(mean)…)
例如
> x = rmvnorm(5,rep(0,3),diag(3))
> x
[,1] [,2] [,3]
[1,] 0.2319038 -0.22077218 0.5232609
[2,] 0.2374806 -2.33881461 -0.6796371
[3,] 0.3614914 -0.07900088 0.2848471
[4,] 1.4700363 -1.26901580 -0.4505197
[5,] -1.5877496 0.01318156 -2.1257104
mean=rep(0,nrow(sigma)，也可以改变均值向量的大小，譬如：rep(1,2)，或者rep(2,3).....甚至是负数
更多生物信息课程：
1. 文章越来越难发？是你没发现新思路，基因家族分析发2-4分文章简单快速，学习链接：基因家族分析实操课程、基因家族文献思路解读
2. 转录组数据理解不深入？图表看不懂？点击链接学习深入解读数据结果文件，学习链接：转录组（有参）结果解读；转录组（无参）结果解读
3. 转录组数据深入挖掘技能-WGCNA，提升你的文章档次，学习链接：WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘？学习链接：转录组标准分析后的数据挖掘、转录组文献解读
5. 微生物16S/ITS/18S分析原理及结果解读、OTU网络图绘制、cytoscape与网络图绘制课程
6. 生物信息入门到精通必修基础课：linux系统使用、perl入门到精通、perl语言高级、R语言入门、R语言画图
7. 医学相关数据挖掘课程，不用做实验也能发文章：TCGA-差异基因分析、GEO芯片数据挖掘、GEO芯片数据标准化、GSEA富集分析课程、TCGA临床数据生存分析、TCGA-转录因子分析、TCGA-ceRNA调控网络分析
8.其他，二代测序转录组数据自主分析、NCBI数据上传、二代测序数据解读
