Could someone please give me a working example of obtaining a PPI interaction network from a list of genes in R or Python.

The input is a list of genes such as : 6416, 2318, 351, 409, 5392, 1956 .......

And the output should be two columns that represent two interacting genes, like below:

6416    2318
6416    5371
6416    351
6416    409
6416    5932
6416    1956

Also,there are thousands of gene and I want to control the number of output interaction so my laptop can handle the computations and not run out of computational power. :P

(If possible, please sow how to use BioGRID and PICKLE using R to solve the above)

Thanks.

Hi, this forum is for questions related to Bioconductor packages. You may consider a more general bioinformatics forum like Biostars.