我有一个名为“unique_tissues”的向量列表,其中每个向量(列表中的元素)表示一种组织类型,向量的元素是样本名称。
plot1<-function(a, b, c){
a$x <- seq_len(nrow(a))
ggplot(a, aes_q(y=substitute(b), x=quote(x))) +
geom_point()+
geom_hline(yintercept = c)
}
plot1(s, t, 10)
文件看起来都一样,我想下载所有文件,并将它们合并为每个组织类型的一个数据框。
所有文件都存在于名为“directory_t”的文件夹中,并且所有文件都使用两列“target_id”和“TPM”共享相同的格式。
我试过用:
$liver
[1] "liver_10a_P1973_515.tsv" "liver_10b_P1973_519.tsv" "liver_10c_P1973_520.tsv" "liver_9a_P1263_202.tsv" "liver_9b_P1263_210.tsv" "liver_9c_P1263_218.tsv"
[7] "liver_9d_P1263_226.tsv" "liver_a_P000_023.tsv" "liver_c_P000_025.tsv" "liver_d_P000_026.tsv"
$lung
[1] "lung_10a_P1973_510.tsv" "lung_11a_P2952_102.tsv" "lung_11c_P2952_104.tsv" "lung_11e_P2952_106.tsv" "lung_3e_P262_148.tsv" "lung_3f_P262_149.tsv" "lung_4a_P282_111.tsv"
[8] "lung_4b_P282_112.tsv" "lung_4d_P282_114.tsv"
$prostate
[1] "prostate_10a_P1973_502.tsv" "prostate_10b_P1973_503.tsv" "prostate_10c_P1973_504.tsv" "prostate_10d_P1973_507.tsv" "prostate_10f_P1973_506.tsv" "prostate_4a_P282_108.tsv"
[7] "prostate_4b_P282_109.tsv" "prostate_4c_P282_110.tsv" "prostate_9b_P1263_242.tsv" "prostate_9d_P1263_301.tsv" "prostate_a_P189_111.tsv"
有没有办法以更好的方式做到这一点?