Entering edit mode

lohith madireddy
▴
40

@lohith-madireddy-3767
Last seen 9.5 years ago

Dear Audrey,
the labeling works good now. Thanks for the new version..
There are certainly memory issues performing spatial plots. I just
have 69 arrays with around 9000 features on each array. I can imagine
its impossible to perform these plots for large number of arrays of
larger sizes.The software mostly never produces spatial plots since it
can not allocate some 'xyz' MB. The error is as show below: I have a
max of 4095MB allocated to R on a 64 bit system.
Error : cannot allocate vector of size 150.1 Mb
In addition: There were 50 or more warnings (use warnings() to see the
first 50)
In addition: Warning messages:
1: In arrayQualityMetrics(ratDTH.ab, do.logtransform = TRUE, outdir =
"QA_BeforeRMA_69Chips2", ?:
?Cannot draw spatial distribution of intensities
2: In rep.int(1L:ncol(x), rep.int(nrow(x), ncol(x))) :
?Reached total allocation of 3583Mb: see help(memory.size)
3: In rep.int(1L:ncol(x), rep.int(nrow(x), ncol(x))) :
?Reached total allocation of 3583Mb: see help(memory.size)
4: In inherits(x, "data.frame") :
?Reached total allocation of 3583Mb: see help(memory.size)
5: In inherits(x, "data.frame") :
?Reached total allocation of 3583Mb: see help(memory.size)
So, I tried to do spatial plots and MA plots separately creating 'aqm'
objects and generating report for these two separately (performing
garbage collection in between). I was successful with all except the
MAplots. It gives the following error. I can not exactly see what's
wrong with the error it gave. Are there any NULL values?? I do not
think so. There were no errors when I used the 'arrayQualityMetrics'
wrapper function though. Also, as a suggestion it would be very useful
also to have a loess curve in each of the MAplots as we get using the
affy or affyPLM packages.
> maplot.aqm<-aqm.maplot(ratDTH.ab)
Error in limits.and.aspect(prepanel.default.xyplot, prepanel =
prepanel, ?:
?need at least one panel
In addition: Warning messages:
1: In is.na(x) : is.na() applied to non-(list or vector) of type
'NULL'
2: In is.na(x) : is.na() applied to non-(list or vector) of type
'NULL'
One more request is with the Wrapper function. I felt the option of
customizing the report is a good thing as the user can plot only what
he needs. But, it would be easier if the wrapper function itself has
an option(s) to not plot some of the plots the user is not interested
in. In my case, I would not want the 'spatial distribution of feature
intensities' plot as it almost always fails to draw all or few of the
69 arrays and I would prefer to do it separately.
Awaiting your reply,
Thank you
Lohith
On Fri, Oct 30, 2009 at 1:51 AM, Audrey Kauffmann <audrey at="" ebi.ac.uk=""> wrote:
>
> Dear Jeremy,
>
> Sorry I did not think about this, I think it takes a couple of days
between the moment the developer modifies the version of the package
and the moment you can actually install it using biocLite. You should
be able to install version 2.4.1 after the week-end.
>
> Best wishes,
> Audrey
>
> jeremy wilson wrote:
>>
>> Dear Audrey,
>>
>> thanks for your reply. I tried re-installing the package from bioc.
It gave the same results. I later installed 2.5 version of the
software. This also gave the same plots. There is still the same
problem with ordering of colors.
>> I see you asked me to use version 2.4.1 but I could not find it in
the release or development version. Am I looking at the wrong places?
Please let me know and thanks for your time...
>>
>> Thank you
>> Lohith
>>
>> P.S: I know this might be just a simple bug in the program and I
wanted to see if I can successfully tweak the R script for the
software but I did not find the R script in both the release and
development versions..It would be great to see the code..
>>
>>
>> On Thu, Oct 29, 2009 at 2:41 AM, Audrey Kauffmann <audrey at="" ebi.ac.uk="" <mailto:audrey="" at="" ebi.ac.uk="">> wrote:
>>
>> ? ?Dear Jeremy,
>>
>> ? ?Thank you for reporting this. You were right there was a problem
>> ? ?in the code for the ordering of the colours.
>> ? ?I have now corrected it, so arrayQualityMetrics version 2.4.1
>> ? ?should be working fine.
>> ? ?Please let me know if you keep experiencing problems.
>>
>> ? ?Best wishes,
>> ? ?Audrey
>>
>> ? ?jeremy wilson wrote:
>>
>> ? ? ? ?Dear BioConductors,
>>
>> ? ? ? ?I am using arrayQualityMetrics software package to do the
QA.
>> ? ? ? ?I very much
>> ? ? ? ?like all the QA plots it produces including the latest PCA
>> ? ? ? ?plot in one
>> ? ? ? ?report. I am using the latest version of R(2.10) and BioC
>> ? ? ? ?(2.5). I am facing
>> ? ? ? ?a problem with inappropriate labeling of covariates. Please
>> ? ? ? ?see the attached
>> ? ? ? ?plots.
>> ? ? ? ?I have total of 69 samples. I asked the software to label
the
>> ? ? ? ?plots based on
>> ? ? ? ?a covariate, the summary of which is as below:
>>
>> ? ? ? ?Notation ? ?total number of samples in order
>> ? ? ? ?V ? ? ? ? ? ? ? ? ? ? ? ? 23
>> ? ? ? ?G ? ? ? ? ? ? ? ? ? ? ? ? 38
>> ? ? ? ?N ? ? ? ? ? ? ? ? ? ? ? ? ?8
>>
>> ? ? ? ?The problem is that it does not label the graphs correctly
in
>> ? ? ? ?the plots!! In
>> ? ? ? ?the plots, the first set of samples correspond to V, the
next
>> ? ? ? ?38 samples for
>> ? ? ? ?G and the last 8 for N. In the boxplot for example, the G
was
>> ? ? ? ?labelled RED
>> ? ? ? ?and this should correspond to 38 samples but in the plot,
the
>> ? ? ? ?boxplots
>> ? ? ? ?corresponding to these 38 samples were colored GREEN. The
ones
>> ? ? ? ?that were
>> ? ? ? ?labeled RED correspond to 8 'N' samples which should be
>> ? ? ? ?represented in blue
>> ? ? ? ?according to the label at the top of plot. ?Likewise for the
>> ? ? ? ?other two. I
>> ? ? ? ?think it can be better understood at the first glance of the
>> ? ? ? ?plots.
>>
>> ? ? ? ?The phenotype of the samples was given correctly.
>> ? ? ? ?Following is the code I used for generating the aqm report.
>>
>> ? ? ? ?library(arrayQualityMetrics)
>> ? ? ? ?arrayQualityMetrics(rat.ab, do.logtransform=TRUE,
outdir="QA",
>> ? ? ? ?intgroup="Material", grouprep=TRUE)
>>
>> ? ? ? ?I have attached only a few plots. I can produce all of them
if
>> ? ? ? ?it helps...
>> ? ? ? ?It might take some time for the post to reach the mailing
list
>> ? ? ? ?as the
>> ? ? ? ?earlier one I posted with a few more plots was held for
>> ? ? ? ?moderator's approval
>> ? ? ? ?as the size of the message body is over the limit!
>>
>> ? ? ? ?Please let me know what is wrong.
>> ? ? ? ?Thanks in advance!
>> ? ? ? ? ? ? ? ?----------------------------------------------------
--------------------
>>
>>
>> ? ? ? ?------------------------------------------------------------
------------
>>
>> ? ? ? ?------------------------------------------------------------
------------
>>
>>
>>
>> ? ? ? ?_______________________________________________
>> ? ? ? ?Bioconductor mailing list
>> ? ? ? ?Bioconductor at stat.math.ethz.ch
>> ? ? ? ?<mailto:bioconductor at="" stat.math.ethz.ch="">
>> ? ? ? ?https://stat.ethz.ch/mailman/listinfo/bioconductor
>> ? ? ? ?Search the archives:
>> ? ? ?
?http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>>
>> ? ?-- ? ?Audrey Kauffmann
>> ? ?EMBL - EBI
>> ? ?Cambridge UK
>> ? ?+44 (0) 1223 492 631 http://www.ebi.ac.uk/~audrey
>> ? ?<http: www.ebi.ac.uk="" %7eaudrey="">
>>
>>
>
> --
> Audrey Kauffmann
> EMBL - EBI
> Cambridge UK
> +44 (0) 1223 492 631 http://www.ebi.ac.uk/~audrey
>