[R] probability from different values

Matthias Kohl Matthias.Kohl at stamats.de
Thu Dec 27 15:19:29 CET 2007


one addition ...
> Hello Ricardo,
>
> another solution could be using package distr:
>
> library(distr)
> A <- c(18,18,18,19,20,21,22,23,24,25,26,27,28)
> DA <- DiscreteDistribution(A)
> # maybe
> # support(DA)
> # plot (DA)
>
> B <- c(82,83,84,85,85,86,87,88,89,90,91,91,92)
> DB <- DiscreteDistribution(B)
> # support(DB)
> # plot(DB)
>
> DC <- DB - DA# convolution with of DB with (-DA)
> # support(DC)
> # plot(DC)
>   
to get the probabilities for the C-values:
p(DC)(support(DC))

> hth,
> Matthias
>
> Richard.Cotton at hsl.gov.uk wrote:
>   
>>> i have two vectors (A and B) and i need create another vector (C) 
>>> from the subtraction of A's values with the B's values. How can i 
>>> estimate the probability of C's values if i have differents values 
>>> combinations of A and B that can result in the same value? like 90 
>>> -20 = 70 and 91 - 21 = 70.
>>>
>>> example:
>>> B = {18      18    18     19      20      21      22      23      24
>>> 25      26      27   28....} 
>>> A = {82     83     84     85      85     86       87     88       89
>>> 90      91     91    92 ....}
>>>     
>>>       
>> Try this:
>> A = c(18,18,18,19,20,21,22,23,24,25,26,27,28)
>> B = c(82,83,84,85,85,86,87,88,89,90,91,91,92)
>> lenA = length(A)
>> lenB = length(B)
>> AA = rep(A, lenB)
>> BB = rep(B, each=lenA)
>> table(AA-BB)/(lenA*lenB)
>>
>> Regards,
>> Richie.
>>
>> Mathematical Sciences Unit
>> HSL
>>
>>
>> ------------------------------------------------------------------------
>> ATTENTION:
>>
>> This message contains privileged and confidential inform...{{dropped:20}}
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>     
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list