Re: [aroma.affymetrix] mufColumns on genes with more than 500 probes across the gene

2015-06-05 Thread Mark Robinson
. To unsubscribe from this group and stop receiving emails from it, send an email to aroma-affymetrix+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout. -- Prof. Dr. Mark Robinson Statistical Bioinformatics Group, UZH http://ow.ly/riRea -- -- When

Re: [aroma.affymetrix] which version cdf should used when apply FIRMAGene

2013-01-06 Thread Mark Robinson
Hi Zaiwei, On Sun, Dec 16, 2012 at 2:22 AM, zhouzaiwei zhouzai...@163.com wrote: Hi , I want to apply FIRMAGene to analysis differential splicing events of hugene_1.0_st array and have red article (Robinson Speed, 2007) and

[aroma.affymetrix] Re: How to extract raw probe intensity from .CEL file

2011-08-09 Thread Mark Robinson
Perhaps something like this is what you want (note: different chip to what you are using)? df - readDataFrame(getCdf(cs), verbose=-80) [...snip...] head(df) unit unitName unitType unitDirection unitNbrOfAtoms group groupName 11 7892501 expression sense 4 1

Re: [aroma.affymetrix] Re: How to extract raw probe intensity from .CEL file

2011-08-09 Thread Mark Robinson
of intensities, but I don't know which probe ids they correspond to. Hope this clarifies. Thanks On Aug 9, 6:03 am, Mark Robinson markrobinson@gmail.com wrote: Perhaps something like this is what you want (note: different chip to what you are using)? df - readDataFrame(getCdf(cs), verbose=-80

[aroma.affymetrix] Re: Question about Firma

2011-06-26 Thread Mark Robinson
--- Dr Florence Jaffrezic INRA Bat 211 78352 Jouy-en-Josas Cedex France Tel: (+33) 1 34 65 21 94 Fax: (+33) 1 34 65 22 10 --- Firma_Exon_question.r -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

Re: [aroma.affymetrix] Re: gcRMA for Gene ST Arrays

2011-06-22 Thread Mark Robinson
Mark Robinson Sun, 29 Mar 2009 23:12:09 -0700 Hi Mario. I have made some modifications to the reading of probe_tab files and to the computing of affinities so that this procedure can run now, either as you have done below by choosing lowly expressed probes, or (perhaps preferably

Re: [aroma.affymetrix] Extract raw data for core transcripts from EXON arrays

2011-02-10 Thread Mark Robinson
/ -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: mrobin...@wehi.edu.au e: m.robin...@garvan.org.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] exon array analysis

2011-01-23 Thread Mark Robinson
aroma.affymetrix group with website http://www.aroma-project.org/. To post to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe and other options, go to http://www.aroma-project.org/forum/ -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

Re: [aroma.affymetrix] Firma score help

2011-01-14 Thread Mark Robinson
aroma.affymetrix group with website http://www.aroma-project.org/. To post to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe and other options, go to http://www.aroma-project.org/forum/ -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

Re: [aroma.affymetrix] Installing FIRMAGene

2010-10-27 Thread Mark Robinson
to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe and other options, go to http://www.aroma-project.org/forum/ -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: mrobin...@wehi.edu.au e: m.robin

Re: [aroma.affymetrix] Missing file bpmapCluster2Cdf.R in Create a CDF from a BpMap

2010-10-20 Thread Mark Robinson
Apologies folks. Inadvertent file permissions change. It should work now. Cheers, Mark -- Forwarded message -- From: vegard vegard.nyga...@medisin.uio.no Date: Thu, Oct 21, 2010 at 2:55 AM Subject: [aroma.affymetrix] Missing file bpmapCluster2Cdf.R in Create a CDF from a BpMap

[aroma.affymetrix] Re: run FirmaGene on exon array ?

2010-08-29 Thread Mark Robinson
? Thanks, Qicheng -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: sorry, I can not reproduce the table 2 in FIRMAGene paper

