Questions tagged [bioinformatics]
For programming-related questions related to Bioinformatics. Other questions do not belong here, but might be on-topic at https://bioinformatics.stackexchange.com/.
4,081
questions
0
votes
0
answers
17
views
Machine learning - classification of data
I have a ML task that might sound simple, but i am so new to this and i've just started learning so excuse my ignorance..
I have a training dataset with 104 columns and 250000 raws. I have a test ...
0
votes
0
answers
23
views
Error in using GADEM function from rGADEM package
I have big peak list in the "Bed" format and I converted it to GenomicRange for use as an input for the GADEM package to find denovo motifs. But when I try the GADEM function always I face ...
0
votes
0
answers
12
views
I was trying to run a docking simulation using rdkit library. But while running im getting a value error stating 'bad conformer id'
from rdkit import Chem
from rdkit.Chem import AllChem
from rdkit.Chem import rdMolAlign
receptor= Chem.MolFromPDBFile('4POJ.pdb')
ligand= Chem.MolFromSmiles('CC(C(=O)O)N')
rdMolAlign.AlignMol(ligand,...
-1
votes
0
answers
24
views
Invalid Syntax for Directory Pull [closed]
This is probably a very simple fix that I'm not seeing. I've tried looking around here, and on other websites to figure out what I am doing wrong. The main idea behind this project is to take a folder ...
0
votes
0
answers
9
views
How to format multiqc config file to get two different part from the same module in yaml config file?
How to format multiqc config file to get two part in HTML corresponding to two cutadapt different reports with yaml config file?
I tried different ways especially the use of files_list option. Here is ...
2
votes
2
answers
87
views
Fastest way to generate a kmer count vector from a nucleotide sequence (Julia)
Given a nucleotide sequence, I'm writing some Julia code to generate a sparse vector of (masked) kmer counts, and I would like it to run as fast as possible.
Here is my current implementation,
using ...
1
vote
0
answers
17
views
How do I download a large number of GenBank sequences using entrez_fetch in R?
I am trying to download sequence data from 1283 records in GenBank using rentrez. I'm using the following code, first to search for records fitting my criteria, then linking across databases, and ...
0
votes
1
answer
29
views
Strange adjusted (pseudo) p-values when running pairwise adonis
I really need some input on an R bioinformatics issue. I suspect this might be a case of me not completely understanding the adjustment step when doing a pairwise anova (using this wrapper https://...
1
vote
1
answer
28
views
Rentrez is pulling the wrong data from NCBI in R?
I am trying to download sequence data from E. coli samples within the state of Washington - it's about 1283 sequences, which I know is a lot. The problem that I am running into is that entrez_search ...
0
votes
0
answers
14
views
Downloading multiple sequences from GenBank in R and writing to separate FASTA files
I'm trying to download sequence data from GenBank using rentrez in R. I have performed my search using the following code:
search <- entrez_search(db = "biosample",
...
0
votes
0
answers
26
views
Error parsing strand (?) from GFF line - Happening in various programs
I'm working with various genomic data and while trying to use gffread, stringtie and cufflinks I went through the same error:
Error parsing strand (?) from GFF line: NC_037304.1 RefSeq gene ...
0
votes
0
answers
16
views
I can't install NetAct
I can´t install NetAct in R, I use got errors
I used this command:
library(devtools)
install_github("lusystemsbio/NetAct", dependencies=T, build_vignettes = T)
And got this:
library(...
0
votes
2
answers
44
views
Loop in R not working, generating single Value
I have some metabolomics data I am trying to process (validate the compounds that are actually present).
`'data.frame': 544 obs. of 48 variables:
$ X : int 1 2 3 4 5 6 7 8 9 10 ...
...
1
vote
1
answer
40
views
How to pass the one output from one process to another process in Nextflow?
So I am new to nextflow, I am trying to pass the bam file produced by ALIGNMENT process to the MKDUP process but its throwing a null value error.
I kind of understood the error, its due to the ${...
0
votes
0
answers
13
views
Call local_auth in flutter resume state will request authentication continuously
I'm using the local_auth package to lock app services after login if the user closes the app
what I want is to call the authentication function in the resume state of AppLifecycleState (when the user ...
0
votes
1
answer
52
views
Kaplan-Meier Survival curve with 3 Plots
I am very new to R studio, and I am currently learning how to do Kaplan-Meier survival curves.
Here are the columns that are needed and the information the columns contain:
Group: “normal”, “high”
...
0
votes
1
answer
29
views
How can I grep patterns in a fasta file, then print all those cooresponding sequences to a fasta file?
I've had so much help with code errors when I've used this page. I really appreciate how much everyone wants to help out!!! I am brand new to any kind of coding and whooooooh this is a learning curve.
...
0
votes
0
answers
28
views
How can I count the number of occurences in a .fasta file using a .txt file? [duplicate]
I am really new to coding generally. I'm trying to figure out how to write a bash script that can help me find the occurrences of patterns in a fasta file. The occurrences are in a .txt file, and are ...
0
votes
0
answers
28
views
Passthrough not supported error in Python
I am trying to run a python script in cmd for analyzing large Next Gen Sequencing data files, however it crashes after writing the following error to the screen:
[8640:0124/145024.529:ERROR:gpu_init....
0
votes
0
answers
13
views
Problem getting PlotMAXY to work in R, packages loaded: beadarray and limma
Can someone tell me what is wrong in an as simple as possible way, I'm a complete novice in R.
1 library (beadarray)
2 library (limma)
3 dataFile <- "/Users/neophytoskouphou/Desktop/...
0
votes
0
answers
26
views
Py script not running properly on Linux/Ubuntu
I have the following files:
SRR15013283_1.fastq SRR15013283_2.fastq SRR15013284_1.fastq SRR15013284_2.fastq
I was trying to run my python script name trimScript_Denovo.py but keeps returning blank,...
0
votes
1
answer
33
views
Downloading DNA sequence data in R using entrez_fetch: cannot retrieve query
I'm trying to download DNA sequence data from NCBI using entrez_fetch. With the following code, I perform a search for the IDs of the sequences I need with entrez_search, and then I attempt to ...
-1
votes
0
answers
53
views
Monte-Carlo simulation of protein structure in Python
I'm doing a project at a university where I have to use a monte carlo simulation to determine the structure of a protein, so far my code looks like this:
from Bio.PDB import PDBIO
import random
import ...
2
votes
2
answers
61
views
Nextflow script to process all files in given directory
I have a nextflow script that runs a couple of processes on a single vcf file. The name of the file is 'bos_taurus.vcf' and it is located in the directory /input_files/bos_taurus.vcf. The directory ...
3
votes
0
answers
34
views
Trimmomatic-0.39 error message "Unknown option -baseout ..."
I have used Trimmomatic-0.39 few times already for trimming some sequencing data. This time I tried the option -baseout as stated in the manual but it does not recognise it as a valid option and the ...
-2
votes
0
answers
22
views
Unable to paralleled running VG bioinformatic tool in parallel computing (if possible with GUI), when working on server?
I am running a server at my workplace with specs as follows:
Architecture: x86_64
CPU op-mode(s): 32-bit, 64-bit
Byte Order: Little Endian
CPU(s): 56
On-line ...
2
votes
2
answers
43
views
Extracting information from a file using another file with list of specific text; Ubuntu/ Linux
I have one file which has list of IDs, index/header, let's call it list.txt
TRINITY_DN10002_c0_g1_i1.p1
TRINITY_DN10002_c0_g1_i2.p1
TRINITY_DN10002_c1_g1_i1.p1
TRINITY_DN10006_c0_g1_i1.p1
...
2
votes
1
answer
28
views
BWA fail to locate the index files
I'm currently working on trying to analyze a dataset. I'm new to the field of bioinformatics and was trying to use BWA tools, however, as soon as I reach bwa mem, I keep running into the same error:
...
0
votes
1
answer
23
views
Iterating through a seq.Date object with a for loop to fill a Matrix with R to calculate person years
I am trying to calculate person-years for my dataset (that contains fake patients data), stratifyng them by age group. My age groups have 5 years incremental ("0-4", "5-9", "...
1
vote
2
answers
51
views
Extracting ID from file using python script?
I'm trying to extract just the sequence ID from the the file in Linux server.
To give you few examples TRINITY_DN0_c0_g1_i1.p1 and TRINITY_DN0_c0_g1_i3.p1 and sequence IDs. The sequence ID lengths are ...
0
votes
1
answer
85
views
Reorder dendrogram in R
I want to plot a heatmap where the x-axis is clustered by "normal" and "KIRP" from left to right.
Currently, my code clusters by dendrogram/similarity, and unfortunately I have one ...
0
votes
0
answers
50
views
multiGSEA package in R not appropriately assigning metabolites to pathways
I've been trying to use the multiGSEA package Vignette for multiGSEA to generate a combined p-value for pathways merging transcriptomics and metabolomics.
Even in their vignette, you can see the ...
0
votes
0
answers
16
views
Problem with creating kinship matrix in R (kinship2) without sex
I need to create a kinship matrix. For this purpose, I wanted to use the kinship2 library in R, but the sex variable is required, which I don't have. From the documentation I read that you can use the ...
0
votes
1
answer
55
views
Some operations on pandas DataFrames take hundreds of times longer on slightly larger data
I have DataFrames that have species indicators as columns, in these examples 30, and a variable amount of rows, the examples here are 450_631 (low) and 948_084 (high). In reality I have tens of ...
3
votes
3
answers
81
views
Extracting certain locus from multiple samples from text file
After profiling STR locus in a population, the output gave me 122 files each of which contains about unique 800,000 locus.
There are 2 examples of my files:
SAMPLE CHROM POS Allele_1 Allele_2 ...
0
votes
1
answer
29
views
merge multiples sequences in a single variable using Biostrings in R
I have a R list with multiples Biostrings sequences (genomes), and each element has 4000 genes (sequences) or more
class(myseq)
[1] "list"
names(myseq)
[1] "Genome-01" "...
1
vote
2
answers
31
views
String search with dynamic lengths against a file
The task I am trying to perform is: e.g. I have a string "Simultaneously", I would like to perform a search of this string with dynamic lengths against strings in a file (contains some ...
0
votes
0
answers
37
views
BWA alignment "fail to locate the index files"
This question has been asked previously, but unfortunately for me the solutions posted did not resolve my issue. I am trying to use BWA to align my ddradseq paired end reads to a reference genome, and ...
0
votes
0
answers
50
views
How to subset my Seurat object to only look at variable features while also keeping specific genes
I am relatively new to bioinformatics and am trying to work with an absolutely massive dataset that comes as a single Seurat object of about 90gb. My computer keeps crashing every time I try to even ...
1
vote
1
answer
53
views
Group rows by a column, then keep only rows with alternating values in different specified
From my previous query, I have obtained the following data using R -
info
s1
s2
s3
gene
a
Het
Het
Ref
AA
b
Het
Ref
Het
AA
c
Het
Ref
Ref
AB
d
Het
Ref
Het
AB
e
Het
Ref
Het
AC
f
Het
Ref
Het
AC
g
...
0
votes
1
answer
21
views
Error when trying to calculate mean and SD of environmental dataset with loop from .nc data
I was trying to calculate mean and SD per month of a variable from an environmental dataset (.nc file of Sea surface temp/day during 2 years) and the loop I used gives me the following error
Error in ...
0
votes
0
answers
32
views
How can I access all PyMol commands through a Python script?
Currently trying to use a Python script to automate protein modelling in PyMol, but I cannot access PyMol commands through the Python script.
Any ideas on how to solve this?
I tried doing it as below, ...
1
vote
1
answer
24
views
Convert Enrichlist to Dataframe
Is their a way to convert an enrichlist from CBNplot into a dataframe, so I can use the the function bngeneplotCustom ?
https://github.com/noriakis/CBNplot
0
votes
1
answer
94
views
How can I calculate the geometric center of a protein in Perl?
I have a PDB file which contains information about a specific protein. One of the information it holds is the positions of the different atoms (xyz coordinates).
The file is the following https://...
2
votes
3
answers
98
views
Can anyone spot the problem in my for loop in Perl?
I have a text file which contains a protein sequence in FASTA format. FASTA files have the first line which is a header and the rest is the sequence of interest. Each letter is one amino acid. I want ...
2
votes
2
answers
52
views
Snakemake/miniforge: Nothing provides libgcc-ng >=12 needed by freebayes-1.3.6-hbfe0e7f_2
I am running a Snakemake NGS data analysis pipeline, one rule uses Freebayes, it's env looks like this:
name: freebayes
channels:
- bioconda
dependencies:
- freebayes=1.3.6
While creating the env ...
0
votes
0
answers
59
views
Directory files as input in Snakemake
I plan to implement a pipeline where I search for specific transcripts at three different genomes to align the best hits and estimate some statistics by each.
To automatize the task in Snakemake, I ...
0
votes
0
answers
23
views
Getting error while installing pyseer using miniconda
I am a new Miniconda user that needs to use pyseer to perform a GWAS study in regards to a bioinformatic pipeline. After downloading 64-bit Miniconda on Windows 11, I first created a new environment (...
-2
votes
0
answers
31
views
Molecular Similarity with USR descriptors
I want to calculate the similarity between molecules using USR descriptors. I calculate USR descriptors using "https://github.com/eoc21/UltraFast-Shape-Recognition" and get 12 values against ...
0
votes
0
answers
35
views
What's wrong with my "Greedy Motif Search" code?
I'm studying bioinformatics at coursera. Here's me in "Greedy Motif Search" code.
Input ans Output is like this.
# Input: A list of kmers Dna, and integers k and t (where t is the number of ...