similar to: R summary (and quantiles)

Dear all i use R only a few days and don't understand the difference between fivenum(x) und summary(x). > x [1] 20.77 22.56 22.71 22.99 26.39 27.08 27.32 27.33 27.57 27.81 28.69 29.36 [13] 30.25 31.89 32.88 33.23 33.28 33.40 33.52 33.83 33.95 34.82 > fivenum(x) [1] 20.770 27.080 29.025 33.280 34.820 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 20.77 27.14

Last Friday, Gregory Chaitin (http://www.umcs.maine.edu/~chaitin/lm.html) mentioned that there can be no proof that a given code is the shortest for a problem, even within a language. Still, the script below, a replacement of the "TDT", one of the most frequently used tests in genetics (http://mustat.rockefeller.edu under "downloads") may get close. It contains a few

Why does the assignment of a 3178x93 object to another 3178x93 object remove the dimension attribute? > GT <- array(dim = c(6,nrow(InData),ncol(InSNPs))) > dim(GT) [1] 6 3178 93 > SNP1 <- InSNPs[InData[,"C1"],] > dim(SNP1) [1] 3178 93 > SNP2 <- InSNPs[InData[,"C2"],] > dim(SNP2) [1] 3178 93 > dim(pmin(SNP1,SNP2)) [1] 3178 93

The function 'wilcox.test' in R and S gives (almost) identical results (see below). 'qwilcox' however, does not: > qwilcox(p,5,5) p: 0.025 0.975 -------------------- R> 3 22 S> 18 37 I originally wanted to ask a questions, but then I found the answer. Given the confusion I run into, I wonder if this experience is worth reporting. The

Dear R users, I've written the following functions to implement Silverman's modality test ("Using kernel density estimates to investigate multimodality", J.R. Stat. Soc. B, 43, (1981), 97-99.), and I tested them on the Chondrite data set (Good & Gaskins, J. Amer. Stat. Ass., 75, (1980), 42-56). Values for the critical window width seem OK, which is not the case for the

Dear List, i am trying to produce a 3d plot using wireframe using the code: wireframe(Residuals_FD ~ Elevation * Temperature, data = data2, scales = list(arrows = FALSE), drape = TRUE, colorkey = TRUE) As you can see when the code (using the data below) is run the plot area is set-up correctly but the actual surface is missing? Any help would be greatly appreciated. Chris #data Elevation

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Hi, I'm trying to receive faxes using hylafax / iaxmodem but I just can't get it to work. We're using Sangoma E1 cards and have calls coming in over PSTN. I've tried turning hardware echo cancellation off but it makes no difference. This is what I get in /var/spool/hylafax/log: [root at faxhost log]# cat c000000003 Jan 14 12:44:43.39: [ 3403]: SESSION BEGIN 000000003 18005551212

I have a new Sun Fire AMD64 server running Solaris 10 12/06 and Samba 3.0.21b. I have migrated a configuration file from an older RHEL machine running Samba 2.2.x. On this new machine, if one tries to copy a large number of files from a Windows XP client to the server, the copy starts and runs for a few minutes, but then stops. The error message on the XP machine is "Cannot copy

Dear useRs, how can I construct a box&whisker plot based on the vector "fivenum"? The challenge I face is as follow: I have a table such as x | fivenum --------------- ... | (.....) ... | (.....) and so forth.... For each observation x I have generated a vector containing the fivenum estimates. The first challenge is to group my fivenum vectors into groups based on a

Hi everyone: I'll try to be both succinct and accurate here: Attempting to get a directory listing from a DFS export fails using "Samba 3.0", but works fine on windows (XP, 2K, NT4). This happens via smbclient and via "smbmount"ed directories. I had no problem joining this domain, and the user that I am connecting as has permissions to open these folders (I checked

Full_Name: Richard Huggins Version: 1.7.1 OS: windows 2000 Submission from: (NULL) (131.172.4.44) > x<-rnorm(100,2,1) > mean(x) [1] 1.73299 > summary(fivenum(x)) Min. 1st Qu. Median Mean 3rd Qu. Max. -0.3655 1.1070 1.7430 1.7320 2.3840 3.7910 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. -0.3655 1.1070 1.7430 1.7330 2.3830 3.7910 >

You can get anything you want, at Alice's Restaurant: > quantile(d,.75,type=1) 75% 913 > quantile(d,.75,type=2) 75% 981.5 > quantile(d,.75,type=3) 75% 913 > quantile(d,.75,type=4) 75% 913 > quantile(d,.75,type=5) 75% 981.5 > quantile(d,.75,type=6) 75% 1015.75 > quantile(d,.75,type=7) 75% 947.25 > quantile(d,.75,type=8) 75% 992.9167 >

Hi everyone: I'd appreciate any help here... Using smbclient I can change into a directory within a DFS share, but when I attempt to list it's contents, I receive the subject error. I can view the directory (folder) from a win client using the same credentials. I already posted much of my configuration, sans the "smb.conf", although if anyone suspects that there are

Hi. I'm going to introduce the R-package for a group of medical doctors later this week and is a little confused about there use of a test named "willcoxon-pratt" for testing if the clinical and biochemical markers has decreased significantly after the use of some drugs for a group of patients. Looking into the R-functions I would in R recommand using a matched-pairs Wilcoxon

I've just started using R and am still a neophyte, but I found the following curious result. I'm using the current version of R (2.5.1 (2007-06-27) ). Why are the results for the third quartile different in the output from the summary and fivenum commands? For the following data set 457 514 530 530 538 560 687 745 745 778 786 790 792

Dear R, I'm trying to plot the first principle component of an analysis vs the first variable but am having trouble. I have no trouble doing the initial plot but have difficulty thereafter. First I want to highlight some points of the following data set list(running) [[1]] X100m X200m X400m X800m X1500m X5K X10K Marathon Argentina 10.39 20.81 46.84 1.81

Full_Name: David Leong Version: 2.7.2 OS: Windows XP Submission from: (NULL) (12.187.86.2) It appears for the fivenum function, there is a strait average between two ranked samples. This should be a linear interpolation between ranks. For the following data set #data

Has anyone considered using Markdown.pl or PHPME as an Apache handler for Markdown files? It would be very slick to post straight .text files and have them rendered behind the scenes, although for high traffic sites pre-rendering would obviously be much more efficient. I suppose I'd also want an escape hatch like daringfireball.net's, perhaps disabling rendering if the URL includes

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