[Bioc-devel] pval in mas5calls()

cstrato cstrato at aon.at
Mon May 28 21:33:39 CEST 2007


Dear Martin

Thank you for your immediate reply, it works great.
(Maybe this could be mentioned in the help file for ?mas5calls)

Best regards
Christian

Martin Morgan wrote:
> Christian:
>
> The values are there, just not accessible with se.exprs. Try
> 'assayData(obj)[["se.exprs"]]':
>
>   
>> example(mas5calls)
>>     
>
> ms5cll> data(affybatch.example)
>
> ms5cll> PACalls <- mas5calls(affybatch.example)
> Getting probe level data...
> Computing p-values
> Making P/M/A Calls
>   
>> dim(assayData(PACalls)[["se.exprs"]])
>>     
> [1] 150   3
>
> Martin
>
> cstrato <cstrato at aon.at> writes:
>
>   
>> Dear all,
>>
>> Could it be that function mas5calls() no longer returns the p-values, 
>> since se.exprs does no longer exist?
>>
>> I get:
>> pr.call <- mas5calls(pr)
>> call <-exprs(pr.call)
>> pval <-se.exprs(pr.call) ||
>> Error in function (classes, fdef, mtable)  :
>>         unable to find an inherited method for function "se.exprs", for 
>> signature "ExpressionSet"
>>
>> Interestingly, the source code in "mas5.R" still creates an 
>> ExpressionSet with "se.exprs=p".
>>
>> sessionInfo()
>> R version 2.5.0 (2007-04-23)
>> i386-apple-darwin8.9.1
>>
>> locale:
>> C
>>
>> attached base packages:
>> [1] "tools"     "stats"     "graphics"  "grDevices" "utils"     "datasets"
>> [7] "methods"   "base"    
>>
>> other attached packages:
>> hgu95av2cdf        affy      affyio     Biobase
>>    "1.16.0"    "1.14.0"     "1.4.0"    "1.14.0"
>>
>> Best regards
>> Christian
>>
>> _______________________________________________
>> Bioc-devel at stat.math.ethz.ch mailing list
>> https://stat.ethz.ch/mailman/listinfo/bioc-devel
>>     
>
>



More information about the Bioc-devel mailing list