search for: exmbvs10

Hello, I want to add text annotation about correlation on "pairs" plots. I found that I could pass a function to the "panel" argument of pairs : panel.annot <- function(x, y, ...) { points(x, y, ...) c <- cor.test(x, y) legend("topleft", legend=substitute(rho == r, list(r=sprintf("%.2f", c$estimate))), bty="n") } And then :

Dear Friends, Here is an algorithm for finding concave hulls: http://get.dsi.uminho.pt/local/ Has anyone implemented such an algorithm in R? RSiteSearch('concave hull') didn't reveal one (I think). _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400 Express Parcels

hie all i am trying to carry out a categorical data analysis but my problem is that when in i use the chi squared test some of my expected values are less than 5. is there a test that can handle this situation. the data is not a 2*2 table. its more from the social sciences where you have from strongly agree to strongly disagree. i know i can collapse vthe tables but there is a loss of

...force two regression coefficients to be equal > but opposite in sign? > To: "Woolner, Keith" <kwoolner at indians.com>, "r-help at r-project.org" > <r-help at r-project.org> > Message-ID: > <B37C0A15B8FB3C468B5BC7EBC7DA14CC60F685895B at LP- > EXMBVS10.CO.IHC.COM> > Content-Type: text/plain; charset=us-ascii > > One simple way is to do something like: > > > fit <- lm(y ~ I(x1-x2) + x3, data=mydata) > > The first coeficient (after the intercept) will be the slope for x1, > the slope for x2 will be the negative o...

...imail.org> To: "Farley, Robert" <[147]FarleyR at metro.net>, "[148]r-help at r-project.org" <[149]r-help at r-project.org> Subject: Re: [R] Constrain density to 0 at 0? Message-ID: <[150]B37C0A15B8FB3C468B5BC7EBC7DA14CC633A53DF11 at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" Look at the logspline package. This is a different approach to density estimation from the kernel densities used by 'density', but does allow you to set fixed boundaries. -- Gregory (Greg) L....

...Mon, 29 Sep 2008 10:38:01 -0600 From: Greg Snow <Greg.Snow at imail.org> Subject: Re: [R] Trend graph To: "C.H." <chainsawtiney at gmail.com>, "r-help at r-project.org" <r-help at r-project.org> Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC6132C7170E at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" Try this: pre <- c(1,2,3,4,5) post <- c(2,5,7,2,3) plot( 1:2, range(pre,post), type='n', xlab='', ylab='CRP (mg/dl)', xaxt='n') segments(1,pre,2,post) axis(1, at=1:2, labels=c('Pre'...

I want to take the first row of each unique ID value from a data frame. For instance > ddTable <- data.frame(Id=c(1,1,2,2),name=c("Paul","Joe","Bob","Larry")) I want a dataset that is Id Name 1 Paul 2 Bob > unique(ddTable) Will give me all 4 rows, and > unique(ddTable$Id) Will give me c(1,2), but not accompanied by the name column.

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML

...----------- Message: 45 Date: Wed, 13 Oct 2010 09:32:31 -0600 From: Greg Snow <Greg.Snow at imail.org> To: li li <hannah.hlx at gmail.com>, r-help <r-help at r-project.org> Subject: Re: [R] extract rows of a matrix Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC633E8150FF at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" > newmat <- oldmat[ c(rep(FALSE,19),TRUE), ] Or > newmat <- oldmat[ seq(20, nrow(oldmat), 20), ] -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.snow at imail.org 801.408.81...

...ssage: 55 Date: Wed, 8 Jul 2009 10:09:08 -0600 From: Greg Snow <Greg.Snow@imail.org> Subject: Re: [R] Uncorrelated random vectors     <r-help@r-project.org>, " Luba (AIM SE)Stein" <LUBA.STEIN@allianz.com> Message-ID:     <B37C0A15B8FB3C468B5BC7EBC7DA14CC62128DEB0D@LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="iso-8859-1" The mvrnorm function in the MASS package has an argument to force the generated data to have the exact mean/variance structure as specified which when used with a diagonal variance matrix will generate data that has a 0 (withi...

...Mon, 4 May 2009 11:32:15 -0600 From: Greg Snow <Greg.Snow at imail.org> Subject: Re: [R] normality test for large a large dataset ? To: Martial Sankar <martial100kar at hotmail.com>, rhelp <r-help at r-project.org> Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC61D19A1765 at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="iso-8859-1" See: http://finzi.psych.upenn.edu/R/Rhelp08/archive/136160.html But the issue really comes down to the fact that the questions: "exactly normal?", and "normal enough?" are 2 very different questions (with...

...27 May 2009 13:16:03 -0600 From: Greg Snow <Greg.Snow@imail.org> Subject: Re: [R] alternative to built-in data editor To: Jose Quesada <quesada@gmail.com>, "r-help@r-project.org"     <r-help@r-project.org> Message-ID:     <B37C0A15B8FB3C468B5BC7EBC7DA14CC61D1FC10B2@LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" Have you tried the View function (note the uppercase V). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.snow@imail.org 801.408.8111 > -----Original Message----- > From: r-help-bou...

...Snow@imail.org> To: Michael Costello <michaelavcostello@gmail.com>, "r-help@r-project.org" <r-help@r-project.org> Subject: Re: [R] Looping with incremented object name and increment function Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC6341402BBA@LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" This is FAQ 7.21. The real gem in the answer there is at the end where it tells you that it is easier to just use a list. If your fit1, fit2, fit3, and fit4 were elements in a list then you can just loop through the list eleme...

...1 -0600 From: Greg Snow <Greg.Snow@imail.org> Subject: Re: [R] unicode only works with a second one To: Thomas Steiner <finbref.2006@gmail.com>,     "r-help@stat.math.ethz.ch"    <r-help@stat.math.ethz.ch> Message-ID:     <B37C0A15B8FB3C468B5BC7EBC7DA14CC61CDD76B4D@LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" I don't know how to help with the Unicode issue, but one alternative is the my.symbols function in the TeachingDemos package (see ?ms.male as well as ?my.symbols). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statisti...

...at imail.org> To: Richard and Barbara Males <rbmales at gmail.com> Cc: "r-help at r-project.org" <r-help at r-project.org> Subject: Re: [R] generating correlated random variables from different distributions Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC6335CD080A at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="iso-8859-1" Transforming variables will generally change the correlation, so your method will give you correlated variables, but not exactly at the correlations you specify (though with some trial and error you may be able to get close)....

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

...-0700 From: Greg Snow <Greg.Snow at imail.org> To: Sonja Klein <sonja.klein.07 at aberdeen.ac.uk>, "r-help at r-project.org" <r-help at r-project.org> Subject: Re: [R] How to produce a graph of glms in R? Message-ID: <B37C0A15B8FB3C468B5BC7EBC7DA14CC633FEC591B at LP-EXMBVS10.CO.IHC.COM> Content-Type: text/plain; charset="us-ascii" Look at Predict.Plot (and possibly TkPredict) in the TeachingDemos package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.snow at imail.org 801.408.8111 > -----Original Message-----...