Re: [R] small sample techniques

This is how I use it prop.test(c(30,23),c(300,300)) Cheers../Murli -Original Message- From: Moshe Olshansky [mailto:[EMAIL PROTECTED] Sent: Thursday, August 09, 2007 12:01 AM To: Rolf Turner; r-help@stat.math.ethz.ch Cc: Nair, Murlidharan T; Moshe Olshansky Subject: Re: [R] small sample

Re: [R] small sample techniques

@stat.math.ethz.ch Subject: Re: [R] small sample techniques -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T Sent: Thursday, August 09, 2007 9:19 AM To: Moshe Olshansky; Rolf Turner; r-help@stat.math.ethz.ch Subject: Re: [R] small sample

Re: [R] small sample techniques

size is large it is still OK to use t-test (because of the Central Limit Theorem) but this is not necessarily true for the small sample size. You could use simulation to find the relevant probabilities. --- Nair, Murlidharan T [EMAIL PROTECTED] wrote: If my sample size is small

[R] small sample techniques

If my sample size is small is there a particular switch option that I need to use with t.test so that it calculates the t ratio correctly? Here is a dummy example? á =0.05 Mean pain reduction for A =27; B =31 and SD are SDA=9 SDB=12 drgA.p-rnorm(5,27,9); drgB.p-rnorm(5,31,12)

Re: [R] t-distribution

: Friday, August 03, 2007 11:01 AM To: Nair, Murlidharan T; Stephen Tucker; r-help@stat.math.ethz.ch Subject: RE: [R] t-distribution Look at the power.examp and run.power.examp functions in the TeachingDemos package. Do these do what you want? If not you can look at the code in them to see how to fill

Re: [R] t-distribution

:11 PM To: Nair, Murlidharan T; r-help@stat.math.ethz.ch Subject: RE: [R] t-distribution yes, or p - seq(0.001,0.999,,1000) x - qt(p,df=9) y - dt(x,df=9) plot(x,y,type=l) f - function(x,y,...) { polygon(x=c(x,rev(x)),y=c(y,rep(0,length(y))),...) } with(data.frame(x,y)[x = 2.3,],f(x,y,col=gray90

Re: [R] t-distribution

, Murlidharan T; r-help@stat.math.ethz.ch Subject: Re: [R] t-distribution p - seq(0.001,0.999,,1000) x - qt(p,df=9) y - dt(x,df=9) plot(x,y,type=l) polygon(x=c(x,rev(x)),y=c(y,rep(0,length(y))),col=gray90) Hope this helps. ST --- Nair, Murlidharan T [EMAIL PROTECTED] wrote: Indeed, this is what I

Re: [R] t-distribution

(y.tmp[b],rep(0,length(y.tmp[b]))),col=gray90) Please let me know if I have made any mistakes. Thanks ../Murli -Original Message- From: Richard M. Heiberger [mailto:[EMAIL PROTECTED] Sent: Thu 8/2/2007 10:25 AM To: Nair, Murlidharan T; Stephen Tucker; r-help@stat.math.ethz.ch Subject

[R] t-distribution

If I have a calculated t can I get the probability associated with it using an R function by giving it the df and t? I know I can do the whole calculation using t.test() or get the t-distribution using qt(). If t=1.11 and df =9 can I get the probability? Thanks../Murli

Re: [R] t-distribution

Indeed, this is what I wanted, I figured it from the function you and Mark pointed me. Thank you both. I am trying to plot it to illustrate the point and I tried this plot(function(x) dt(x, df = 9), -5, 5, ylim = c(0, 0.5), main=t - Density, yaxs=i) Is there an easy way to shade the area under

[R] manipulating arrays

Can I insert an element in an array at a particular position without destroying the already existing element? X-c(1,2,3,4,5) I want to insert an element between 1 and 2. Thanks ../Murli [[alternative HTML version deleted]] __

[R] inter-rater agreement index kappa

Is there a function that calculates the inter-rater agreement index (kappa) in R? Thanks ../Murli [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] debug()

I am having problems using debug() when I am running R under linux. It works fine under the windows environment. Is there anything special that need to specify when running under linux. I would appreciate your help. Cheers../Murli [[alternative HTML version deleted]]

Re: [R] multcomp

Thanks for your help, Cheers../Murli From: talepanda [mailto:[EMAIL PROTECTED] Sent: Sun 3/18/2007 11:07 PM To: Nair, Murlidharan T Cc: r-help@stat.math.ethz.ch Subject: Re: [R] multcomp ?glht says with 'print', 'summary', 'confint', 'coef' and 'vcov' methods

[R] multcomp

I used the multcomp package sometime back for doing multiple comparisons. I see that it has been updated and the methods like simint are no longer supported. When I run the program it prompts to me to use glht. How do I get the lower and upper conf int and the pValues using glht? Does anyone have

[R] Rcmd and memory

I am running a recursive feature selection wrapper on svm for a large data set. The routine gets progressively slower. I am deleting the objects not required. Is there anything else you advise me look for? I tried running it from the command line but it does not print anything on the screen. I

[R] slice usage

I am trying to understand the usage of slice while plotting. Can any one provide me with some explanation to it or point me to a resource where I can read it in greater detail. I have an example from the usage page for SVM. library(MASS) library(e1071) data(iris) m2 - svm(Species~., data =

[R] slice usage

I am trying to understand the usage of slice while plotting. Can any one provide me with some explanation to it or point me to a resource where I can read it in greater detail. I have an example from the usage page for SVM. library(MASS) library(e1071) data(iris) m2 - svm(Species~., data =

[R] Recursive-SVM (R-SVM)

I am trying to implement a simple r-svm example using the iris data (only two of the classes are taken and data is within the code). I am running into some errors. I am not an expert on svm's. If any one has used it, I would appreciate their help. I am appending the code below. Thanks../Murli

[R] reading images in R

Are there functions to read image files in jpg, gif or even a pdf file? Thanks ../Murli [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] labeling graphs

I am trying to add text at specific location on my graph. I know this can be done in R but I can't recollect. I was trying using locator() to identify the position and use identify() but I can get it to work. Can someone jog my memory? Thanks ../Murli [[alternative HTML version

[R] biplot label size

Which is the parameter that is used to decrease the size of ylabs plotted in biplot? I tried playing with cex and cex.lab I am not getting it right pc - princomp(USArrests) biplot(pc, xlabs = rep(, nrow(USArrests)),ylabs=(colnames(USArrests))) Thanks../Murli

Re: [R] princomp/biplot

/2006 2:40 AM To: Rick Bischoff Cc: Nair, Murlidharan T; r-help@stat.math.ethz.ch Subject: Re: [R] princomp/biplot On Sat, 2 Sep 2006, Rick Bischoff wrote: Infinite values? Or NA or NaN values On Sep 2, 2006, at 3:48 PM, Nair, Murlidharan T wrote: I am getting the following error when I

Re: [R] princomp/biplot

I figured it out. I had to specify the first column as my row names. Please ignore my previous mail From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Sun 9/3/2006 2:40 AM To: Rick Bischoff Cc: Nair, Murlidharan T; r-help@stat.math.ethz.ch Subject: Re: [R

[R] labels in biplot

I am unable to give % sign in the label. The percent sign is stripped off in the plot. Is there a way to force this in? Also can I color the labels differently? pc - princomp(USArrests) biplot(pc, xlabs = rep(, nrow(USArrests))) Thanks ../Murli [[alternative HTML version

[R] princomp/biplot

I am getting the following error when I an trying to use princomp princomp(unique.data) Error in cov.wt(z) : 'x' must contain finite values only What do I look for? __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] histograms

I am interested in plotting histograms for the following data Isoform Tumor_65_198 Tumor_50_192 Tumor_80_167 Tumor_80_204 Tumor_95_197 Tumor_70_189 Tumor_90_202 Tumor_40_177 Tumor_60_21 Tumor_70_174 Tumor_70_147 Tumor_50_5 ABCC4-2007 1 1 1 6 1 9 10 1 2 0 10 1

Re: [R] histograms

Thanks ../Murli -Original Message- From: hadley wickham [mailto:[EMAIL PROTECTED] Sent: Friday, September 01, 2006 5:08 PM To: Gabor Grothendieck Cc: Nair, Murlidharan T; r-help@stat.math.ethz.ch Subject: Re: Re: [R] histograms Your data seems to have come through messed up but lets

[R] spectral clustering

Is there a function in R that does spectral clustering? Thanks ../Murli [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] mult comp significance

This has a stats question and a R question. I am sure there are many core statisticians here how would know the answer to this simple question. In determining the significant comparisons using the methods in multcomp, the ones that are designated as significant are the ones that do not intersect

Re: [R] Multcomp

to analyze them and combine significant ones into one plot. -Original Message- From: Greg Snow [mailto:[EMAIL PROTECTED] Sent: Tuesday, July 25, 2006 12:12 PM To: Nair, Murlidharan T Subject: RE: [R] Multcomp Doing: str( fungus.cirec ) Suggests that fungus.cirec$conf.int contains

Re: [R] Multcomp

to analyze them and combine significant ones into one plot. -Original Message- From: Greg Snow [mailto:[EMAIL PROTECTED] Sent: Tuesday, July 25, 2006 12:12 PM To: Nair, Murlidharan T Subject: RE: [R] Multcomp Doing: str( fungus.cirec ) Suggests that fungus.cirec$conf.int contains

Re: [R] Saving R objects

] Sent: Sun 7/23/2006 10:11 PM To: Nair, Murlidharan T Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Saving R objects It depends on what information you want to save and how the program on the other end needs it. For the save version I would at least use ascii = TRUE to get it in a more readable

[R] running jobs in background

Can I run jobs in the background and the check the status of it from time to time in Windows version of R? Thanks ../Murli [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

[R] Saving R objects

I am trying to find the best way to save the follwoing object I am creating library(multcomp) data(recovery) Dcirec-simint(minutes~blanket, data=recovery, conf.level=0.9, alternative=less) I am probably not doing it the most efficient way I think. Here is what I am doing a-print(Dcirec)

[R] multcomp plotting Please help :)

I REALLY NEED HELP WITH THIS PLEASE I am using the multcomp package for doing multiple comparisons (R 2.3.1 windows). Since the data I am handling is huge the number of comparisons are also large. I am interested in: 1 Breaking down my plots to get rid of the clutter that happens when

[R] Multcomp

Here it is again, hope this is more clear I am using the following data (only a small subset is given): Habitat Fungus.yield Birch 20.83829053 Birch 22.9718181 Birch 22.28216829 Birch 24.23136797 Birch 22.32147961 Birch 20.30783598 Oak 27.24047258 Oak 29.7730014 Oak 30.12608508 Oak 25.76088669

Re: [R] Multcomp

-Original Message- From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Sent: Saturday, July 22, 2006 1:37 PM To: Nair, Murlidharan T Cc: R-help@stat.math.ethz.ch Subject: Re: [R] Multcomp On 7/22/06, Nair, Murlidharan T [EMAIL PROTECTED] wrote: Here it is again, hope this is more clear

[R] Multcomp plotting

I am using the multcomp package for doing multiple comparisons. Since the data I am handling is huge the number of comparisons are also large. I am interested in: 1 Breaking down my plots to get rid of the clutter that happens when plotting the entire data set. How do I pass only part of the

[R] multicomp plotting

I am using the multcomp package for doing multiple comparisons. Since the data I am handling is huge the number of comparisons are also large. I am interested in: 1 Breaking down my plots to get rid of the clutter that happens when plotting the entire data set. How do I pass only part of the

[R] multicomp plotting

I am using the multcomp package for doing multiple comparisons. Since the data I am handling is huge the number of comparisons are also large. I am interested in: 1 Breaking down my plots to get rid of the clutter that happens when plotting the entire data set. How do I pass only part of the

[R] multicomp

Is there a multicomp() function in R as in S plus? Thanks ../Murli __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] minimum adequate model using contrasts

I am trying to model a set of data that I am working on using contrasts. Since the data set is too large, I am interested in scripting the entire process to arrive at the minimum adequate model. Has anyone in the group done this before? If you I would appreciate their comments. Cheers ../Murli

[R] minimal adequate model

I am trying to model a set of data that I am working on using contrasts. Since the data set is too large, I am interested in scripting the entire process to arrive at the minimum adequate model. Has anyone in the group done this before? If you I would appreciate their comments. Cheers ../Murli