x <- iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] y <- iris$Petal.Length[1:50]/iris$Petal.Width[1:50] I want to check whether (x,y) follows a bivariate normal distribution or not, using density plot or scatter plot. Is it possible to plot a bivariate density in R. I cant find any. Arindam Fadikar M.Stat Indian Statistical Institute. New Delhi, India [[alternative HTML version deleted]]
On Apr 6, 2010, at 7:48 AM, arindam fadikar wrote:> x <- iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] > y <- iris$Petal.Length[1:50]/iris$Petal.Width[1:50] > > I want to check whether (x,y) follows a bivariate normal > distribution or > not, using density plot or scatter plot. Is it possible to plot a > bivariate > density in R. I cant find any.wireframe and package functions for mv-normal distribution?> >require(MASS) ?kde2d OR ??kde which on my machine pulls up quite a 2 alternatives in my collection of installed packages: ade4::s.kde2d, KernSmooth::bkde2D, and maybe ks::kde.> > Arindam Fadikar > M.Stat > Indian Statistical Institute. > New Delhi, India-- David Winsemius, MD West Hartford, CT
On Tue, Apr 6, 2010 at 12:48 PM, arindam fadikar <arindam.fadikar at gmail.com> wrote:> x <- iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] > y <- iris$Petal.Length[1:50]/iris$Petal.Width[1:50] > > I want to check whether (x,y) follows a bivariate normal distribution or > not, ?using density plot or scatter plot. Is it possible to plot a bivariate > density in R. I cant find any. >Perhaps there is something in mvnormtest? Liviu