similar to: how to send email with R

Hi I would like to plot the variation of some mean values with time, and have the standard deviation around the mean shaded on the plot. I could not find a way to have the shaded area on the curve with the default R commands, do I need a special package to do that? Or any idea of a way with the default R commands? Many thanks Thomas -- Thomas Loridan King's College email: thomas.loridan

Is it possible to produce the same graphic in different formats in one shot? Now if I want to produce the same graphic in pdf, ps and jpeg I run the producing code, but I change the graphic device, and I choose successively : pdf(...) ps(...) jpeg(...) So in this case I run the same code 3 times, or I have to choose a format and use converters... Is it possible to produce the 3 graphics in once

I all, I want to create a graph layout in a 3x3 matrix like this: ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| ___ ___ ___ ylab |__| |__| |__| xl xl xl With this layout, then I'll insert the 9 plots. How ca I create it? -- CdeB

Hello, I calculate splines from messured points(x,y) of an unknown function f(x). e.g. x <- c(0.004115, 0.012345, 0.037037, 0.111110, 0.333330, 1.000000) y <- c(37, 50, 45, 60, 50, 66) w <- c(0.8540541, 0.8320000, 0.8822222, 0.7983333, 0.8220000, 0.8151515) as weights f <- smooth.spline(x,y,w) Now I have y-values and want to calculate the x-value(s) from this spline. There is no

Hi I have been using RGoogleDocs successfully for some time now but something seems to have happened which is preventing me from accessing my data in google spreadsheets. I get the message: "problems connecting to get the list of documents" when I use getDocs, despite being logged in e.g. sheets.con = getGoogleDocsConnection(getGoogleAuth("username", "password",

Does anybory knows any work comparing R with other (charged) statistical softwares (like Minitab, SPSS, SAS)? I work in a brasilian government bureau and I intend to use R as our preferable statistical software, but I have to show it's as good as the others. I also intend to use Weka, and for this one I have the same problem. Can anyone help me? Thanks René M. Raupp e-mail: rener@mpdft.gov.br

Dear all, (Sorry if the question has already been answered.) Could someone please suggest a good text editor for writing R sources ? (I know emacs exists ... but I find it a bit heavy). I use crimson (http://www.crimsoneditor.com) which is small and simple, but the R syntax seems not to be supported. Thanks for any advice

Dear R-help readers, I'm sure this problem has been answered but I can't find the solution. I have two vectors v1 <- c("a","b") v2 <- c(1,2,3) I want an easy way to produce every possible combination of v1, v2 elements Ie I want to produce c("a1","a2","a3", "b1","b2","b3") regards Desmond Desmond

Dear list, I would like to start some R workshops at King's College London, and to do so, I would like to use the "Use R!" logo at http://www.agrocampus-ouest.fr/math/useR-2009//useR%21%202008_fichiers/useR-middle.png Since it seems to be difficult to get a shell account at KCL, I also went ahead and registered use-r.org.uk and am starting to put together a website at

Dear all, One can use ks.test(x,y) for a one-sample kolmogorov-smirnov test: x being the data sample y being a string specifying a distribution I notice the help on ks.test does not tell you how to get such a list. Is this a hole in my R knowledge? Where can I get a list of the strings specifying the possible distributions? and more specifically What would be the string and following

I appreciate that this is OT, but I'd be grateful for pointers to examples of where Sweave has been used for web-based applications. In particular, examples of where reports/analyses are produced automatically through submission of data to a web-sever. I am mostly interested in situations where pdf reports have been produced rather than, say, a plot/table etc shown on a web page.

Hi all, I am currently trying to read, write and append data between R and MS access using the RODBC library functions. I have no problems reading in the data but when using sqlSave and sqlAppend it doesn't seem to work. I have made sure that all the column names are sensible and there are no gaps etc in the variables. My call looks like this: sqlSave(channel,treatlist,test=F) I've

About two years ago there was a thread about this which suggested that at that time nobody had these coefficients ready to go. (a) has anyone in the meanwhile programmed them? (b) I think I can see how to do the tetrachoric one with mvtnorm on similar lines to an example on the help page so will try that if nobody else already has (c) looking at the polychoric one makes me realise yet again

I have a list of groups of xy positions I want to set to 0 in an array full of 1s. When the assignments are done directly from the list, they are incorrect, while if I use a temporary array derived from the list the assignments are correct. The following example will hopefully make my problem clearer. The matrices z and zz are initialised with 1. The z and zz values at coordinates in

Hello, Is it possible to draw level lines on the surface of a bivariate function? In the following example, to draw surface and levels lines for a multivariate normal law, I use persp, trans3d, contourLines and lines, but if the lines are correctly drawn, some parts of them are, of course, visible even if they are drawn on a non visible "face". Any suggestion to avoid this problem

Timothy, I believe you are mistaken. Fisher's exact test give the correct answer even in the face of small expected values for the cell counts. Pearson's Chi-square approximates Fisher's exact test and can give the wrong answer when expected cell counts are low. Chi-square was developed because it is computationally "simple". Fisher's exact test, particularly with tables

Folks: I consider my reply below rather clumsy: One has to keep track of index numbers other than that which is inserted and must separately change column names. Is there as "essentially better" way to do this, either via base R or via an R package. I leave it to you to define "essentially better." Thanks. Cheers, Bert On Mon, Aug 1, 2011 at 10:17 AM, Bert Gunter

I generate a groupedData object library(nlme) obj <- groupedData(mg10 ~ time | gp, data = common, outer = ~pct) gp has 101 levels, and pct has 3. There are 38, 25, 38 gps in each of the levels of pct respectively. I fit my model fit.rtg <- lme(mg10 ~ time * group, data = obj, random = ~time * group | gp) Now I try to plot the results. I would like to print 40 panels on each

Début du message réexpédié : > De: Depiereux Constant <constant.depiereux@aqte.be> > Date: 16 mars 2005 00:11:11 GMT+01:00 > À: spjgmwn@iop.kcl.ac.uk > Objet: Rép : [R] RODBC, sqlSave and sqlAppend > > > Le 15-mars-05, à 16:19, Matthew W Nash a écrit : > >> Hi all, >> >> I am currently trying to read, write and append data between R and MS

Full_Name: Cuvelier Etienne Version: R version 2.2.1 OS: Windows XP Submission from: (NULL) (81.240.71.204) If we use the "density" function "far" from the data, some values are negatives. I know that these value of the density are "near equal" to zero, but the change of sign can involve some "false" decision. Here is an example of this fact: > x =