2010-08-09 Thread Mark Robinson
,TisMix_mix9,28.2833712932546,ACTA1 7982090,TisMix_mix1,27.4880562975114,SNORD115-42 8103789,TisMix_mix1,27.399249112706,GPM6A 7982008,TisMix_mix1,25.0039468447955,SNORD115-1 Could you please tell me where I am wrong ? Thanks, Qichengm -- Mark Robinson, PhD (Melb

Re: [aroma.affymetrix] Firma Scores in the range of 700-1500

2010-07-22 Thread Mark Robinson
are subscribed to the Google Groups aroma.affymetrix group with website http://www.aroma-project.org/. To post to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe and other options, go to http://www.aroma-project.org/forum/ -- Mark Robinson, PhD (Melb

Re: [aroma.affymetrix] FIRMA and GenomeGraphs

2010-05-13 Thread Mark Robinson
and other options, go to http://www.aroma-project.org/forum/ -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] Re: what kind of file

2010-04-01 Thread Mark Robinson
, comment.char=,stringsAsFactors=FALSE) Error in read.table(MoGene-1_0-st-v1.probe.tab, sep = \t, header = TRUE, : no lines available in input I double checked the file I downloaded from Affy and it seems to be fine! Thx, Daniela On Mar 16, 4:56 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi

Re: [aroma.affymetrix] Re: FIRMAGene command

2010-03-18 Thread Mark Robinson
into a txt file? Many thanks for your help. I really appreciate it. Daniela On Mar 16, 4:57 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Daniela. You haven't told us what inputs you've used for 'plm' and 'cls' ... and what is stored in 'u'? Have you read the docs

Re: [aroma.affymetrix] getUniqueCdf inflates dimensions of original cdf

2010-03-16 Thread Mark Robinson
unique CDF to the original CDF, which obviously is curious in itself. So, a full explanation of what you've done upstream of this would be useful. Cheers, Mark On 16-Mar-10, at 1:09 AM, Henrik Bengtsson wrote: Hi, I leave this one to Mark Robinson who is designed createUniqueCdf

Re: [aroma.affymetrix] FIRMAGene command

2010-03-16 Thread Mark Robinson
@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan

Re: [aroma.affymetrix] Is the PdInfo2Cdf.R script working

2010-03-16 Thread Mark Robinson
Hi Peter. I haven't looked at this since early 2009 and our motive in making this available as a script (instead of within the aroma.affymetrix package) was simply as an FYI and the code is readily available for others to modify for their needs. You will of course understand that underlying

Re: [aroma.affymetrix] can't load CDF file

2010-03-09 Thread Mark Robinson
Hi Daniela. Is your CDF in the: annotationData/chipTypes/MoEx-1_0-st-v1/ directory? (http://aroma-project.org/node/66) Cheers, Mark Hi, I stored my CDF file in annotationData/chipTypes; nevertheless I cannot upload the file. Can anyone please tel me what I am doing wrong:

Re: [aroma.affymetrix] custom CDF and flat file problems

2010-01-17 Thread Mark Robinson
email to aroma-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics

[aroma.affymetrix] Re: Gene-Level Summarization of Expression Data

2010-01-17 Thread Mark Robinson
--is it the right one for this? -Randy On Jun 19 2009, 10:08 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Steve. I don't know how common this is. Basically, a colleague found agene that was very differentially expressed when analyzing using the Affymetrix probesets definition and found

Re: [aroma.affymetrix] Re: a question about FIRMA score

2010-01-10 Thread Mark Robinson
at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] FIRMA SCORE from different test

2010-01-10 Thread Mark Robinson
-- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] Question for custom CDF of ST-Array

2010-01-10 Thread Mark Robinson
group. To post to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark

Re: [aroma.affymetrix] questiones about annotations for exon arrays: no gene symbol or refseq id for majority of core probe sets?

2010-01-10 Thread Mark Robinson
://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] Re: FIRMA score for each transcript

2009-12-16 Thread Mark Robinson
? Thanks! /Libing On Fri, Nov 6, 2009 at 6:35 AM, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Libing. Are you after the probe IDs from the probe.tab file? For example: Probe IDProbe Set IDprobe x probe y assembly seqname start stop strand probe sequence

Re: [aroma.affymetrix] probe_id in cdf file

2009-12-16 Thread Mark Robinson
-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

[aroma.affymetrix] Re: annotation of ST gene Arrays

2009-12-16 Thread Mark Robinson
dealing with ST arrays. Thanks, Wade On Feb 10, 3:13 am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Simon. See comments below. I am using the mouse gene ST arrays and am having problems with annotation. When i write a csv file, theannotationis only the probeset_id, no gene names

Re: [aroma.affymetrix] Re: custom CDF

2009-12-11 Thread Mark Robinson
Hi Zaid. There is an example flat file (well, 10 lines of it) at that page that Pierre suggested. You'll want to filter the missing lines and make sure all the data in a column is of the same type. Cheers, Mark No. My flat file has some lines with missing values and non integer values. Is

Re: [aroma.affymetrix] Problem generating CDF file (Arabidopsis)

2009-12-06 Thread Mark Robinson
not allow it? Thanks, Renyi On Sat, Dec 5, 2009 at 2:43 PM, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Renyi. One thing to do is check that the genome positions in your BPMAP file for chr5 are all 0. To do this, try: library(affxparser) bp - readBpmap(At35b_MR_v04

[aroma.affymetrix] Re: more question on tiling array

2009-12-06 Thread Mark Robinson
of California, Riverside Riverside, CA 92521 Email: renyi@ucr.edu Phone: (951)827-3987 Fax: (951)827-4437 -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3

Re: [aroma.affymetrix] Debate sobre gene-1-0-st-array-analysis

2009-12-04 Thread Mark Robinson
-- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 -- -- When reporting problems on aroma.affymetrix, make sure

Re: [aroma.affymetrix] Re: exon array analysis errors

2009-12-04 Thread Mark Robinson
$`20091119_UHR_Exon2`$out numeric(0) $`20091119_UHR_Exon3` $`20091119_UHR_Exon3`$stats [1] 0 0 0 0 0 $`20091119_UHR_Exon3`$n [1] 18705 $`20091119_UHR_Exon3`$conf [1] 0 0 $`20091119_UHR_Exon3`$out numeric(0) attr(,type) [1] NUSE On Nov 30, 1:28 am, Mark Robinson mrobin...@wehi.edu.au

Re: [aroma.affymetrix] Re: apt- affymetrix power tool

2009-11-30 Thread Mark Robinson
@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics

Re: [aroma.affymetrix] Re: exon array analysis errors

2009-11-30 Thread Mark Robinson
: plotBoxplot.ChipEffectSet(ces, type = RLE, ...) 3: plotBoxplot(ces, type = RLE, ...) 2: plotRle.QualityAssessmentModel(qamTr) 1: plotRle(qamTr) On Nov 24, 1:50 am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Yu Chuan. Comments below. On 24-Nov-09, at 1:00 PM, Yu Chuan wrote: Hi, I

Re: [aroma.affymetrix] exon array analysis errors

2009-11-24 Thread Mark Robinson
-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

Re: [aroma.affymetrix] how to analysis the FIRMA score

2009-11-24 Thread Mark Robinson
to this group, send email to aroma-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb

extracting FIRMA scores (Was: [aroma.affymetrix] a question)

2009-11-17 Thread Mark Robinson
need replicates. Therefore, it does not need this information. Cheers, Mark -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: [aroma.affymetrix] FIRMAGene with masked CEL files

2009-11-16 Thread Mark Robinson
options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Discussion on using-the-genomegraphs-package-with-firma

2009-10-26 Thread Mark Robinson
wrote: Thanks for sharing your code. I don't really understand - when getting the PLM for this instruction are you supposed to set mergeGroups to true like this: plm - ExonRmaPlm(celSet, mergeGroups=TRUE) -L -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

[aroma.affymetrix] Re: error in exon array analysis: fit(plmTr, verbose=verbose)

2009-10-26 Thread Mark Robinson
R.huge_0.2.0 aroma.core_1.2.0 aroma.light_1.13.6 [8] matrixStats_0.1.6 R.rsp_0.3.6R.filesets_0.5.3 digest_0.4.1 R.cache_0.2.0 R.utils_1.2.2 R.oo_1.6.2 [15] R.methodsS3_1.0.3 projectManager_1.0 XML_2.6-0 Mark Robinson wrote

[aroma.affymetrix] Re: selection of CDF and PA call for exon array

2009-10-20 Thread Mark Robinson
-- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 -- --~--~-~--~~~---~--~~ When

[aroma.affymetrix] Re: FIRMA score and limma analysis

2009-10-18 Thread Mark Robinson
0.4552856 20.4449589 3 -0.1899892 4 -0.6681926 50.3545261 6 -0.5645628 Thanks Hailei -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345

[aroma.affymetrix] Re: error in exon array analysis: fit(plmTr, verbose=verbose)

2009-10-15 Thread Mark Robinson
Hi Hailei. For starters, can you give the *full* output of your sessionInfo()? The error you are getting has something to do with the 'preprocessCore' package and I first want to check whether it is a package version mismatch error. I haven't used aroma.affymetrix on R 2.10 and I don't know if

[aroma.affymetrix] Re: error in exon analysis

2009-10-15 Thread Mark Robinson
AromaPlatformInterface AromaMicroarrayDataSet [7] GenericDataFileSet Object Thanks. Hailei -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628

[aroma.affymetrix] Re: Exception: Unknown arguments: cdf, checkChipType

2009-10-14 Thread Mark Robinson
(staticMethod, args = args) -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: AffymetrixCelSet, Could not locate a file for this chip type

2009-10-11 Thread Mark Robinson
Hi Carol. This throws an error because your CDF has a tag (r3) and there is no way to send a tag for the CDF file through the AffymetrixCelSet$byName function. I recommend that you just do this in 2 commands, as you've done successfully. Alternatively, you can remove the ,r3 from your CDF file

[aroma.affymetrix] Re: exon array analysis

2009-10-05 Thread Mark Robinson
,] 8.038540 3.542727 3.808023 7.923751 [8,] 13.260725 6.238900 4.660217 4.200828 I have not sure what I'm doing wrong and would really appreciate any help. Thank you very much in advance, Enid -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan

