[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