similar to: About normality tests...

Hi, this is a made-up example. Function "myfun" returns two arguments. Can "apply" be used so that "myfun" is called only once? Thanks Serguei mat<-matrix(runif(50),nrow=10,ncol=5) myfun<-function(x) { mymean<-mean(x) mysd<-sd(x) return(mymean,mysd) } out1<-t(apply(mat,1,function(x) myfun(x)$mymean)) out2<-t(apply(mat,1,function(x)

Dear R-users, please consider the following minimal example: \documentclass[a4paper,titlepage,onecolumn,12pt]{article} \usepackage[italian]{babel} \usepackage{amssymb} \usepackage[utf8x]{inputenc} \usepackage[pdftex]{graphicx} \begin{document} <<label=test, echo=FALSE, results=tex>>= df.data1 <- cbind.data.frame(A = rnorm(18), B =factor(rep(LETTERS[1:6],

I have a dataframe ("mydf") that contains "differences of means". I wish to test whether these differences are significantly different from zero. Below, I calculate the t-statistic for each column. What is a "good" method to calculate/look-up the p-value for each column? mydf=data.frame(a=c(1,-22,3,-4),b=c(5,-6,-7,9)) mymean=mean(mydf) mysd=sd(mydf)

Dear R-users Idea: Analysing tree height frequency with hist(), normal distribution (ks.test & shapiro.test) and skewness (package e1071 - thanks a lot for this useful package)as an indication of possible self-thinning in an experimental tree stand. Problem: Results from the ks.test and the shapiro.test are not comparable (see example of both tests). Tree height is a nice continuous

Hi all T gurus, I would like to test if my dataset is indeed from N(0, 0.011908969). K.S. test gives following result: > ks.test(data, "pnorm", 0, 0.011908969) One-sample Kolmogorov-Smirnov test data: data D = 0.1092, p-value = 1.318e-05 alternative hypothesis: two-sided How ever "Shapiro-Wilk" test give following : >

I have a large amount of data which I break down into a collection of vectors of 100-125 values each. I would like to test the normality of the vectors and compare them. In the interactive mode I can test any one vector using the Shapiro-Wilk test or the Kolmogorov-Smirnov test. My problem is that when I try to write out the results to a file, the term output is a mixture of alpha characters

Hello I have used recently the kolmogorov smirnov test, which is a test of normality. This test is named ks.test() in ctest library of R. I wonder if the results of ks.test () are true, because the results are strange, time to time. thank you for help meriema -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Good Night I made different test to check normality and multinormality in my dataset, but I don´t know which test is better. To verify univariate normality I checked: shapiro.test, cvm.test, ad.test, lillie.test, sf.test or jaque.bera.test and To verify multivariate normal distribution I use mardia, mvShapiro.Test, mvsf, mshapiro.test, mvnorm.e. I have a dataset with almost 1000 data and 9

Hi, does the shapiro wilk test in R-1.3.0 work correctly? Maybe it does, but can anybody tell me why the following sample doesn't give "W = 1" and "p-value = 1": R> x<-1:9/10;x [1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 R> shapiro.test(qnorm(x)) Shapiro-Wilk normality test data: qnorm(x) W = 0.9925, p-value = 0.9986 I can't imagine a sample being

Hi all, I have a dataset of 2000 numbers ( it's noise measured with a scoop ) Now i want to know of my data is normal distributed (Gaussian distribution). I did already: - 68-95-99.7 test - Q-Q-plot and now i used "nortest library" and the Lilli.test() However i don't understad the output? lillie.test(z) Lilliefors (Kolmogorov-Smirnov) normality test data: z D =

Hi, this is my first help request so please bear with me. I've been running some normality tests using the nortest package. For some of my datasets the Cramer-von Mises normality test generates an extremely high probability (e.g., 1.637e+31) and indicates normality when the other tests do not. Is there something I'm misunderstanding or potentially a bug in the code? Below are the

Hello, While doing test of normality under R and SAS, in order to prove the efficiency of R to my company, I notice that Anderson Darling, Cramer Van Mises and Shapiro-Wilk tests results are quite the same under the two environnements, but the Kolmogorov-smirnov p-value really is different. Here is what I do: > ks.test(w,pnorm,mean(w),sd(w)) One-sample Kolmogorov-Smirnov test data: w D

When building R-devel with gcc 5.2.0 (mingw-w64 v4) on Windows, make check fails reg-tests-1b.R at the following check: x <- c(1:2, NA) sx <- sd(x) !is.nan(sx) Here 'sx' should be 'NA' but it is 'NaN'. It turns out this problem only appears when the function is byte compiled with optimization level 3: mysd <- function (x, na.rm = FALSE) sqrt(var(if

I'm using R-2.6.2 on Fedora Linux 9. I've been experimenting with plotmath. I wish it were easier to combine expressions in plotmath with values from the R program itself. There are two parameters in the following example, the mean "mymean" and standard deviation "mystd". I am able to use bquote to write elements into the graph title like mu = mymean and R will

Hi, I couldn't figure out how to use the functions from the ctest library. I'm using the r-base package that comes with debian potato. library("ctest") told me that no such package existed. I checked the CRAN, but no such package was availiable, instead I was told that it would be part of the standard installation. But functions from ctest like shapiro-wilk don't work. The

Hi, I have a problem with Kolmogorov-Smirnov test fit. I try fit distribution to my data. Actualy I create two test: - # First Kolmogorov-Smirnov Tests fit - # Second Kolmogorov-Smirnov Tests fit see below. This two test return difrent result and i don't know which is properly. Which result is properly? The first test return lower D = 0.0234 and lower p-value = 0.00304. The lower 'D'

Hello, I am analysing several samples whose sizes are from 9 to 110. I would like to test their distribution with R, whether they are normal or not. I wonder which test for normality from R should I use . Thank you for help. Samuel Samuel BERTRAND Doctorant Laboratoire de Biomecanique LBM - ENSAM - CNRS UMR 8005 151, bd de l'Hopital 75013 PARIS Tel. +33 (0) 1 44 24 64 53 Fax +33 (0) 1

Hello, I want to do normality test on my data I write this but I don't understand the display of the results ks.test(data,"pnorm") In fact I want to know if my data is a normal distribution. I have to check the p-value or D? Thanks. _____________________________________________________________________________ l [[alternative HTML version deleted]]

Hello. I have tested such form: <form method=''POST''> <input type=''hidden'' name=''myfield'' value=''value1''> <input type=''hidden'' name=''myfield'' value=''value2''> <input type=''submit'' name=''_submit''

Hi All, I have a problem with actionmailer corrupting my html emails. It is adding the character ''3D'' before any quotes (single or double) in my source and smattering ''='' characters throughout as well. Very strange. Here''s an example from one of the emails: ________________ <body> <div id=3D''wrapper''> <div