ℍost bionimbus HostName bionimbus-pdc.opensciencedatacloud.org User username IdentityFile ~/.ssh/id_rsa ForwardAgent yes Host bionimbusvm HostName instanceIPaddress User ubuntu IdentityFile ~/.ssh/id_rsa ProxyCommand ssh -q -A bionimbus -W %h:%p ForwardAgent yes Load your key into the ssh-agent by … Read more →

𝕋he error message– INFO - 0 % of model's snps used– can typically be traced to inconsistencies between variant IDs in prediction models and input GWAS files. Our GTEx v8 mashr and elastic net models are defined in hg38. For GWAS defined in hg19, we recommend first processing the GWAS with our … Read more →

𝔾andal et al analyzed autism spectrum disorder, schizophrenia, and bipolar disorder across multiple levels of transcriptomic organization—gene expression, local splicing, transcript isoform expression, and coexpression networks for both protein-coding and noncoding genes to produce a quantitative, … Read more →

𝕎ant to write a function that takes in two data frames with expression data and returns […] vector of correlations of matched genes […] summary of the correlation vect (min,1Q,median,3Q,max) […] returns a list of aligned data frames […] plots correlation for a sample of … Read more →

𝕄any users had difficulties matching the genotype variant id to the prediction model variant id. Here is one example added to the PrediXcan tutorial where the matching was failing because of the on the fly option not taking into account that in the GTEx v8 vcf file, chromosomes are names as chr# … Read more →

𝔸 p-value of zero should be interpreted as an extremely small positive value. S-PrediXcan or PrediXcan will provide the zscore as well as the p-value. You can calculate the p-value corresponding to the Zscore using the formula below. For example a Zscore of 30 will give you a (natural) log p-value … Read more →

𝕆perating System Upgrade - The operating system will be upgraded from Red Hat Linux 6.7 to 7.6. This will provide a kernel that will allow for a more modern software ecosystem. For example, software such as tensorflow will not run on Red Hat 6. […] GPFS Upgrade - GPFS storage clients will be … Read more →

/gpfs/data/im-lab/nas40t2/Data/References/mappings/UKB2GTEx_mapping.txt.gz contains information for variants in UK Biobank genotypes. The columns are variant chromosome position non_effect_allele effect_allele rsid zscore panel_variant_id 1:692794:CA:C chr1 757414 CA C 1:692794_CA_C 1 NA … Read more →

𝕎hen working with large datasets, only the files with code should be pushed to Github repositories, not the data itself. The raw data inputs or analysis output should either be kept in a local directory that is never committed, or for best practices, they should be stored in Box (download). Once … Read more →

ℕote: You can download tensorqtl using pip install. However, there seems to be a bug that makes tensorqtl incompatible with pandas plink 2.2.2. If you want to download tensorqtl using pip, then you have to downgrade pandas plink to 2.1.x so that it matches the version of pandas 1.1.x I installed … Read more →