问题描述
我正在R中创建PCA绘图函数,该函数使用来自2个数据集的数据:RNA数据和表型数据。 RNA数据包括ID,样品和值
ID样本1样本2样本3
和表型数据包含样本和年龄和细胞类型等因素
采样因子1因子2因子3因子4
我要创建的是一个自定义PCA图,该图使用数据作为PCA分数
PCA_plot <- function(data,phenotype){
datapca1 <-prcomp(scale(data[-1]),center = T,scale = T)
contributions<-summary(datapca1)$importance[2,]*100
#select the scores and plot them,including the contributions
dataScores12 <-data.frame(datapca1$x[,1:2])
ggplot(dataScores12,aes(x=PC1,y=PC2))+
geom_point(#I wanna group it by Phenotype$sample)+
theme_bw(base_size=16)+theme(legend.position = "bottom")+
xlab(paste("PC1",contributions[1],"%"))+
ylab(paste("PC2",contributions[2],"%"))
}
我想按Phenotype $ sample或Phenotype $ factor4对其进行分组,如何通过ggplot创建此函数? 预先谢谢你
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)