[BioC] get parameters of estimateLogicle

Samuel GRANJEAUD IR/INSERM samuel.granjeaud at inserm.fr
Wed Feb 26 21:01:56 CET 2014


Perfect as usual.

Thanks,
Samuel

> object = lgcl
> sapply(object at transforms, function(transMap) summary(new("transform", 
> .Data = transMap at f))  ,simplify = FALSE)
$`FL3-H`
$`FL3-H`$a
[1] 0

$`FL3-H`$k
transform object 'FL3-H_logicleTransform'

$`FL3-H`$m
[1] 4.5

$`FL3-H`$t
[1] 10000

$`FL3-H`$transformationId
[1] "FL3-H_logicleTransform"

$`FL3-H`$w
[1] 0


$`FL4-H`
$`FL4-H`$a
[1] 0

$`FL4-H`$k
transform object 'FL4-H_logicleTransform'

$`FL4-H`$m
[1] 4.5

$`FL4-H`$t
[1] 10000

$`FL4-H`$transformationId
[1] "FL4-H_logicleTransform"

$`FL4-H`$w
[1] 0




On 26-02-2014 20:07, Mike wrote:
> There isn't such function available for transformList. but you can
> write something like this:
> 
>  setMethod("summary",
>  signature=signature(object="transformList"),
>  definition=function(object, ...)
>  {
>  sapply(object at transforms, function(transMap)
>  summary(new("transform", .Data = transMap at f))
>  ,simplify = FALSE)
>  })
> 
> and then summary(lgcl)
> 
> Mike
> 
> On 02/26/2014 03:00 AM, bioconductor-request at r-project.org wrote:
> 
>> Hi,
>> 
>> I would to get the parameters a, w, m, t that the estimateLogicle
>> computes. The question may be extended to how to get parameters from
>> a transformList or a transformMap. What I would like, is a function
>> that returns the same information as summary applied to a user
>> defined transform.
>> 
>> Cheers.
>> 
>>>> data(GvHD)
>>>> samp <- GvHD[[1]]
>>>> lgcl <- estimateLogicle(samp, channels = c("FL3-H", "FL4-H"))
>>>> lgcl <- logicleTransform(w = 0.74, t=2^18, m =4.5)
>>>> summary(lgcl)
>> 
>> $a
>> [1] 0
>> 
>> $k
>> transform object 'defaultLogicleTransform'
>> 
>> $m
>> [1] 4.5
>> 
>> $t
>> [1] 262144
>> 
>> $transformationId
>> [1] "defaultLogicleTransform"
>> 
>> $w
>> [1] 0.74



More information about the Bioconductor mailing list