similar to: bivariate normal cdf and rho

Hi all, I'm looking for an analytic method to obtain the PDF & CDF of the sampling distribution of a given correlation (rho) at a given sample size (N). I've attached code describing a monte carlo method of achieving this, and while it is relatively fast, an analytic solution would obviously be optimal. get.cors <- function(i, x, y, N){ end=i*N

hello, I'm trying to improve the speed of my calculation but didn't get to a satisfying result. It's about the numerical Integration of a bivariate normal distribution. The code I'm currently using x <- qnorm(seq(.Machine$double.xmin,c(1-2*.Machine$double.eps),by=0.01), mean=0,sd=1) rho <- 0.5 integral <- function(rho,x1){

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

Suppose I know the value of cumulative bivariate standard normal distribution. How can I solve correlation between variables? Pekka --------------------------------- [[alternative HTML version deleted]]

-- R Help List -- I am looking for a bivariate normal cdf routine in R. I have some fortran routines for this, which appear to be based on 15-point quadrature. Any guidance/suggestions on making these in loadable R-functions would be appreciated. Thanks, Dan =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Daniel A. Powers, Ph.D. Department of Sociology University of Texas at Austin

R-Help I am attempting to create a series of bivariate normal distributions. So using the mvtnorm library I have created the following code ... # Standard deviations and correlation sig_x <- 1 sig_y <- 1 rho_xy <- 0.0 # Covariance between X and Y sig_xy <- rho_xy * sig_x *sig_y # Covariance matrix Sigma_xy <- matrix(c(sig_x ^ 2, sig_xy, sig_xy, sig_y ^ 2), nrow = 2, ncol = 2)

Hi, Let f(rho) = E[F_1(x) F_2(y)], i.e f(rho) is the expectation of F(x) * F(y) with respect to the bivariate Gaussian density with mean 0 and covariance matrix [1 rho; rho 1]. Moreover, assume F_1(x) and F_2(y) to be increasing functions of x and y respectively. I was wondering if it was true that f(rho) is an increasing function of rho. If so, are there any references? Best, Agos

Dear All I want to generate variable with Bivariate Normal Distribution by use mean1 = a, variance1 = b, mean2 = c, variance2 = d, rho = e. How I can do this. Many Thanks. IRD [[alternative HTML version deleted]]

Dear R-Users I would like to estimate a constrained bivariate normal density, the constraint being that the means are of equal magnitude but of opposite signs. So I need to estimate four parameters: mu (meanvector (mu,-mu)) sigma_1 and sigma_2 (two sd deviations) rho (correlation coefficient) I have looked at several packages, including Gaussian mixture models in Mclust, but I am not sure

Hi, there. Does anyone know the R function for calculating the cdf of bivariate normal distribution function? Thanks. Yulei [[alternative HTML version deleted]]

Dear all, I am calculating the bivariate skew normal cdf in "sn" package using "pmsn" function. Although it is quite convenient ( thanks to prof. Azzalini) but it seems to be slow. For example, it takes about 1 minute in calculation of 100k of such cdf values. I am thinking to write a c++ code for this although not very familiar with it. Any other idea?    Thanks in advance,

Hi R-users,   I have this code here: library(numDeriv)   fprime <- function(z) { alp  <- 2.0165;   rho  <- 0.868;   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))^a   bes1   <- besselI(z*c2,a)   t1bes1 <- t1*bes1   c1*t1bes1*t2 }   ## Newton

Hi R-users,   I have this code below and I understand the error message but do not know how to correct it.  My question is how do I get rid of “with absolute error < 7.5e-06” attach to value of cdf so that I can carry out the calculation.   integrand <- function(z) { alp  <- 2.0165   rho  <- 0.868   # simplified expressions   a      <- alp-0.5   c1     <-

HI Guys, I know that this forum is not for homework but I am trying to interpret R output code. I was just wondering if someone might be able to help. I have been given the following. For (X1,X2) distributed bivariate normal with parameters mu1 = 5.8 mu2 = 5.3 sd1 = sd2 = 0.2 and p = 0.6 The r-code and inpit/output are as follows input m <- 5.3 + 0.6*(6.3 - 5.8) s <-

Haj i try to perform a principal component analysis by using a tetrachoric correlation matrix as data input tetra <- tetrachoric (image_na, correct=TRUE) t_matrix <- tetra$rho pca.tetra <- principal(t_matrix, nfactors = 10, n.obs = nrow(image_na), rotate="varimax", scores=TRUE) the problem i have is to compute the individual factor scores from the pca. the code runs perfect

Hi, Does anyone know what the distribution for the product of two correlated normal? Say I have X~N(a, \sigma1^2) and Y~N(b, \sigma2^2), and the \rou(X,Y) is not equal to 0, I want to know the pdf or cdf of XY. Thanks a lot in advance. yu [[alternative HTML version deleted]]

Hi, I draw a surface corresponding to bivariate density of independent variables (rho=0) using persp(). Then I add a contour line (i.e., circle in my case) at a particular density. Below is a minimal example of what I have so far. # Bivariate density dnorm2d <- function(x, y, rho = 0) { xoy = (x^2 - 2 * rho * x * y + y^2)/(2 * (1 - rho^2)) density = exp(-xoy)/(2 * pi * sqrt(1 - rho^2))

About two years ago there was a thread about this which suggested that at that time nobody had these coefficients ready to go. (a) has anyone in the meanwhile programmed them? (b) I think I can see how to do the tetrachoric one with mvtnorm on similar lines to an example on the help page so will try that if nobody else already has (c) looking at the polychoric one makes me realise yet again

Hi, I would like to know if there is a program written in R to get the CDF (cumulative distribution function) of a bivariate non-central chi-square distribution. Hope someone will reply. Thank you, Rossita M Yunus yunus@usq.edu.au This email (including any attached files) is confidentia...{{dropped:19}}

Hi, I have a dataset with four categories of data, the number of samples are not the same in each category. I want to find the Spearaman's Rho. Let me give an example. x=(14.22770439,26.49420624,46.7277932,19.02550707,23.37379361,16.97789862,19.77100085,23.11270162,13.72929843,33.54430621,14.4756979,70.15811106,11.22789833,NA,NA,NA)