[aroma.affymetrix] Re: QC of probe level model.

2009-10-03 Thread Mark Robinson
-- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 -- --~--~-~--~~~---~--~~ When reporting problems

[aroma.affymetrix] Re: Mat implementation - comparing MAT (pure) vs MAT aroma.affymetrix

2009-09-01 Thread Mark Robinson
Lavinia. On Sep 1, 1:24 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Lavinia. Yes. Bandwidth(MAT)=probeWindow(aroma.affymetrix MAT) Also, MinProbe(MAT)=nProbes(aroma.affymetrix MAT) Cheers, Mark On 1-Sep-09, at 11:31 AM, Lavinia Gordon wrote: Hi, I have some older MAT pure

[aroma.affymetrix] Re: Mat implementation - comparing MAT (pure) vs MAT aroma.affymetrix

2009-08-31 Thread Mark Robinson
. Can I just check, from the MAT .tag (parameters file), does probeWindow correspond directly to Bandwidth? many thanks Lavinia. -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin

[aroma.affymetrix] Re: Discussion on gene-1-0-st-array-analysis

2009-08-30 Thread Mark Robinson
to proceed after this step Diya Love Cricket? Check out live scores, photos, video highlights and more. Click here -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au

[aroma.affymetrix] Re: installation fails

2009-08-13 Thread Mark Robinson
by: source(http://www.braju.com/R/hbLite.R;); hbLite(R.filesets); library(R.filesets); Q1) Does this work? hbLite(aroma.core); library(aroma.core); Q2) Does this work? hbLite(aroma.affymetrix); library(aroma.affymetrix); Q3) Does this work? Mark (Robinson), you mentioned this problem

[aroma.affymetrix] Re: Re-run aroma.affymetrix

2009-08-09 Thread Mark Robinson
(in plmData, probeData, and reports) that are created by aroma.affymetrix? Thanks in advance. Kind regards, Anbarasu -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin

[aroma.affymetrix] Re: RMA of aroma.affymetrix, affyPLM and affy

2009-08-09 Thread Mark Robinson
Hi Yiwen. I think this thread will answer your question in detail: http://groups.google.com/group/aroma-affymetrix/browse_thread/thread/1b0ab11fad9b4df3 In brief, the main difference is how the probe-level linear model is fit -- median polish OR iteratively reweighted least squares with a

[aroma.affymetrix] Re: Re-run aroma.affymetrix

2009-08-06 Thread Mark Robinson
, and reports) that are created by aroma.affymetrix? Thanks in advance. Kind regards, Anbarasu -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61

