[R] multivariate integration and partial differentiation

[Ravi Varadhan]

Hi,

The "adapt" package might work, but note that it cannot handle infinite limits.  So, integrate() is your best bet.  Since you need to integrate out only one dimension, integrate() would work just fine. 

As for differentiation, you might try the grad() fucntion in the "numDeriv" package.

Ravi.

____________________________________________________________________

Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins University

Ph. (410) 502-2619
email: rvaradhan at jhmi.edu


----- Original Message -----
From: Wei-han Liu <weihanliu2002 at yahoo.com>
Date: Saturday, March 7, 2009 9:38 am
Subject: [R] multivariate integration and partial differentiation
To: r-help at r-project.org


> Hi R Users:
>  
>  Could somebody share some tips on implementing multivariate 
> integration and partial differentiation in R? 
>  
>  For example, for a trivariate joint distribution (cumulative density 
> function) of F(x,y,z), how to differentiate with respect to x and get 
> the bivariate distribution (probability density function) of f(y,z). 
> Or integrate f(x,y,z) with respect to x to get bivariate distribution 
> of (y,z).
>  
>  Your sharing is appreciated.
>  
>  Wei-han Liu
>  
>  
>        
>  	[[alternative HTML version deleted]]
>  
>  ______________________________________________
>  R-help at r-project.org mailing list
>  
>  PLEASE do read the posting guide 
>  and provide commented, minimal, self-contained, reproducible code.