[R] probability from different values

Matthias Kohl Matthias.Kohl at stamats.de
Thu Dec 27 15:16:39 CET 2007


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)

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.
>



More information about the R-help mailing list