[aroma.affymetrix] Re: MedianNormalization

2009-07-31 Thread Mark Robinson
this with extract(). For example, say 'cs' is an AffymetrixCelSet object. You can extract the first 2 samples with: csSubset - extract(cs, 1:2) Hope that helps. Mark Thank you very much. -- Cathrine Mitchell -- Mark Robinson, PhD (Melb) Epigenetics Laboratory

[aroma.affymetrix] Re: Differential expression analysis

2009-07-31 Thread Mark Robinson
-- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: CDF creation

2009-07-24 Thread Mark Robinson
Hi Naresh. That particular script (flat2Cdf.R) creates a PM-only CDF. But, if you understand the format of a CDF file -- you might read in the your original CDF with readCdf() from the 'affxparser' package to understand it -- then you should be able to modify the script to include MM probes.

[aroma.affymetrix] Re: Discussion on using-the-genomegraphs-package-with-firma

2009-07-16 Thread Mark Robinson
by using windows, you could alternatively install the cygwin tools and then have access to the standard battery of command line tools, such as grep/awk from the example. Cheers, Mark regards, John -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan

[aroma.affymetrix] Re: Calling enriched regions in tiling array experiments

2009-07-06 Thread Mark Robinson
Hi Lars. Apologies for the slow response. I do have some scripts for calling enriched regions, but they are not really ready for public consumption. And, they are geared more towards ChIP-chip of modified histones / DNA methylation than either TF binding or transcript expression. We have

[aroma.affymetrix] Re: Gene-Level Summarization of Expression Data

2009-06-20 Thread Mark Robinson
. Regards, -Steve On Jun 17, 9:56 am, Steve Piccolo steve.picc...@gmail.com wrote: Yesterday I posted this question to the list, but the spam blocker didn't let it through. Below my question is a response from Mark Robinson

[aroma.affymetrix] Re: FIRMA score for each transcript

2009-06-17 Thread Mark Robinson
, Apr 11, 2009 at 5:48 PM, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Libing. As the error message suggests, there are no degrees of freedom for the fit, meaning you have no replicates. It appears you only have 2 total samples, one for each group. You wouldn't be able to use limma to do

[aroma.affymetrix] Re: FIRMAGene

2009-06-11 Thread Mark Robinson
, Nick -- Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: SNPs affecting EXon splicing detection

2009-06-11 Thread Mark Robinson
, 12:32 am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Sabrina. How about you try and create a 'flat' file like the one described at:http://groups.google.com/group/aroma-affymetrix/web/creating-cdf-file ... Presumably, you will be comfortable with the Exon Array's 'probetab' file

[aroma.affymetrix] Re: CDF files for Affymetrix whole transcript arrays (Gene 1.0, Exon 1.0)

2009-06-11 Thread Mark Robinson
. Or, if you could direct me to some document that would tell me how to build custom cdfs, that would be good too. Thanks very much, Dick On Aug 4 2008, 10:05 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi folks. I've been updating/creating CDF files for some of the recent Affymetrix whole

[aroma.affymetrix] Re: SNPs affecting EXon splicing detection

