Handling < 1 values in NormalyzerDE
1
0
Entering edit mode
@jakob-willforss-12061
Last seen 2.0 years ago
Lund University, Sweden

Dear developer,

I have problem when running NormalyzerDE as values below 1 leads to warnings that this will cause negative values during processing. But my data is not log transform, so I don't want to run using the "noLogTransform=TRUE" option assigned - I want NormalyzerDE to log-transform my data.

normalyzerde • 522 views
ADD COMMENT
0
Entering edit mode
@jakob-willforss-12061
Last seen 2.0 years ago
Lund University, Sweden

The reason values below 1 leads to issues in NormalyzerDE is that they after log transformation will become negative.

If these values are rare you could filter out these entries by doing the following. I am using the package readr for loading and writing the updated data - it could also be done with read.table but then you need to carefully look through what options to use.

> library(readr)  # Load the library readr I use to easily read and write tsv-files
> df <- read_tsv("dataset.tsv")  # Load the dataset
> rows_above_one <- apply(df[, -1], 1, min) > 1  # Calculate which rows have only values above one, the df[,-1] is assuming that your first column contains annotation - if more you need to trim them too
> df_above_one <- df[rows_above_one, ]  # Select the rows that passed the threshold
> write_tsv(df_above_one, path="dataset_above_one.tsv")  # Write a new file

This new file can now be used as input to NormalyzerDE.

library(NormalyzerDE)
normalyzer(jobName="my_job", designPath="design.tsv", dataPath="dataset_above_one.tsv")
ADD COMMENT

Login before adding your answer.

Traffic: 946 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6