Problems with limma package in R
2
0
Entering edit mode
d.leung • 0
@dleung-9129
Last seen 9.1 years ago

Dear all,

I would be very grateful for some desperately needed help.

I am trying to run backgroundCorrect command in R after I loade the limma package.  However, I came up with an error code saying: Could not find function "backgroundCorrect"

I am sure that I have installed the limma packages and the statmod.

Any help will be much appreciated

Regards

Dominic Leung

limma • 1.6k views
ADD COMMENT
0
Entering edit mode

Dear Efstathios,

Thank you for your help.

I am using Mac and running R 3.2.2

R version 3.2.2 (2015-08-14) -- "Fire Safety"
Copyright (C) 2015 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin13.4.0 (64-bit)


How do I find session info?

Thanks

Dominic

ADD REPLY
0
Entering edit mode

Im using Windows, but i think sessionInfo() should work

ADD REPLY
2
Entering edit mode
@gordon-smyth
Last seen 10 minutes ago
WEHI, Melbourne, Australia

Dear Dominic,

I am guessing, since you say that you're not sure how to type sessionInfo() at the R prompt, that this may be one of the first times that you are running R.

You have to understand that, so far, you're not giving us much information to help you. First you said that the error message was 'Could not find function "backgroundCorrect"'. You could only get this message if you have not in fact loaded limma. From your comments to Efstathios, I guess you later realized that you have to type library(limma) before you can start using limma functions.

Then you said that you got the error message 'Array 1Error in (1 - h) * qs[i] : non-numeric argument to binary operator'.

It is hard to diagnose problems without knowing what commands you have used. Usually people who ask for help have to show all the R code that they have run up to the point that they get the first error message. People also need to show us the output from typing

sessionInfo()

at the R prompt. Have a read of the posting guide at https://www.bioconductor.org/help/support/posting-guide/

I am going to take a wild stab and guess that you have not read the microarray expression data into R correctly. I am guessing that you have actually read a data.frame with columns of different types, and that you have tried to background correct data that is not actually numerical. I am guessing that you may have input a character matrix or some such to backgroundCorrect() instead of a numeric matrix. This could happen if you have read annotation and intensities as a data.frame, then left a column of annotation in when you intended to subset the intensities.

It would be a good idea for you to run your R code line by line and check that you are in fact inputting the correction quantities to each function.

ADD COMMENT
0
Entering edit mode
svlachavas ▴ 840
@svlachavas-7225
Last seen 13 months ago
Germany/Heidelberg/German Cancer Resear…

Dear d.leung,

what kind of platform do you want to perform background ? Also could you provide your current sessionInfo() ? It seems weird to say that you have loaded the limma package and then the above error 

Best,

Efstathios

ADD COMMENT
0
Entering edit mode

Dear Efstathios,

I have made some progress with this

I typed in the command library (limma)

and this error message comes out

Array 1Error in (1 - h) * qs[i] : non-numeric argument to binary operator

How do I fix this problem?

THanks

Dominic

ADD REPLY
0
Entering edit mode

You mean the error appeard after the function library(limma) ?

Or when you used the fuction backgroundCorrect ?

ADD REPLY
0
Entering edit mode

When I use the function backgroundCorrect

ADD REPLY
0
Entering edit mode

Well if you type ?backgroundCorrect:

backgroundCorrect(RG, method="auto", offset=0, printer=RG$printer, normexp.method="saddle", verbose=TRUE)

RG : a numeric matrix, EListRaw or RGList object.

That's why i asked you earlier on which platform do you want to perform background correction ? also, what is your input ? a microarray expression set ? a matrix of expression intensities, etc.? 

ADD REPLY
0
Entering edit mode

I used read.table to read in a .txt file(tab delimited txt file)  The text file contains column with ID and intensities and other columns containing string data.

I suspect it is the way that I put the command in that is the problem

data = read.table("datafile.txt", header=TRUE, sep="\t", as.is=TRUE, quote="", fill=TRUE, check.names=FALSE)

 

ADD REPLY

Login before adding your answer.

Traffic: 699 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