similar to: R and Wordpress

Dear R experts, I have plot the PDF of t distribution with df = 74. curve(dt(x,df=74),from=-4, to=4) how can I shade the area under curve (for example, col="red") from t=+- 1.996? Thank you. Regards, CH -- CH Chan Research Assistant - KWH http://www.macgrass.com

I am getting an unexpected gap in a simple plot of monthly data series. I have created a csv file of monthly climate observations that I store on-line. When I download the csv file and plot one of the series, I get a gap even though there is data for the missing point. Here is a snippet to show the problem. ## Strange plot results link <-

I am trying to build an easy to use climate data analysis tool kit that will let non-R users run my detailed r script with minimum R learning curve effort. Here's an example: link <- "http://chartsgraphs.wordpress.com/files/2010/01/nsidc_trend_plot_2.doc" source(link) This lets user run my R script stored at site 1 and work with data stored at site 2 without having to worry

Dear R experts, Suppose I have a data.frame recording the date and test results of some subjects like this: Name Date results John 01/01/1991 2 John 02/01/1991 3 John 09/0101991 4 Micheal 02/01/1991 4 Micheal 04/01/1991 5 .... How to select the earliest (or latest) test result from all subjects? Thank you. Regards, CH -- "The scientists of today think deeply instead of clearly. One

Hi all, A friend just sent me this: http://www.mathworks.com/company/events/webinars/index.html?id=&language=en <http://www.mathworks.com/company/events/webinars/index.html?id=&language=en>And asked me if there is something of the like in the R community. Does anyone know of such a think ? Cheers, Tal ----------------Contact

Dear R Gurus, I have a problem related to plot. For example, I have two variables, pre and post. pre <- c(1,2,3,4,5) post <- c(2,5,7,2,3) How can I plot a line graph similar to this one? http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1847566&rendertype=figure&id=F1 Would you please provide me a keyword so that I can search. I don't know what is the name of this

Hi everybody, somehow i dont get the shapiro wilk test for normality. i just can?t find what the H0 is . i tried : shapiro.test(rnorm(5000)) Shapiro-Wilk normality test data: rnorm(5000) W = 0.9997, p-value = 0.6205 If normality is the H0, the test says it?s probably not normal, doesn ?t it ? 5000 is the biggest n allowed by the test... are there any other test ? ( i know qqnorm

Dear experts, I would like to know how to plot the log-minus-log plot for survival analysis (to check the proportional assumption) in R. Using the AML example. fit <- survfit(Surv(time, status) ~ x, data=aml) length(fit$surv) #20 as the length of fit$surv is shorter than aml$x and aml$time. I don't know how to plot. Thank you. Regards, CH -- CH Chan Research Assistant - KWH

Hello dear R help members (and also Yihui and Romain), There are currently 28 R bloggers (out of the 117 R-bloggers<http://www.r-bloggers.com/>I know of) that are using wordpress.com for publishing their R code (and I suspect this number will increase with time). WordPress.com doesn't support R syntax highlighting, nor can it be embedded from other services (like gist

I'm sure this is a really simple problem, but I've spent hours digging and I keep running into roadblocks. I'm trying to get a simple chart with three time series. Similar to the attached example http://www.nabble.com/file/p25398419/Excel%2Bchart%2Bexample.pdf Excel+chart+example.pdf , something that was quite easy to do in Excel, except that I need a log y-axis: something that R

Hi all, LLVM now is using Sphinx to generate it's documentation. For syntax hightlighting, Sphinx uses Pygments [2]. Pygments already support LLVM IR syntax, however it doesn't support LLVM tblgen. I wrote up a Pygments extension that can be incorpreated in Sphinx [3]. After installing the extension, you can write the code below and get hightlight tblgen syntax. .. code-block:: tblgen

Dear R Users, I would like to know is there any package to create a plot like this? http://dl.dropbox.com/u/5409929/cs1160521f01.gif X axis is categorical. And the positions of the points are corresponding to the frequency. (similar to violinplot) Thank you. Regards, CH -- CH Chan

Hi all, In R, is there some functions or ways to create a Clustered-Stacked Column Chart as the example in the following page http://peltiertech.com/Excel/ChartsHowTo/ClusterStack.html? I have browsed the R Graph Gallery (http://addictedtor.free.fr/graphiques/) and searched the R site, and didnot find an appropriate method to do it. Anybody has met this problem before? Thanks a lot.

Hello, I am trying to do a logistic regression and have one predictor variable (x) that is ratio and two predictor variables (y and z) that are categorical. These have three levels each which I have called "High", "Medium" and "Low". My question: do I need to use a numerical coding scheme for the categorical variables as required by some statistical software

Dear R experts, Suppose I have an data frame likes this: > example <- data.frame(age=c(1,2,3, 4,5,6), height=c(100,110,120,130,140,150), disease=c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE)) > example age height disease 1 1 100 TRUE 2 2 110 TRUE 3 3 120 TRUE 4 4 130 FALSE 5 5 140 FALSE 6 6 150 FALSE Is there anyway to compare the age and

I am trying to download and open an on-line netcdf file. I'm using Windows XP and R 2.11.1 Here's my script library(ncdf) link <- "http://ibis.grdl.noaa.gov/SAT/SeaLevelRise/slr/slr_sla_gbl_free_all_66.nc" dest <- "C:/temp/slr_sla_gbl_free_all_66.nc" download.file(url=link,destfile=dest) nc1 <- open.ncdf(dest) The file appears in my C:/temp

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1,

Dear R users, I have a very simple question and I've tried to search for the answer. (But failed.) there should be a function (func) that work like > abc <- c(1,2,3,4) > func(abc) "abc" I would like to know the name of that function. Thank you very much for your help. Regards, CH -- CH Chan

Hi R Users,   Thanks in advance.   I am using R-2.11.1 on Windows XP.   May I request you to assist me for the following please.   1.       Is there any R-package or if any to compose music?   2.       Is there any R-package or if any to analyse music?   Once again, thank you very much for the time you have given.   Regards,   Deb   [[alternative HTML version deleted]]

What's the best way around "firstboot" ? "firstboot" hung my machine and had to blown it up to get past it. If you bypass firstboot, what other config settings is one missing? I found the following info.... ==================================================================== # Booting without Internet Connection / Stuck on RedHat Registration during First Bootup This