Hi,
I know this question has been asked before a long time ago and I don't see an answer of that question in the mailing list or in the vignette of GOsemsim package. So I was wondering what is the easiest possible way of calculating GO semantic similarity value for orthologus gene pairs between two species using the above R package or any other package you know of. I am not doing this for less annotated species I need to calculate that for orthologus genes between Human and Mouse (both of which are well annotated IMHO). So I would much appreciate it if anyone who has already done this before can point me to a resource which already has pre-calculated semantic similarity values for Mouse and Human orthologues or has inbuilt code to do that.
Thanks & regards
But mgoSim also takes semData as input. How do I bypass that?
I meant what should I specify as semData. In the vignette it uses hsGO as semdata. Since here the data is from two different species what should be semData here.
As I said, merge two semData from two different species will be in TODO.
Currently, what you can do is mgoSim(measure="Wang"). In this case, see the following example:
But in this example you still have to give semData as an argument. Should I give hsGO or mmGO as semData or would both yeild similar results?
In this case, mgoSim only need the information of which ontology is using.
using hsGO or mmGO as semData will generate identical result, since species information will not used in mgoSim(measure="Wang").
Thanks that helps.
Some of the genes where number of GO terms are not equal always return 1. Is that an accurate value for functional similarity to use? what can be possible turnaround for the same other than not considering those genes altogether?
number of go terms is not a factor in semantic similarity calculation.
see https://bioconductor.org/packages/devel/bioc/vignettes/GOSemSim/inst/doc/GOSemSim.html#combine-methods
So whats explains the reason for so many gene pairs having value 1?
as you are comparing orthologus genes, it is expected.
did you read the document? If you did, you may want to try 'avg' combine method.