search for: mvnormtest

Dear all, I want to save the p-value from mshapiro.test(mvnormtest). But mshapiro.test(mvnormtest) gives a list with class "htest" containing statistic, p.value, method and data.name as a whole. How do I save only p-value from the result? Thanks a lot!

Dear all, I have a multivariate dataset containing 100,000 or more points. I want find the p-value for the dataset of points coming from a particular multivariate normal distribution With mean vector u Covariance matrix s2 So H0: points ~ MVN( u, s2) H1: points not ~ MVN( u, s2) How do I find the p-value in R? To me this is a likelihood ratio test problem. In H0 the parameters are

Hi, there. I have two questions about the diagnostics in multivarite statistics. 1. Is there any diagnostics tool to check if a multivariate sample is from multivariate normal distribution? If there is one, is there any function doing it in R? 2. Is there any function of testing if two multivariate distribution are same, i.e. the multivariate extension of Kolomogrov-Smirnov test? Thanks for

Hello, Could someone help me to explain the VERY big difference in applying two tests on multivariate normality: library(mvnormtest) data(EuStockMarkets) mshapiro.test(t(EuStockMarkets[15:29,1:4])) Shapiro-Wilk normality test data: Z W = 0.8161, p-value = 0.005955 and library(energy) mvnorm.etest( EuStockMarkets[15:29,1:4] ) Energy test of multivariate normality: estimated parameters data: x, sample si...

is there any function under R that allows me to test the normality of my 92 sumples? -- View this message in context: http://www.nabble.com/multinormality-tp17717230p17717230.html Sent from the R help mailing list archive at Nabble.com.

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

...imulated data from bivariate normal 1000 times and used mshapiro test to see how many times the test would reject the null hypothesis when the p-value is 0.05. The answer should be around 50 since the p-value is 0.05. But I got a much higher value. Here is the R code I used and the result. library(mvnormtest,lib.loc="~/mshapiro") library(MASS) n=50 dim=2 ntrial = 1000 x<-matrix(1:(dim*n),ncol=dim) count = 0 for( trial in 1:ntrial) { x<-mvrnorm(n,rep(0,dim),diag(1,dim,dim)) data_hn<-x p<-mshapiro.test(t(...

My stat textbook tells me that using Shapiro-Wilk test for each variable one by one is not equal to a test for multivariate normality as a whole. Does R have a function of testing for multivariate normality? Thanks. Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571

Hola! Estoy teniendo un problema para testar la normalidad y la homogeneidad de varianzas para mi modelo de manova Es este mi modelo: m1<-manova(cbind(pred,fit,oniv)~tratamento/parcela+bloco,data.vegetacao) Aparte de tener tres variables respuestas (predador, fitófago, omnívoro), también es un diseño anidado. Alguien sabría responder esto? Gracias! Carol [[alternative HTML version

Dear all, I got this error message > library(dprep) > mardia(Savg) Error in cov(data) : 'x' is empty But with the same data, I got > library(mvnormtest) > mshapiro.test(Savg) Shapiro-Wilk normality test data: Z W = 0.9411, p-value = 0.6739 What does the error message "Error in cov(data) : 'x' is empty" mean? Thanks a lot! Jiao

Dear All,   I often have to work with certain models in which I try to "reproduce" a distribution the best I can with very little known information avaible. Is there a package or function in R that could best reproduce a probability distribution using only the mean, median and SD values availble without knowing the actual distribution type to begin with and/or the covariance matrix (for

...odeltools","moments","monoProc","monreg","mprobit","MPV","MSBVAR","msm","muhaz","multcomp","multilevel","multinomRob","multtest","mvbutils","mvnmle","mvnormtest","mvoutlier","mvpart","mvtnorm","NADA","ncomplete","ncvar","negenes","network","neural","nFDR","NISTnls","nlme","nlmeODE","nlreg","nnet","n...

...0-1), ggmap (2.0), growcurves (0.1.1), gWidgetsRGtk2 (0.0-81), hash (2.1.0.1), HiddenMarkov (1.7-0), ISIPTA (0.1-7), lavaan (0.4-13), lcmm (1.5.1), lifecontingencies (0.9.2), LIHNPSD (0.2.1), LMERConvenienceFunctions (1.6.8.2), markdown (0.4), medAdherence (1.03), MKmisc (0.91), modTempEff (1.5.1), mvnormtest (0.1-9), ncdf4 (1.6), ncdf4 (1.6.1), negenes (1.0-1), oposSOM (0.1-4), pcaPP (1.9-46), phenology (3.29), polySegratioMM (0.6-2), pomp (0.41-3), psgp (0.3-1), QCA (1.0-1), QCA3 (0.0-5), quantreg (4.79), Rcmdr (1.8-4), RcppArmadillo (0.3.0.1), RcppArmadillo (0.3.0), RcppDE (0.1.1), RDieHarder (0.1.2)...

CRAN (and crantastic) updates this week New packages ------------ * allan (1.0) Alan Lee http://crantastic.org/packages/allan Automates Large Linear Analysis Model Fitting * andrews (1.0) Jaroslav Myslivec http://crantastic.org/packages/andrews Andrews curves for visualization of multidimensional data * anesrake (0.3) Josh Pasek http://crantastic.org/packages/anesrake This