2009-06-09 Thread Mark Robinson
. Cheers, Mark Sabrina On Apr 30, 3:46 am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Sabrina. I have not had to deal with this myself, but I do know that it exists and I can at least suggest a possible route to exclude affected exons. Presumably, there is a database (dbSNP

[aroma.affymetrix] Re: SNPs affecting EXon splicing detection

2009-06-06 Thread Mark Robinson
am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Sabrina. I have not had to deal with this myself, but I do know that it exists and I can at least suggest a possible route to exclude affected exons. Presumably, there is a database (dbSNP?) that tells you the genome locations of each SNP

[aroma.affymetrix] Re: FIRMAGene

2009-06-04 Thread Mark Robinson
Hi Ettore. Comments below. I have the following question. In the sup3.R file the probe level model fitting is realised using the instructions: plm - RmaPlm(csNU) fit(plm, verbose=verbose) where csNU is an object obtained after background correction, quantile normalisation and

[aroma.affymetrix] Re: Discussion on affymetrix-defined-transcript-clusters

2009-05-08 Thread Mark Robinson
special mouse Exon array CDF files... any advice on how I can get hold of them? Thanks -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Compare the splicing pattern of two samples

2009-05-06 Thread Mark Robinson
23153762 35 6 2315373 23153772 46 What does each column, especially the unitName and groupName mean ? And how can I correlate unitName and groupName to gene name and exon number? Thanks in advance! Xinjun On Thu, Apr 30, 2009 at 3:31 PM, Mark Robinson

[aroma.affymetrix] Re: Compare the splicing pattern of two samples

2009-04-30 Thread Mark Robinson
row corresponds to. I'm guessing here, but you are probably interested in the YRIvsCEU parameter, so you'd probably sort on the p.value.YRIvsCEU column ... Cheers, Mark On Tue, Apr 28, 2009 at 4:29 PM, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Xinjun. Comments below

[aroma.affymetrix] Re: SNPs affecting EXon splicing detection

2009-04-30 Thread Mark Robinson
with SNP data before, can anyone give me a hint? Thanks a lot! Sabrina -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Compare the splicing pattern of two samples

2009-04-28 Thread Mark Robinson
here. The probesets where the Grp2 coefficient is significantly different from 0 may highlight differentially spliced exons. Does that help? Mark Thanks in advance! Xinjun On Mon, Apr 27, 2009 at 6:19 AM, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Xinjun. Quick comments

[aroma.affymetrix] Re: Quality assessment of Gene ST Array

2009-04-27 Thread Mark Robinson
Mitchell -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Error Occured When Convertting CDF File

2009-04-25 Thread Mark Robinson
: system.time({ writeCdf(outFilename, cdfheader = cdfHeader, cdf = cdfUnits, cdfqc = cdfQcUnits, overwrite = TRUE, verbose = verbose2) }) 1: convertCdf(files[kk], outFiles[kk], version = 4, force = TRUE, verbose = TRUE) On Sat, Apr 25, 2009 at 4:37 PM, Mark

[aroma.affymetrix] Re: question for using custom Affy Exon array CDF

2009-03-31 Thread Mark Robinson
, Mark On 01/04/2009, at 6:25 AM, jing ma wrote: Dear Mark, When doing background subtraction for the Affy exon array, do you think it is sufficient to use the CORE CDF (I usually use HuEx-1_0- st-v2,coreR3,A20071112,EP.cdf)? Thanks, Jing On Mon, Mar 16, 2009 at 5:35 PM, Mark Robinson

[aroma.affymetrix] Re: gcRMA for Gene ST Arrays

2009-03-25 Thread Mark Robinson
where TRUE/FALSE needed -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Write normalized intensities as a CEL file

2009-03-23 Thread Mark Robinson
wondering if writeCel() is the one I should use. There are no documents for its usage. Can anyone provide more details? Or some links. Thanks! -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin

[aroma.affymetrix] Re: Custom CDF Creation

2009-03-19 Thread Mark Robinson
Hi Jake. As a starting point, you might have a look at: http://groups.google.com/group/aroma-affymetrix/web/creating-cdf-files-from-scratch In there is a script called 'flat2Cdf.R' that takes a flat file of probe information (X/Y location on the chip, probe sequence, identifiers) and creates a

[aroma.affymetrix] Re: Problem understanding FIRMA

2009-03-17 Thread Mark Robinson
/1b0ab11fad9b4df3/f745ed0860546313 But, I'd be interested to hear more details if you do look into it more. Cheers, Mark Best regards Christian On Mar 16, 9:13 am, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Christian. From what I can tell looking at your code (rather quickly, i must admit

[aroma.affymetrix] Re: Using GenomeGraphs with FIRMA (Was: FIRMA score)

2009-03-04 Thread Mark Robinson
? GenomeGraphs_ENSG0060237.pdfGenomeGraphs_ENSG0060237.R -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: to use array list file

2009-02-12 Thread Mark Robinson
. The reason is that I want to process these CEL files together, however I don't want to copy them to one folder. Thanks. -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345

[aroma.affymetrix] Re: probes missing from mogene10st.db

2009-02-11 Thread Mark Robinson
this to the BioC mailing list instead? thanks for your help, Sebastien -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: annotation of ST gene Arrays

2009-02-10 Thread Mark Robinson
-- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 -- --~--~-~--~~~---~--~~ When reporting

Ensembl-centric CDFs from scratch (Was: Re: [aroma.affymetrix] Error: cannot allocate vector of size 692.8 Mb)

2009-02-05 Thread Mark Robinson
-- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 -- --~--~-~--~~~---~--~~ When reporting problems

[aroma.affymetrix] Re: changing colours in plotRle

2009-02-03 Thread Mark Robinson
Hi Sebastien. Note that plotRle() eventually makes a call to 'bxp' (in the graphics package that is loaded by default) and any/all arguments are passed on. Have a look at ?bxp for what you can specify. For example: [...] qamTr - QualityAssessmentModel(plmTr) plotRle(qamTr,

[aroma.affymetrix] Re: Error: cannot allocate vector of size 692.8 Mb

2009-01-28 Thread Mark Robinson
() is: 1535.875 Does anyone have any suggestion to solve this problem? Thanks!!! Sabrina -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: exon array with technical replicates

2009-01-23 Thread Mark Robinson
(nProbesPerExon) Thanks and Have a great weekend Sabrina On Jan 22, 5:18 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Sabrina. See below. On 23/01/2009, at 3:17 AM, sabrina wrote: Hi, Mark: Thanks for the suggestions. I think I will go with one of the replicates for now

[aroma.affymetrix] Re: exon array with technical replicates

2009-01-22 Thread Mark Robinson
these rows before you use lmFit(). Cheers, Mark Thanks! Sabrina On Jan 20, 7:48 pm, Mark Robinson mrobin...@wehi.edu.au wrote: Hi Sabrina. Do you have biological replicates of some samples and technical replicates of others? Or, just technical replicates of everything? My experiment

[aroma.affymetrix] Re: Can't read the CDF file

2009-01-21 Thread Mark Robinson
would be greatly appreciated. Thanks. Anbarasu -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: exon array with technical replicates

2009-01-20 Thread Mark Robinson
Hi Sabrina. Do you have biological replicates of some samples and technical replicates of others? Or, just technical replicates of everything? My experiment has two groups, one has 5 samples (biological replicates, that is 5 mice from one strain), and the other has 4 samples. Among 5

[aroma.affymetrix] Re: exon array with technical replicates

2009-01-19 Thread Mark Robinson
! Sabrina -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

Re: Reproducing RMA with Gene ST data (Was: Re: [aroma.affymetrix] Re: How do you analyze Gene ST Data?)

2009-01-12 Thread Mark Robinson
looks good. Pairs plots are attached. Thanks, Peter On Thu, Dec 4, 2008 at 5:41 PM, Mark Robinson mrobin...@wehi.edu.au wrote: Thanks Peter. Perhaps you can repeat this comparison after the next release (this will be very soon!) and split the aroma.affymetrix comparison

[aroma.affymetrix] Re: SNP 6.0 processing

2008-12-16 Thread Mark Robinson
. Thanks, Anguraj -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin...@garvan.org.au e: mrobin...@wehi.edu.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852

[aroma.affymetrix] Re: Batch Effects

2008-12-15 Thread Mark Robinson
get decent normalization between batches. Is there a way to adjust for batch effects in aroma.affy that can be incorporated into an exon array analysis? Thank you, Sarah -- Mark Robinson Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: m.robin

[aroma.affymetrix] Re: FIRMA score

2008-12-10 Thread Mark Robinson
of samples ... Thanks for all these questions Sabrina. This will help make a nice vignette ... when I get time :) Unless you'd be interested in summarizing your approach! :) Cheers, Mark Sabrina On Dec 9, 5:02 pm, Mark Robinson [EMAIL PROTECTED] wrote: Hi Sabrina. Great work! See below

Re: Reproducing RMA with Gene ST data (Was: Re: [aroma.affymetrix] Re: How do you analyze Gene ST Data?)

2008-12-03 Thread Mark Robinson
need to create or download your mogene10stv1cdf library from the Affy unsupported CDF file (https://stat.ethz.ch/pipermail/bioc- devel/2007-October/001403.html has some detail on how to do this). However, as Mark Robinson pointed out there are potential issues with using the Affy unsupported

  1   2   >