Entering edit mode
hello, I am trying to extract the fold change of several genes from a list resFilt and I tried the use the code below, but I got an error. is the code for extracting the genes is wrong?
myterms <-read.table("GO-0004553.txt", header= FALSE)
head(myterms)
typeof(resFilt)
head(resFilt[[1]])
de<- resFilt[resFilt$table %in% myterms$V1]
#
>myterms <-read.table("GO-0004553.txt", header= FALSE)
> head(myterms)
V1
1 ATCC64974_100070-T1
2 ATCC64974_10330-T1
3 ATCC64974_103600-T1
4 ATCC64974_104830-T1
5 ATCC64974_104980-T1
6 ATCC64974_105530-T1
> typeof(resFilt)
[1] "list"
> head(resFilt[[1]])
logFC logCPM PValue FDR
ATCC64974_37230-T1 7.884895 8.768304 6.165094e-129 6.172492e-125
ATCC64974_22790-T1 8.445245 8.551614 6.821145e-123 3.414665e-119
ATCC64974_1640-T1 7.064251 8.514810 4.084580e-122 1.363161e-118
ATCC64974_83150-T1 11.827867 9.097794 2.269826e-108 5.681375e-105
ATCC64974_17880-T1 -10.040251 5.026719 5.238859e-98 1.049029e-94
ATCC64974_17840-T1 -7.017230 5.483155 4.549345e-97 7.591341e-94
> de<- resFilt[resFilt$table %in% myterms$V1]
Error: Two subscripts required
sessionInfo( )
thank you a lot and thank you for your advice.
If there are further questions for general bioinformatics things I recommend to post from now on at
biostars.org
. This forum here is for specific technical help with Bioconductor packages. Please also try to not open multiple questions on the exact same topic. You can always comment on your existing ones and ask for clarification if answers are not sufficient.