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

Gabor Grothendieck ggrothendieck at gmail.com writes: Microsoft Word's graphics editor can edit R graphics saved in metafile format, wmf. That includes x axis labels, etc. In theory, yes, in practice: not really. Vertical texts are horizontal after editing, and clipped points may

[R] graphical user interface

Hi there, I would like to design a graphical user interface to provide the functionality of my R-scripts to other users (not knowing anything about R...) Can anybody give me some hints which possibilities I have (advantages / disadvantages)? What do you use to build applications? Ciao, Antje

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

On Tue, 5 Aug 2008, Dieter Menne wrote: Gabor Grothendieck ggrothendieck at gmail.com writes: Microsoft Word's graphics editor can edit R graphics saved in metafile format, wmf. That includes x axis labels, etc. In theory, yes, in practice: not really. Vertical texts are horizontal after

Re: [R] Extracting variable names of final model in stepAIC

The 'final model' returned by stepAIC is just a model fit, so you do this the same way as any fit. I'll assume you want to know for lm() fits, but this is fairly general. library(MASS) example(stepAIC) formula(quine.stp) attr(terms(quine.stp), term.labels) (You are making a habit of asking

Re: [R] graphical user interface

There is quite a few discussions on this in the list archives, particularly on the R-SIG-GUI mailing list. It would help to know something about your application... * how complex is it? * how are the results to be presented? * is it intended to be used be only a few people in one location, or

[R] About Creating a List by Parsing Text

Hi all, I have the following data in which I want to parse and store them in a list __DATA__ print(comp.ll) [1] \tGene 11340 211952_at RANBP5 k= 1 LL= -970.692 [2] \tGene 11340 211952_at RANBP5 k= 2 LL= -965.35 [3] \tGene 11340 211952_at RANBP5 k= 3 LL= -963.669 [4]

Re: [R] graphical user interface

Hi Felix, thanks a lot for the hint with the mailing list. I was not aware of its existence... I'll give it a closer look. To answer your questions: - Complexity: I would prefer to be able to design on each level of complexity. - Results: TextFiles and Graphs, ... also open - Usage: Only for a

Re: [R] graphical user interface

On Tue, 5 Aug 2008, Antje wrote: Hi there, I would like to design a graphical user interface to provide the functionality of my R-scripts to other users (not knowing anything about R...) Can anybody give me some hints which possibilities I have (advantages / disadvantages)? Can you first

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

Prof Brian Ripley ripley at stats.ox.ac.uk writes: On Tue, 5 Aug 2008, Dieter Menne wrote: Gabor Grothendieck ggrothendieck at gmail.com writes: Microsoft Word's graphics editor can edit R graphics saved in metafile format, wmf. That includes x axis labels, etc. In theory,

[R] boxplot with average instead of median

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. Is there a way to do this? Chad Junkermeier, Graduate Student Dept. of Physics West Virginia

Re: [R] simulate data based on partial correlation matrix

Hi Benjamin, Creating 0 correlations is easier and always possible, but creating arbitrary correlations can be done as well (when possible - see below). Suppose that x1,x2,x3,x4 have mean 0 and suppose that the desired correlations are r = (r1,r2,r3,r4). Let A be an orthogonal 4x4 matrix such

Re: [R] Trouble when I call a function within another

Yep, sorry I forgot it :wistle: so the data are data frames with dim(8,1000) called DATA.ValX and DATA.ValY (part of the files attached as a xls file) http://www.nabble.com/file/p18826528/ValX.xls ValX.xls http://www.nabble.com/file/p18826528/ValY.xls ValY.xls When I call the function Test

Re: [R] boxplot with average instead of median

Use bxp() and feed it with data that is not from boxplot.stats(). Therefore you might want to invent some alternative function along the lines of boxplot.stats(). Uwe Ligges Chad Junkermeier wrote: I really like the ease of use with the boxplot command in R. I would rather have a boxplot

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

Does this get you close to what you want: x - read.table(textConnection(Gene 11340 211952_at RANBP5 k= 1 LL= -970.692 + Gene 11340 211952_at RANBP5 k= 2 LL= -965.35 + Gene 11340 211952_at RANBP5 k= 3 LL= -963.669 + Gene 12682 213301_x_at TRIM24 k= 1 LL= -948.527 + Gene 12682

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

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 - \tGene 11340 211952_at RANBP5 k= 1 LL= -970.692 y - lapply(split(p, p[3]), [[, 8) y list() - Gundala Viswanath Jakarta - Indonesia On Tue, Aug

[R] Confidence interval for the coefficient of variation

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. The R-code: library(MBESS) n = 10 ci.cv(mean = 100, sd =

[R] qqline function doesn't plot

I have a data vector x. When I try qqline(x) I get the following error: Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet And a blank plot appears. Can anybody help? What am I doing wrong? Thanks, Scotty

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

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 - \tGene 11340

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 me the

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 other

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: Excellent! I

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

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 of the plot, as it is implemented in ggplot2. Is there a way to

Re: [R] qqline function doesn't plot

Hi qqline does not do the plot. It adds a line to an existing plot. So you have to do something else first, like: qqnorm(x) qqline(x) JS --- -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Scotty Nelson Sent: 05 August 2008 09:58 To:

Re: [R] Long Range Dependence: Hurst exponent estimation

tiu == tolga i uzuner [EMAIL PROTECTED] on Mon, 4 Aug 2008 20:04:15 +0100 writes: tiu Dear R Users, Can anyone point me to a package for R tiu vrsion 2.7.1 which implements some Hurst exponent tiu estimation methods ? o fracdiff -- has been the first package to do so, o

Re: [R] backslash in character string?

Please read the R for Windows FAQs. Uwe Ligges zack holden wrote: Dear list, After searching many old posts, I can't find the solution to a simple problem. can someone tell me how to create a character string with multiple backslashes, as in: file_dir - c(C:\files\data\) I need to

[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

[R] unexpected problem

Dear R users, I have run into a very unexpected problem and I was hoping someone could explain it to me. I have a 650 000 by 12 matrix and I want to perform a rolling regression on it, width 36 or 48, using the package performanceAnalytics. ie:

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 AM,

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

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

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

Katrien Baert katrien.baert at 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.

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 ripley at stats.ox.ac.uk writes: On Tue, 5 Aug 2008, Dieter Menne wrote: Gabor Grothendieck ggrothendieck at gmail.com writes: Microsoft Word's graphics editor can edit R graphics saved in

[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

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

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

[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

[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

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 will

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

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 -954 -954 -949

[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.

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 the

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] 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 inspect for

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 - -),

[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

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.org

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] Subject:

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 k l m n o p q r s t u v w

[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] 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

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] 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

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

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 of the

[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

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 (using the max

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 file some

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 function to

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,

[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

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.

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: A)

[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:

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] 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

[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

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

[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

[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

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. Can

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 is on both

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

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

[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:

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

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

[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 http://usa.autodesk.com/adsk/servlet/ps/item?siteID=123112id=9583842linkID=9242018 id=9583842linkID=9242018 Da: Adam Terando [mailto:[EMAIL PROTECTED] Inviato:

[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

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

[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] 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

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), type

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,

[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

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

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] 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

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

[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 have

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 available

[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

[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 **

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

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 location

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 the

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

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 an

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:

  1   2   >