[R] p value for joint probability

[moleps]

My terminology is probably way off. I´ll try again in plain english. 

I´d like  to generate a scatter plot of r1 & r2 and color code each pair according to the probability of observing the pair given that the two samples (r1 & r2) are drawn from two independent normal distributions.

rr<-data.frame(r1=-rnorm(1000,10,5),r2=-rnorm(1000,220,5))

with(rr,plot(r1,r2))

Best,
//M




On 31. jan. 2011, at 23.13, Peter Ehlers wrote:

> On 2011-01-31 12:42, moleps wrote:
>> Dear all,
>> 
>> Given
>> 
>> rr<-data.frame(r1<-rnorm(1000,10,5),r2<-rnorm(1000,220,5))
>> 
>> 
>> How can I add a column (rr$p) for the joint probability of each r1&  r2 pair?
> 
> If you take the values in each pair to be observations
> from two independent Normal distributions, it's easy:
> The "joint probability" of those values is zero.
> 
> But I suspect you mean something else by "joint probability".
> Can you elaborate?
> 
> Peter Ehlers
> 
>>  I know how to add the column.. I just dont know how to compute the p value for joint probabilities given the two samples.
>> 
>> 
>> //M
>> 
>> ______________________________________________
>> 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.
>