Re: [R] How to calculate GLM least square means?

Hi Mark, What I always do is to produce a prediction data set and use > pv <- predict(model, newdata = pdata, se = TRUE) This is pretty straightforward. You just build the data set like > pdata <- with(oldData, expand.grid(group = levels(group), value = mean(value)) You probably ne

Re: [R] How to calculate GLM least square means?

Hi Bill, Thanks for your reply. I do work with SAS people, so I may have picked up some of their jargon. I've also heard least square means referred to as marginal means. I understand them to be group means after having controlled for a covariate (i.e. holding it constant at some typical value of

[R] Replacing double-quote with backtick

Hi, How can I change the back quoted strings below > print(x) [1] "foo" [1] "bar" into [1]`foo` [2]`bar` Because later I want to access a named list with this string: mylist$`foo` mylist$`bar` I can't do it with: mylist$"foo" mylist$"bar" - Gundala Viswanath Jakarta - Indonesia

[R] R CMD check message that I can't figure out (.onLoad failed)

I have been able to track and fix everything so far with R CMD check. How do I make this reproducible and I will do it. this is the message that I got: Error in fun(...) : couldn't connect to display ":0.0" Error : .onLoad failed in 'loadNamespace' for 'tcltk' -- Let's not spend our time and r

Re: [R] boxplot with average instead of median

Have a look at forestplot in the rmeta package. Its not a boxplot but does allow you to plot this sort of info and would probably be less confusing than a non-standard boxplot. On Tue, Aug 5, 2008 at 12:36 AM, Chad Junkermeier <[EMAIL PROTECTED]> wrote: > I really like the ease of use with the bo

[R] About colours in violin and simple violin plots

Dear R users Let’s assume I have the following batches of data a <- rnorm(20,200,100) b <- rnorm(20,250,100) c <- rnorm(20,300,100) # I plot them as violin plots require(vioplot) vioplot(a, b, c) # I plot them as simple violin plots require(UsingR) simple.violinplot(a, b, c) # I plot them as b

Re: [R] How to calculate GLM least square means?

Can you explain what is a "least square mean"? It sounds like SAS talk to me. Also, care to tell us who you really are? Bill Venables http://www.cmis.csiro.au/bill.venables/ "A propaganda é a alma do negocio" -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On B

[R] How to calculate GLM least square means?

Hello R-helpers, I would like to calculate least square means after having built a GLM with quasipoisson errors. In my model the dependent variable is continuous, I have one continuous independent variable and one categorical independent variable (that is the variable for which I would like to ca

Re: [R] boxplot with average instead of median

On Mon, Aug 4, 2008 at 10:36 PM, Chad Junkermeier <[EMAIL PROTECTED]> wrote: > I really like the ease of use with the boxplot command in R. I would rather > have a boxplot that shows the average value and the standard deviation then > the median value and the quartiles. I would suggest that you d

Re: [R] [R-sig-Geo] LIDAR Problem in R (THANKS for HELP)

Let's not cross-post! I'll pick up only one point of general relevance on this list (R-help). On Tue, 5 Aug 2008, Matt Oliver wrote: you can try memory.limit(size=4000) only if you have 4GB of memory on the system This is not guaranteed to solve your problem though Assuming this is Windo

Re: [R] Matching a period in grep...

on 08/05/2008 08:55 PM [EMAIL PROTECTED] wrote: Hi folks, Can anyone enlighten me as to why I get the following when I search for ".csv" at the end of a string? grep("\.csv$","Blah.csv",value=TRUE) [1] "Blah.csv" Warning messages: 1: '\.' is an unrecognized escape in a character string 2:

Re: [R] creating overall title for plots made with par(mfrow=c(2, 2))

Ouch! I had searched the archives and read over ?par and ?plot, but sure missed the post of today. Thanks, Mark On Tue, Aug 5, 2008 at 7:27 PM, Marc Schwartz <[EMAIL PROTECTED]>wrote: > on 08/05/2008 06:05 PM Mark Kimpel wrote: > >> I'm making some plots on the same page and would like to include

[R] Matching a period in grep...

Hi folks, Can anyone enlighten me as to why I get the following when I search for ".csv" at the end of a string? > grep("\.csv$","Blah.csv",value=TRUE) [1] "Blah.csv" Warning messages: 1: '\.' is an unrecognized escape in a character string 2: unrecognized escape removed from "[\.]csv$" R re

Re: [R] qqline function doesn't plot

Thanks! -- View this message in context: http://www.nabble.com/qqline-function-doesn%27t-plot-tp18827175p18842548.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] boxplot with average instead of median

Look at this disscussion from two weeks ago: http://www.nabble.com/adding-the-mean-and-standard-deviation-to-boxplots-td15271398.html Chad Junkermeier wrote: > > I really like the ease of use with the boxplot command in R. I would > rather have a boxplot that shows the average value and t

Re: [R] Opening R from Tinn without setting directory each time

That doesn't make sense. Which version of Tinn-R are you using? >From Options->Main->Application do you get a dialog box with a tab "R", and within that a tab "General"? At the bottom of the "General" tab there is a button under "Rgui". That should bring up a file selection box. Proceed to the

Re: [R] creating overall title for plots made with par(mfrow=c(2, 2))

on 08/05/2008 06:05 PM Mark Kimpel wrote: I'm making some plots on the same page and would like to include an overall title instead of individual main titles as they are similar and their x and y axis labels are sufficient to distinguish them. Is there a way to assign an overall "main" to this p

[R] creating overall title for plots made with par(mfrow=c(2,2))

I'm making some plots on the same page and would like to include an overall title instead of individual main titles as they are similar and their x and y axis labels are sufficient to distinguish them. Is there a way to assign an overall "main" to this page of plots? Mark -- Mark W. Kimpel MD *

[R] predict for garch

Dear R People: Is there a predict method for garch methods, please? I tried the usual "predict(d1.garch,n.ahead=3) but only got values for the original data. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown ma

Re: [R] R: LIDAR Problem in R (THANKS for HELP)

On Tuesday 05 August 2008, Alessandro wrote: > Thank you Dylan > > Do you think Kriging method is not valid to processing a large data-set as > LiDAR data? > > About RST interpolation and NNI do you know a code in R? > > ale Please reply to the list next time. I would suggest going over the avail

[R] write .table

hi there Can i write two objects in one file. one is a list and another is a vector. The vector here is the names of the list.so can i write for(i in 1:60) # i have 60 list elements. { write.table(names(mylist[[i]],mylist[[i]],filename=test.txt,sep="\t",append=TRUE) } my testfile.txt should hav

Re: [R] Time series, least squares line

The zoo package can represent time and plot such time series using both classic and lattice graphics There are three vignettes that come with the package that you can read for more info. Also see ?zoo ?plot.zoo ?xyplot.zoo Regarding ols see ?lm or the dyn or dynlm packages both of which work with

Re: [R] Time series, least squares line

For the first part: Create a two-column data frame of which one contains your data (Y) and the other the time index (time). Then do: plot(Y~time,type="l") ##lower case L if you want a least-squares line, run a linear regression of Y on time and add that line together with lines for the confidence

Re: [R] Time series, least squares line

do you want just a date time stamp or do you want other things. There are many ways to skin this cat, but we are going to need more information... as a first start try looking at the zoo package. stephen On Tue, Aug 5, 2008 at 5:47 PM, Gareth Campbell <[EMAIL PROTECTED]> wrote: > Hello, > > I

Re: [R] RESHAPE cast help.

On Tue, Aug 5, 2008 at 11:59 AM, <[EMAIL PROTECTED]> wrote: > I have a set of data that is basically sales figures for a given year. It has > columns for Yeaqr, Day Of Year, Sku, SubCatetory, and Category. The first few > lines of data look like: > Year DayOfYearSku Quantity CatId

[R] Time series, least squares line

Hello, I have a time-series of standards measured for Refractive index. They are daily standards, however, I didn't run one everyday so some days have no data. I can plot the values, but the x-axis does not represent the correct time series (i.e. it's just an evenly spaced 1,2,3 type axis). I w

Re: [R] LIDAR Problem in R (THANKS for HELP)

On Tuesday 05 August 2008, Alessandro wrote: > Hi All, > > > > I am a PhD student in forestry science and I am working with LiDAR data set > (huge data set). I am a brand-new in R and geostatistic (SORRY, my > background it’s in forestry) but I wish improve my skill for improve > myself. I wish to

[R] qgamma inaccuracy

Hello, I have been working with various probability distributions in R, and it seems the gamma distribution is inaccurate for some inputs. For example, qgamma(1e-100, 5e-101, lower.tail=FALSE) gives: 1.0. However, it seems this is incorrect; I think the correct answer should be 0.0823720296207172

Re: [R] [R-sig-Geo] LIDAR Problem in R (THANKS for HELP)

you can try memory.limit(size=4000) only if you have 4GB of memory on the system This is not guaranteed to solve your problem though With big datasets like lidar, you are much better off getting access to a 64bit system with a ton of RAM (>64GB). Cheers Matt On Tue, Aug 5, 2008 at 1:47 PM,

Re: [R] xyplot key issue - line colors

On Tue, Aug 5, 2008 at 1:53 AM, Karin Lagesen <[EMAIL PROTECTED]> wrote: > > > I have a problem regarding the colors assigned to the lines in the key > to an xy plot. I specify the plot like this: > > xyplot(numbers~sqrt(breaks)|moltype+disttype, groups = type, data = alldata, > layout = c(3,2

Re: [R] Legality Question about R's Open Source GNU GPL License

I just wanted to let everyone know the good news, that I was approved to use R at work by legal and compliance. Now to begin learning here. Speaking of, is there anyone who is an expert at R and would be willing to have our company pay them as a consultant to come to my place of work (in Connect

[R] R/Splus Advanced Programming Course in Seattle: August 14-15

R/Splus Advanced Programming Course in Seattle: August 14-15 Payment due after the class! Email Sue Turner: [EMAIL PROTECTED] Cheers, Elvis Miller, PhD Manager Training. XLSolutions Corporation 206 686 1578 www.xlsolutions-corp.com [EMAIL PROTECTED] __

Re: [R] problem with nested loop for regression

Hi Jim, Thanks very much for your reply and suggestions. Although the statement Preg[,k]<-coef(lm(tt~sel_col)) is applying lm() between the same vectors, the contents of one of the vectors (sel_col) should be changing with each cycle through the outer loop. My mistake is that I didn't specify inf

[R] mixed model hypothesis testing

The following is a summary of what I have gathered about hypothesis testing with mixed-effects models. I would appreciate it if someone can clarify or correct this, or make any further comments on the topic. To test a single fixed effect: 1) Likelihood-ratio test (anova) using ML (not REML) is app

[R] xyplot key issue - line colors

I have a problem regarding the colors assigned to the lines in the key to an xy plot. I specify the plot like this: xyplot(numbers~sqrt(breaks)|moltype+disttype, groups = type, data = alldata, layout = c(3,2), type = "l" , lwd = 2, col = c("gray", "skyblue"), key = simpleKey(levels(

[R] I: R: [R-sig-Geo] R: LIDAR Problem in R (THANKS for HELP)

I found this link to resolve: '3 GB Switch Windows Vista' http://usa.autodesk.com/adsk/servlet/ps/item?siteID=123112 &id=9583842&linkID=9242018 Da: Adam Terando [mailto:[EMAIL PROTECTED] Inviato: ma

Re: [R] Mixed model with multiple response variables?

Gang If you want to analyze all variables simultaneously and account for some correlational structure among the different response variables, then the best strategy is to pre-whiten the data and then use lmer. The methods for pre-whitening are described in detail in Pinhiero and Bates in the GL

Re: [R] Greek characters in plots

You need R >= 2.7 on Windows, but it works for me on both XP and Vista. As the posting guide says, we assume that you have an up-to-date R. On Tue, 5 Aug 2008, michael watson (IAH-C) wrote: Hi Brian Thanks again for the advice, it is appreciated. I had come across U0394 code through the char

[R] R: [R-sig-Geo] LIDAR Problem in R (THANKS for HELP)

My notebook is: Hp Pavilion dv6700 notebook PC Intel() Core ™2 Duo CPU T9300 °2.50gHz 2,50GHz RAM: 4.00 GB OS: 32bit Windows (TERRIBLE!) VISTA Da: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Per conto di Matt Oliver Inviato: martedì 5 agosto 2008 12.02 A: Alessandro Cc: r-hel

Re: [R] Labeling lines in xyplot

Deepayan, Thanks for your help, works great! Cheers -- John Deepayan Sarkar wrote: On Tue, Aug 5, 2008 at 11:15 AM, John Poulsen <[EMAIL PROTECTED]> wrote: Hello, I apologize in advance as this example seems really elementary. Below I have created a simple scatterplot with lines. I would l

Re: [R] Greek characters in plots

Excellent, many thanks, that works a treat! But I don't understand how... >bquote(.(tag) ~ Delta * .(suffix)) "A)" ~ Delta * "stbA::cat" Yet when I use: plot(1, main = bquote(.(tag) ~ Delta * .(suffix))) It works fine in the plot. Must be magic ;) -Original Message- From: Gabor Groth

Re: [R] Greek characters in plots

Hi Brian Thanks again for the advice, it is appreciated. I had come across U0394 code through the character map, but I get: > "\u0394" Error: invalid \u sequence So when I execute your code, I get: > a <- "A)" > b <- "stbA::cat" > paste(a, "\u0394", b) Error: invalid \u sequence This i

Re: [R] Labeling lines in xyplot

On Tue, Aug 5, 2008 at 11:15 AM, John Poulsen <[EMAIL PROTECTED]> wrote: > Hello, > > I apologize in advance as this example seems really elementary. Below I > have created a simple scatterplot with lines. I would like to label each > line with the name of the village, instead of using a legend.

Re: [R] Greek characters in plots

On Tue, 5 Aug 2008, michael watson (IAH-C) wrote: OK, I am definitely struggling here. The text I am trying to plot isn't a mathematical expression, it's: A) ΔstbA::cat OK, so you want a Greek character and not a Greek symbol. As I said they are different (different fonts, different shapes

[R] Mixed model with multiple response variables?

Hi, I have a data set collected from 10 measurements (response variables) on two groups (healthy and patient) of subjects performing 4 different tasks. In other words there are two fixed factors (group and task), and 10 response variables. I could analyze the data with aov() or lme() in package nl

[R] Labeling lines in xyplot

Hello, I apologize in advance as this example seems really elementary. Below I have created a simple scatterplot with lines. I would like to label each line with the name of the village, instead of using a legend. Can someone please let me know how to do it. Thanks in advance -- John Po

Re: [R] Selecting Variables

I think that you have to be a little more explicit with a description of your data. I am not clear as to what this means: > There are lots of variables between each exposure and the values are nominal > with upto 6 values.. Can you provide a more complete description. How many columns of exposu

[R] RESHAPE cast help.

I have a set of data that is basically sales figures for a given year. It has columns for Yeaqr, Day Of Year, Sku, SubCatetory, and Category. The first few lines of data look like: Year DayOfYearSku Quantity CatId Category SubCategory 1 2007 1 1000911 10862

[R] LIDAR Problem in R (THANKS for HELP)

Hi All, I am a PhD student in forestry science and I am working with LiDAR data set (huge data set). I am a brand-new in R and geostatistic (SORRY, my background it’s in forestry) but I wish improve my skill for improve myself. I wish to develop a methodology to processing a large data-set o

Re: [R] Smartest way to evaluate question forms

This is still too wide-open a question. :) However there is any number of packages that deal with binary data. Again it depends on what you want to do with it. If you just want some scoring of binary data you can do this with tables ( ?tables). From there you can do about anything you want (

[R] Programming Contest

Dear all, for those of you who have too much time (nobody) and those who enjoy participating in programming contests (maybe some of you), I found an interesting site linked from reddit's programming section: http://www.reddit.com/r/programming/ The actual contest is at the site: http://selfex

Re: [R] Greek characters in plots

Try bquote / substitute as has been mentioned: tag <- "A)" suffix <- "stbA::cat" plot(1, main = bquote(.(tag) ~ Delta * .(suffix 2008/8/5 michael watson (IAH-C) <[EMAIL PROTECTED]>: > OK, I am definitely struggling here. > > The text I am trying to plot isn't a mathematical expression, it's:

Re: [R] Greek characters in plots

Or, maybe I could just draw it with lines() -Original Message- From: [EMAIL PROTECTED] on behalf of michael watson (IAH-C) Sent: Tue 05/08/2008 1:37 PM To: Prof Brian Ripley Cc: r-help@r-project.org Subject: Re: [R] Greek characters in plots OK, I am definitely struggling here. The

[R] PDF append help

On Tue, 5 Aug 2008, Rajasekaramya wrote: Is there any function to append the pdf file. I want to write in a pdf file some 100 plots(in one single pdf containing 100 plots) while all the plot are created using a for loop. I can create 100 pdf one for each for each plot using a for loop but i wa

Re: [R] PDF append help

Hi, just don't turn if off immediately, see example below: Cheers, Gary # example pdf(file="/home/user/tmp/moreplots.pdf") for (i in 1:5) plot(rnorm(100)) dev.off() -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rajasekaramya Sent: Tuesday, August 05,

Re: [R] PDF append help

Just put all the commands within the pdf()- dev.off() ? --- On Tue, 8/5/08, Rajasekaramya <[EMAIL PROTECTED]> wrote: > From: Rajasekaramya <[EMAIL PROTECTED]> > Subject: [R] PDF append help > To: r-help@r-project.org > Received: Tuesday, August 5, 2008, 12:41 PM > hi there, > > Is there any fu

Re: [R] PDF append help

Hi Ramya, Perhaps pdf("C:/100plots.pdf") for(i in 1:100) plot(rnorm(10), type='b', main='My 100 plots') dev.off() HTH, Jorge On Tue, Aug 5, 2008 at 12:41 PM, Rajasekaramya <[EMAIL PROTECTED]>wrote: > > hi there, > > Is there any function to append the pdf file. > > I want to write in a pdf f

[R] PDF append help

hi there, Is there any function to append the pdf file. I want to write in a pdf file some 100 plots(in one single pdf containing 100 plots) while all the plot are created using a for loop. I can create 100 pdf one for each for each plot using a for loop but i want only one pdf with 100 plots.

Re: [R] request for fine panel axis controls in lattice

On 8/5/08, Chuck Cleland <[EMAIL PROTECTED]> wrote: > On 8/5/2008 8:37 AM, Chosid, David (FWE) wrote: > > > I'm trying to use fine axis controls in lattice for each panel. > > Specifically, within each panel, I want to set the limits for x and y > > equal to each other since it is paired data (usin

[R] Fix for nls bug???

Hi All, I've hit a problem using nls. I think it may be a restriction in the applicability of nls and I may have found a fix, but I've been wrong before. This example is simplified to the essentials. My real application is much more complicated. Take a function of matrix 'x' with additional ar

Re: [R] xyplot strip=function for two conditioning variables

On 8/5/08, Henning Wildhagen <[EMAIL PROTECTED]> wrote: > > Dear Deepayan, > > thanks for the hint with the wrapper in latticeExtra, it is very convenient > and the plot looks fine. Just from intuition i think that it would look even > nicer if the strips at the left side appear at the right site

Re: [R] max size of a matrix

On Tue, 5 Aug 2008, rostam shahname wrote: Hi R users, I am trying to create a matrix, but R has problem with the size of dim, wondering if there is anything that I can do? No. See ?"Memory-limts", and consider a sparse matrix (e.g. package Matrix). Had this worked 'a' would have used 8*14

Re: [R] optimize simultaneously two binomials inequalities using nlm( ) or optim( )

Just a quick follow up to Spencer's post, you might want to look at the AcceptanceSampling package on CRAN: http://cran.r-project.org/web/packages/AcceptanceSampling/index.html HTH, Marc Schwartz on 08/05/2008 09:00 AM Spencer Graves wrote: I saw your post on 7/29, and I have not seen a

Re: [R] Exporting data to a text file

On Tue, 5 Aug 2008, John Kane wrote: Don't worry about the warning. R is just telling you that you're appending something :) I'm not sure why it does that. Perhaps just to make sure that you know that you reallly want to append ? Correct: such a file is close to unreadable (e.g. not by read.

Re: [R] max size of a matrix

Are you aware that that matrix will have 147456^2 elements each of size 8 bytes ("double") resulting in R trying to allocate (147456^2)*8/1024^3 = 162 GB of RAM? If you are aware of this and still trying to allocate a large matrix, it is unfortunately too large due to "technical" limitations in R.

Re: [R] Selecting Variables

Thanks for the help guys, i think i needed to be a bit more explicit however (sorry) There are lots of variables between each exposure and the values are nominal with upto 6 values.. And to add to the problem the datasets i deal with range from anything upto 5G. My guess is that the melt functio

[R] max size of a matrix

Hi R users, I am trying to create a matrix, but R has problem with the size of dim, wondering if there is anything that I can do? > a <- diag(147456) Error in array(0, c(n, p)) : 'dim' specifies too large an array Thanks for your help Rostam [[alternative HTML version deleted]]

Re: [R] Iterating Named List

'lapply' is your friend: lapply(y, function(.element) {...process each element of the list}) > y <- list(a=1:10, b=letters, c=LETTERS) > invisible(lapply(y, function(.ele) {print(.ele); cat("--\n")})) [1] 1 2 3 4 5 6 7 8 9 10 -- [1] "a" "b" "c" "d" "e" "f" "g" "h" "i" "j" "

Re: [R] Exporting data to a text file

Don't worry about the warning. R is just telling you that you're appending something :) I'm not sure why it does that. Perhaps just to make sure that you know that you reallly want to append ? --- On Mon, 8/4/08, pacomet <[EMAIL PROTECTED]> wrote: > From: pacomet <[EMAIL PROTECTED]> > Subje

Re: [R] R: log and Log Histogram

It looks like x is not a numeric. Try str(testground) to see what the data set looks like. --- On Mon, 8/4/08, Alessandro <[EMAIL PROTECTED]> wrote: > From: Alessandro <[EMAIL PROTECTED]> > Subject: [R] R: log and Log Histogram > To: "'David Scott'" <[EMAIL PROTECTED]> > Cc: r-help@r-project

[R] Iterating Named List

Hi all, I have the following named list: > print(y) $`200052_s_at` [1] -1066.975 -1063.893 -1062.815 -1062.121 -1059.004 $`200071_at` [1] -959.823 -953.980 -953.886 -948.781 -974.890 $`200084_at` [1] -1135.804 -1132.863 -1128.197 -1128.633 -1125.890 What I want to do is to iterate this name l

Re: [R] Selecting Variables

I am not sure where the "Max" comes from, but this might be a start for you: > x <- read.table(textConnection("ID Exposure_1 Exposure_2 Exposure_3 + 1y yy + 2y y- + 3y - -"), hea

Re: [R] Selecting Variables

Try this: data.frame(ID = x$ID, Max_Exposure = max(rowSums(x[,-1] == "y")), Unique_Exposure = rowSums(x[,-1] == "y")) On Tue, Aug 5, 2008 at 10:21 AM, Michael Pearmain <[EMAIL PROTECTED]> wrote: > Hi All, > > i have a dataset that i want to dynamically inspec

Re: [R] optimize simultaneously two binomials inequalities using nlm( ) or optim( )

I saw your post on 7/29, and I have not seen a reply, so I will attempt a response to the question at the start of your email: obtain the smallest value of 'n' (sample size) satisfying both inequalities: (1-alpha) <= pbinom(c, n, p1) && pbinom(c, n, p2) <= beta, where alpha, p1,

Re: [R] re quest for fine panel axis controls in lattice

Hi David, >> Specifically, within each panel, I want to set the limits for x and y >> equal to each other since it is paired data (using the max value of the >> two). In addition to the code Chuck Cleland sent you, you may want to "square" things up by adding the argument: aspect = "iso" before

[R] Add arrows to a filled.contour margin

Hi, I try to put an arrow to a filled.contour plot. I make this: filled.contour(volcano, color = terrain.colors,key.axes=F) arrows(0.95,0,0.95,1,lwd=2) This work, but the problem is that I find the coordinates at random try, for other graphic with other scale is other fight with coordinates. An

Re: [R] Extracting Names of Named List (Index of Hash)

names(y) On Tue, Aug 5, 2008 at 10:14 AM, Gundala Viswanath <[EMAIL PROTECTED]> wrote: > Hi, > > How can I extract the names of the following > list: > >> str(y) > List of 815 > $ 200052_s_at : num [1:5] -1067 -1064 -1063 -1062 -1059 > $ 200071_at: num [1:5] -960 -95

[R] Opening R from Tinn without setting directory each time

Hi - I can access R from Tinn-R by going to Options->Main->Application/R and setting the search path, but each time I exit Tinn-R I have to redefine the search path. Is there no way of fixing that directory as default? I have installed R under its default directory C:/Program Files/R/R-2.7.1 and

Re: [R] boxplot with average instead of median

Another option is to modify panel.bpplot in the Hmisc package and specify library(lattice) bwplot(..., panel=mypanel) Note that panel.bpplot will show the mean. It shows more quantiles than a standard box plot so you get more than a 3-number summary. If you show the mean and standard deviati

Re: [R] R init file and source()

On Tue, 5 Aug 2008, Prof Brian Ripley wrote: > This relates to the order of things being done. See ?Startup. > > Note that the default packages are not loaded until after the > Rprofile file has been run, and postscript() is in package > grDevices, not in graphics. Using grDevices::postscript

[R] Selecting Variables

Hi All, i have a dataset that i want to dynamically inspect for the number of variables that start with "Exposure_" and then for these count the entries across each case i.e ID Exposure_1 Exposure_2 Exposure_3 1y yy 2y y

[R] Extracting Names of Named List (Index of Hash)

Hi, How can I extract the names of the following list: > str(y) List of 815 $ 200052_s_at : num [1:5] -1067 -1064 -1063 -1062 -1059 $ 200071_at: num [1:5] -960 -954 -954 -949 -975 $ 200084_at: num [1:5] -1136 -1133 -1128 -1129 -1126 ... Thus in the

Re: [R] Random number generation

Please don't cross-post. I've replied in R-devel. On 8/5/2008 8:47 AM, subramanian R wrote: Hi All, I have a few queries regarding Random Number generation in R. according to the help(Random.User) i defined my own functions for user_unif_rand and user_norm_rand (uniform an

Re: [R] request for fine panel axis controls in lattice

On 8/5/2008 8:37 AM, Chosid, David (FWE) wrote: I'm trying to use fine axis controls in lattice for each panel. Specifically, within each panel, I want to set the limits for x and y equal to each other since it is paired data (using the max value of the two). Of course, I have no problems settin

Re: [R] Greek characters in plots

OK, I am definitely struggling here. The text I am trying to plot isn't a mathematical expression, it's: A) ΔstbA::cat This is biological notation for a particular mutation. The "A)" part is a tag identifying the graph, and the "stbA::cat" bit is the gene that is mutated. Both the tag and g

[R] Random number generation

Hi All, I have a few queries regarding Random Number generation in R. according to the help(Random.User) i defined my own functions for user_unif_rand and user_norm_rand (uniform and normal distribution) But what i figured out was even when i call rexp,rpois,rgeom

Re: [R] Are there any guis out there, which will allow editing of the graph?

If interested in the very very good looking Graphics , try silverlight graphics using this open source package- http://www.visifire.com/silverlight_chart_designer.php On Tue, Aug 5, 2008 at 6:00 PM, Gabor Grothendieck <[EMAIL PROTECTED]>wrote: > On Tue, Aug 5, 2008 at 3:16 AM, Dieter Menne > <[

[R] request for fine panel axis controls in lattice

I'm trying to use fine axis controls in lattice for each panel. Specifically, within each panel, I want to set the limits for x and y equal to each other since it is paired data (using the max value of the two). Of course, I have no problems setting the limits for the entire plot but I am having t

Re: [R] Are there any guis out there, which will allow editing of the graph?

On Tue, Aug 5, 2008 at 3:16 AM, Dieter Menne <[EMAIL PROTECTED]> wrote: > Prof Brian Ripley stats.ox.ac.uk> writes: > >> >> On Tue, 5 Aug 2008, Dieter Menne wrote: >> >> > Gabor Grothendieck gmail.com> writes: >> > >> >> >> >> Microsoft Word's graphics editor can edit R graphics saved in >> >> me

Re: [R] Confidence interval for the coefficient of variation

Katrien Baert gmail.com> writes: > > Dear, > > We are trying to determine the (one-sided) CI for the coefficient of > variation in a small sample (say n = 10), with mean 100 and standard > deviation 21. > It appears though that the R-function ci.cv() and our simulation do not > agree. > You

Re: [R] Greek characters in plots

On Tue, 5 Aug 2008, michael watson (IAH-C) wrote: Excellent! I even tried DELTA but not Delta ?plotmath says 'Alpha' - 'Omega' uppercase Greek symbols How could that be made clearer? (Note it carefully says Greek *symbols* not *character*: they are not the same thin

Re: [R] 95% CI bands on a Lowess smoother

If you use loess instead of lowess, you can get standard errors and hence approximate confidence bands by using predict. Example: plot(cars) plx<-predict(loess(dist~speed, data=cars), se=T) lines(cars$speed,plx$fit) lines(cars$speed,plx$fit+2*plx$s, lty=2) #rough & ready CI lines(cars$speed,plx$

Re: [R] About Creating a List by Parsing Text

You need to read in the line with read.table to parse the string. The solution assumes a dataframe . Sent from my iPhone On Aug 5, 2008, at 7:16, "Henrique Dallazuanna" <[EMAIL PROTECTED]> wrote: I think that you need: p <- scan(textConnection(p), what = "") On Tue, Aug 5, 2008 at 7:41

Re: [R] Greek characters in plots

Where delt should be your variable. On Tue, Aug 5, 2008 at 8:23 AM, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote: > If I understand you need this: > > plot(rnorm(10), main = bquote(Delta == .(delt))) > > On Tue, Aug 5, 2008 at 8:20 AM, michael watson (IAH-C) > <[EMAIL PROTECTED]> wrote: >> Excel

Re: [R] Greek characters in plots

Excellent! I even tried DELTA but not Delta Now, how do I go about pasting that with some text that I have in various other variables? -Original Message- From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Sent: 05 August 2008 12:17 To: michael watson (IAH-C) Cc: r-help@r-project

Re: [R] Greek characters in plots

If I understand you need this: plot(rnorm(10), main = bquote(Delta == .(delt))) On Tue, Aug 5, 2008 at 8:20 AM, michael watson (IAH-C) <[EMAIL PROTECTED]> wrote: > Excellent! I even tried DELTA but not Delta > > Now, how do I go about pasting that with some text that I have in various > oth

Re: [R] Greek characters in plots

Try, expression(Delta) On Tue, Aug 5, 2008 at 8:14 AM, michael watson (IAH-C) <[EMAIL PROTECTED]> wrote: > Hi > > I am running an R script that creates 100s of graphs, and I need to use > the greek CAPITAL letter delta in the mtext() function. > > I got as far as expression(delta) but this gives

Re: [R] About Creating a List by Parsing Text

I think that you need: p <- scan(textConnection(p), what = "") On Tue, Aug 5, 2008 at 7:41 AM, Gundala Viswanath <[EMAIL PROTECTED]> wrote: > Thanks Jim, > > But how can I modify this line of yours > > y <- lapply(split(x, x$V3), "[[", 8) > > to suit my 'comp.ll' > > I tried this but fail: >> p

[R] Greek characters in plots

Hi I am running an R script that creates 100s of graphs, and I need to use the greek CAPITAL letter delta in the mtext() function. I got as far as expression(delta) but this gives me the lowercase version. Can anyone help? I'm using R 2.7 on Windows XP Mick Head of Informatics Institute for A

Re: [R] boxplot with average instead of median

boxplot itself is hardwired to produce the boxplot.stats list, and that is not easy to change. To get a different set of stats, you would need to do things in rwo stages: i) create a boxplot object of the type returned by boxplot, but using your own stats ii) call bxp on that object. That's kind

  1   2   >