[R] from a normal bivariate distribution to the marginal one

From: Elisabetta Petracci <elisabetta.petracci_at_unimi.it>
Date: Wed, 23 Jan 2008 13:18:50 +0100


Hello,

I'm quite new with R and so I would like to know if there is a command to calculate an integral.
In particular I simulated a bivariate normal distribution using these simple lines:

rbivnorm <- function(n, # sample size

                    mux,     # expected value of x

                    muy,     # expected value of Y

                    sigmax,  # standard deviation of X

                    sigmay,  # standard deviation of Y

                    rho){    # correlation coefficient

       
           x <- rnorm(n,mux,sigmax)
           y <- rnorm(n,muy+rho*sigmay*(x-mux)/sigmax, sigmay*sqrt(1-rho2))

           cbind(x,y)
           }

In this way I've sampled from a normal bivariate distribution and I've obtained two vectors of values:x and y..

Now, I would like to find the marginal distribution of y from the bivariate distribution.
Is there a command to do this?
Thank you,

Elisabetta



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. Received on Wed 23 Jan 2008 - 12:19:22 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 23 Jan 2008 - 19:30:18 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive