Re: [R] How to install older version of R?

2009-11-19 Thread Paul Hiemstra
Dear Jia-Chiun Pan, First of all, why install an older version of R? In my experience the backward compatibility of R is quite good. So using R 2.10.0 should be ok. But if you really want this you can build R from source. The source can be downloaded from CRAN:

[R] name of a name

2009-11-19 Thread William Simpson
I have quite a complicated problem that's hard to describe. Suppose I have a dataframe d. I want to access the vector d$var, where var is one of the variables in d. Just take for granted that there's a good reason I want to do this as follows. var-DeOxyA xx-paste(d$,var, sep=) mean(xx) [1] NA

[R] Overlaying multiple ROC curves using ROCR

2009-11-19 Thread Matthew Peterson
Hello, I was following an example on The ROCR Package pdf, learning to overlay ROC curves on the same plot using the add = TRUE statement. I used this one: data(ROCR.hiv) attach(ROCR.hiv) pred.svm - prediction(hiv.svm$predictions, hiv.svm$labels) perf.svm - performance(pred.svm, 'tpr',

Re: [R] Unnecesary code?

2009-11-19 Thread Hun S. Tesatte
On Thu, 19 Nov 2009 00:13:27 +0100 Duncan Murdoch murd...@stats.uwo.ca wrote: hunsynte...@hush.com wrote: Dear R-ers, While browsing the R sources, I found the following piece of code in src\main\memory.c: static void reset_pp_stack(void *data) { R_size_t *poldpps = data;

[R] 3D plot, rotatable and with adjustable symbols

2009-11-19 Thread Joel Fürstenberg-Hägg
Hi all, I've tried to make a 3D plot, but have run into some problems. I'd like to have a plot that I can rotate interactively using the mouse, which is possible using the plots3d {R.basic}. However, I would like to change the symbols used as the points, but there's no pch in plot3d().

Re: [R] name of a name

2009-11-19 Thread baptiste auguie
Hi, Try this, d - data.frame(a=1:4, b=3:6) var - a mean(d[var]) ## or, if you are not aware of ## fortune(parse) xx - paste(d$,var, sep=) mean(eval(parse(text=xx))) HTH, baptiste 2009/11/19 William Simpson william.a.simp...@gmail.com: I have quite a complicated problem that's hard to

Re: [R] Presentation of data in Graphical format

2009-11-19 Thread Jim Lemon
On 11/19/2009 03:13 AM, Sunita Patil wrote: Hello Sir I have got 150 observations, got 10 posts/ 6 departments/ tasks vary from 5 to 10, A few of the variables are crossed specially in case of Office boy, where the tasks are like open the door, put on the lights, Yes time variable I have used

Re: [R] name of a name

2009-11-19 Thread William Simpson
Thanks very much Phil and Baptiste! d[[var]] is exactly what I wanted. Sorry for being so inarticulate. I still couldn't describe my problem if I wanted to! Anyway, now I have the solution. Cheers Bill On Thu, Nov 19, 2009 at 9:39 AM, baptiste auguie baptiste.aug...@googlemail.com wrote: Hi,

[R] PLoS, Arial, R linux

2009-11-19 Thread Ted Harding
Greetings All. According to the PLoS (Public Library of Science) Guidelines for Figure preparation, if a figure is submitted as EPS rather than TIFF then Figure text must be in Arial font -- see: [1] http://www.plosntds.org/static/figureGuidelines.action#text and also other sections in that

Re: [R] 3D plot, rotatable and with adjustable symbols

2009-11-19 Thread Remko Duursma
Try the rgl package. r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia

[R] randomForest: impact of bug fixed in version 4.5-13

2009-11-19 Thread J Sousa
The NEWS of the randomForest R library mention that version 4.5-13 fixed a bug in predict.randomForest() when newdata is a matrix with no rownames. I think it corresponds to the difference in files predict.randomForest.R which is the new line if (is.null(rn)) rn - keep As I've been using version

[R] How to analyse/visualize a city budget ?

2009-11-19 Thread Tal Galili
Hello all, I happen to get a (legitimate) hold of a city budget for the (4) years: 2006,2007,2008,2009 The budget holds over 12,000 rows of budget sections with numbers being Zero's positive and negatives. I would like to find something interesting in this dataset. I don't have a clear

Re: [R] Unnecesary code?

2009-11-19 Thread Duncan Murdoch
On 19/11/2009 4:23 AM, Hun S. Tesatte wrote: On Thu, 19 Nov 2009 00:13:27 +0100 Duncan Murdoch murd...@stats.uwo.ca wrote: hunsynte...@hush.com wrote: Dear R-ers, While browsing the R sources, I found the following piece of code in src\main\memory.c: static void reset_pp_stack(void

Re: [R] PLoS, Arial, R linux

2009-11-19 Thread Liviu Andronic
Hello On Thu, Nov 19, 2009 at 11:00 AM, Ted Harding ted.hard...@manchester.ac.uk wrote: While there is a section (Enable the use of Arial in R) in the Guidelines (URL [2]), the instructions assume the presence of Arial .ttf files, not usually the case with Linux. If you have a Windows

Re: [R] PLoS, Arial, R linux

2009-11-19 Thread Deepayan Sarkar
On Thu, Nov 19, 2009 at 4:30 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: Greetings All. According to the PLoS (Public Library of Science) Guidelines for Figure preparation, if a figure is submitted as EPS rather than TIFF then Figure text must be in Arial font -- see: [1]

[R] Problem with sqlSave

2009-11-19 Thread anna_l
Hello, the sqlSave function is used in order to write a dataframe on excel. This function creates worksheets using the attribute tablename and writes the data.frame in it. What I want to do is to create this data.frame but being able in case this worksheet already exists to delete the former

Re: [R] PLoS, Arial, R linux

2009-11-19 Thread Ted Harding
On 19-Nov-09 11:52:57, Deepayan Sarkar wrote: On Thu, Nov 19, 2009 at 4:30 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: Greetings All. According to the PLoS (Public Library of Science) Guidelines for Figure preparation, if a figure is submitted as EPS rather than TIFF then Figure text

[R] [Follow-up] Re: PLoS, Arial, R linux

2009-11-19 Thread Ted Harding
See Addendum at end. On 19-Nov-09 12:17:54, Ted Harding wrote: On 19-Nov-09 11:52:57, Deepayan Sarkar wrote: On Thu, Nov 19, 2009 at 4:30 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: Greetings All. According to the PLoS (Public Library of Science) Guidelines for Figure preparation,

[R] Accessing list names in lapply

2009-11-19 Thread Bjarke Christensen
Hi, When using lapply (or sapply) to loop over a list, can I somehow access the index of the list from inside the function? A trivial example: df1 - split( x=rnorm(n=100, sd=seq(from=1, to=10, each=10)), f=letters[seq(from=1, to=10, each=10)] ) str(df1) #List of 10 # $ a: num [1:10]

Re: [R] Accessing list names in lapply

2009-11-19 Thread Gabor Grothendieck
lapply over the list names rather than the list itself: junk - lapply(names(df1), function(nm) plot(df1[[nm]], ylab = nm)) On Thu, Nov 19, 2009 at 7:27 AM, Bjarke Christensen bjarke.christen...@sydbank.dk wrote: Hi, When using lapply (or sapply) to loop over a list, can I somehow access the

Re: [R] Accessing list names in lapply

2009-11-19 Thread Duncan Murdoch
On 19/11/2009 7:27 AM, Bjarke Christensen wrote: Hi, When using lapply (or sapply) to loop over a list, can I somehow access the index of the list from inside the function? No, but you can loop over the indices in lapply, not just in a for loop. For example, lapply(names(df1),

[R] [Follow-up] Re: PLoS, Arial, R linux

2009-11-19 Thread Dirk Eddelbuettel
On 19 November 2009 at 13:00, (Ted Harding) wrote: | However, | (a) I do not have the 'ttf2afm' program; | (b) Synaptic cannot find it in the repositories. | | So, I'm off on another hunt! i) You should try to get over the 'Synaptic is the only interface to package' syndrome. It limits

Re: [R] [Follow-up] Re: PLoS, Arial, R linux

2009-11-19 Thread Liviu Andronic
On 11/19/09, Ted Harding ted.hard...@manchester.ac.uk wrote: However, (a) I do not have the 'ttf2afm' program; (b) Synaptic cannot find it in the repositories. http://packages.debian.org/search?searchon=contentskeywords=ttf2afmmode=exactfilenamesuite=testingarch=any Liviu

[R] t-criterion calculation using means and SE

2009-11-19 Thread Антон Морковин
Dear all, is there any functions which allow to calculate Student t-criterion using means, their SE and sample size? I've seek for, but bulit-in t-criterion works only with initial sample... Best regards, A.Morkovin __ R-help@r-project.org mailing

Re: [R] Accessing list names in lapply

2009-11-19 Thread Romain Francois
Maybe this : http://tolstoy.newcastle.edu.au/R/e4/help/08/04/8720.html Romain On 11/19/2009 01:27 PM, Bjarke Christensen wrote: Hi, When using lapply (or sapply) to loop over a list, can I somehow access the index of the list from inside the function? A trivial example: df1- split(

Re: [R] Accessing list names in lapply

2009-11-19 Thread Henrique Dallazuanna
You can try this: par(mfcol=c(5,2)) lapply(df1, function(x){ nm - names(eval(as.list(sys.call(-1))[[2]]))[as.numeric(gsub([^0-9], , deparse(substitute(x] plot(x, main = nm) }) On Thu, Nov 19, 2009 at 10:27 AM, Bjarke

Re: [R] re placing the dates format in R for exporting the data set...

2009-11-19 Thread David Winsemius
On Nov 18, 2009, at 11:00 PM, ychu066 wrote: hey Jim , I have solve the column name problems now. But i am still unable to read the date in R ... toms_dat- replace(toms_dat, toms_dat ==2009-08-24, 6) replace needs its first argument to be a vector, while you have given it a

[R] Acceptance sampling

2009-11-19 Thread Bodil Tufvesson
Hi I'm trying to learn how to use R when it comes to acceptance sampling. What I want to do is to type in a lot size, an AQL level and a LQ level and get the required sample size in accordance with inspection level II (ISO 2859-2). I have tried find.plan (c(0.01,0.95), c(0.07,0.05),

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread Hrishi Mittal
Add col=i in the histogram call. -- View this message in context: http://old.nabble.com/How-do-I-change-the-colour-and-format-for-the-trelli-plot---tp26418382p26421302.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] [Follow-up] Re: PLoS, Arial, R linux

2009-11-19 Thread Frank E Harrell Jr
Also see of the journal will let you use NimbusSan: pdf('my.pdf', onefile=FALSE, pointsize=18, family=NimbusSan,height=6,width=8,paper=special) ... dev.off() embedFonts('my.pdf') Frank Dirk Eddelbuettel wrote: On 19 November 2009 at 13:00, (Ted Harding) wrote: | However, | (a) I do not

Re: [R] name of a name

2009-11-19 Thread David Winsemius
On Nov 19, 2009, at 4:39 AM, baptiste auguie wrote: Hi, Try this, d - data.frame(a=1:4, b=3:6) var - a mean(d[var]) I tried get() on Simpson's pasted argument without success. Assuming the OP still wanted to use only text arguments, then this approach works dfn - d get(dfn)[var]

[R] ddply function nesting problems

2009-11-19 Thread Jason Rupert
While putting my R code into functions, I've encountered a ddply function nesting issue and need a bit of advice on the proper way to fix it.  I've tried several approahces, but neither worked and I need to have the ability to include the cut, range, and fullseq methods within ddply.  (For a

[R] t-criterion calculation using means and SE

2009-11-19 Thread Антон Морковин
Dear all, I need to create n*n table with sums of all possible pair combinations of numbers from n-row column. What function allows it? Best regards, A.Morkovin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Presentation of data in Graphical format

2009-11-19 Thread Sunita Patil
Hello Sir Thanx even I will try to work out on your suggestions, will keep you updated on the progress. Thanx a lot Regards Our Thoughts have the Power to Change our Destiny. Sunita On Thu, Nov 19, 2009 at 3:50 PM, Jim Lemon j...@bitwrit.com.au wrote: On 11/19/2009 03:13 AM, Sunita Patil

Re: [R] Presentation of data in Graphical format

2009-11-19 Thread Sunita Patil
Hello Sir Thanx a lot, will try Pareto chart for my data Regards Our Thoughts have the Power to Change our Destiny. Sunita On Thu, Nov 19, 2009 at 12:29 PM, Petr PIKAL petr.pi...@precheza.cz wrote: Well, from what you say it seems to me that you could also use Pareto charts together with

Re: [R] Accessing list names in lapply

2009-11-19 Thread Bjarke Christensen
Thanks to everybody who replied - I got three distinct, very useful suggestions. Bjarke Christensen Romain Francois romain.francois@

Re: [R] t-criterion calculation using means and SE

2009-11-19 Thread David Winsemius
On Nov 19, 2009, at 9:40 AM, Антон Морковин wrote: Dear all, I need to create n*n table with sums of all possible pair combinations of numbers from n-row column. What function allows it? ?expand.grid Best regards, A.Morkovin __

Re: [R] 3D plot, rotatable and with adjustable symbols

2009-11-19 Thread Greg Snow
Look at the ggobi program and the rggobi package for interactions between R and ggobi. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-

Re: [R] ddply function nesting problems

2009-11-19 Thread baptiste auguie
Hi, I think your ddply call with a calculation inside .( ) is the problem. Are you sure you need to do this? Performing the cut outside ddply seems to work fine, determine_counts-function() { min_range-1 max_range-30 bin_range_size-5 Me_df-data.frame(Data =

[R] advice about R for windows speed

2009-11-19 Thread Carlos Hernandez
Dear All, I appreciate any advice or hints you could provide about the following. We are running R code in a server (running Windows XP and QuadCore Xeon processors, see details below) and we would like to use the server efficiently. Our code takes a bit more than 6 seconds per 25 iterations in

[R] Efficient cbind of elements from two lists

2009-11-19 Thread Stephan Dlugosz
Hi! I have a data.frame data and splitted it. data - split(data, data[,1]) This is a quite slow procedure; and I do not want to do it again. So, any unsplit and resplit is no option for me. But: I have to cbind variables to the splitted data from another list, that contains of vectors with

Re: [R] Efficient cbind of elements from two lists

2009-11-19 Thread Jorge Ivan Velez
Dear Stephan, Here is a suggestion using do.call(): res - do.call(cbind, yourlist) res HTH, Jorge On Thu, Nov 19, 2009 at 10:03 AM, Stephan Dlugosz wrote: Hi! I have a data.frame data and splitted it. data - split(data, data[,1]) This is a quite slow procedure; and I do not want to

Re: [R] advice about R for windows speed

2009-11-19 Thread Marc Schwartz
On Nov 19, 2009, at 9:25 AM, Carlos Hernandez wrote: Dear All, I appreciate any advice or hints you could provide about the following. We are running R code in a server (running Windows XP and QuadCore Xeon processors, see details below) and we would like to use the server efficiently.

Re: [R] Efficient cbind of elements from two lists

2009-11-19 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Stephan Dlugosz Sent: Thursday, November 19, 2009 7:03 AM To: r-help@r-project.org Subject: [R] Efficient cbind of elements from two lists Hi! I have a data.frame data and

[R] Problem on using Rexcel in the vba code

2009-11-19 Thread anna_l
Hello everybody, I started to use Rexcel and I am getting an error with the following code: Sub AutoForma1_Clique() Call RInterface.StartRServer Call RInterface.RRun(setwd(C:/Program Files/R/R-2.10.0/Working Directory)) Call RInterface.RRun(getwd()) End Sub The error is the following

[R] Printing labeled summary to text file ?

2009-11-19 Thread A Singh
Dear List, I am trying to run a mixed model which, on the R console, prints output as follows: [1] Marker [1] perm no. [1] NA Linear mixed model fit by REML Formula: peg.no.prm ~ 1 + (1 | family/f) Data: modeldf AIC BIC logLik deviance REMLdev 3119 3134 -1555 31123111 Random

[R] Splitting massive output into multiple text files

2009-11-19 Thread A Singh
Dear List, I thought it would be much easier to put a second query into a second mail. I need to print 426*1 blocks of variance components data, where 426 is the number of columns of data that have 1 permutations of variance generated for each of them. I have tried printing out a

[R] Problem with zoo and BootPR packages

2009-11-19 Thread Ricardo Gonçalves Silva
Hi, I'm trying to plot the forecasts I generated using the Plot.Fore function of the BootPR package. But I got an error from zoo: My data: Time Series: Start = 1 End = 18 Frequency = 1 [1] 38731 38628 39117 92809 71984 31226 58613 72360 107956 92066 [11] 95208 99098 95848

[R] person-level to person-period xfm

2009-11-19 Thread Robert Terwilliger
Dear R experts, I have a so-called person-level data frame that I need to transform into a person-period data frame. If the lingo is unclear, the data have one row for each subject, with repeated measures data each in a separate column. I need to transform these data so that each subject has

Re: [R] Problem with zoo and BootPR packages

2009-11-19 Thread Gabor Grothendieck
Contact the BootPR maintainer regarding a bug in this line of Plot.Fore: y1 - zooreg(x, start, end, frequency) where x is a ts object but that may not be used in that context. as.zooreg is available for converting ts series (and certain other objects) to zooreg objects. 2009/11/19 Ricardo

Re: [R] advice about R for windows speed

2009-11-19 Thread Carlos Hernandez
Thanks for your reply! I just added some more details below. Our code needs around 1GB of RAM and all machines and R configurations have its default maximum above this number. Our suspicion is that the windows server could run the code in half of its current time (given the apparent factor of 2

Re: [R] Problem with zoo and BootPR packages

2009-11-19 Thread Achim Zeileis
On Thu, 19 Nov 2009, Ricardo Gonçalves Silva wrote: Hi, I'm trying to plot the forecasts I generated using the Plot.Fore function of the BootPR package. But I got an error from zoo: My data: Time Series: Start = 1 End = 18 Frequency = 1 [1] 38731 38628 39117 92809 71984 31226 58613

Re: [R] Problem with zoo and BootPR packages

2009-11-19 Thread Ricardo Gonçalves Silva
Ok, Thanks all. Rick. -- From: Achim Zeileis achim.zeil...@wu-wien.ac.at Sent: Thursday, November 19, 2009 3:06 PM To: Ricardo Gonçalves Silva ricard...@terra.com.br Cc: R-Help r-help@r-project.org; j@latrobe.edu.au Subject: Re: [R] Problem

Re: [R] person-level to person-period xfm

2009-11-19 Thread Richard M. Heiberger
Robert Terwilliger wrote: Dear R experts, I have a so-called person-level data frame that I need to transform into a person-period data frame. If the lingo is unclear, the data have one row for each subject, with repeated measures data each in a separate column. I need to transform these data

[R] plot filled.contour over continent map

2009-11-19 Thread Matthias Demuzere
Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp,

Re: [R] Error system is computationally singular by using function dmvnorm

2009-11-19 Thread Uwe Ligges
Alla Bulashevska wrote: Dear R users, i try to use function dmvnorm(x, mean, sigma, log=FALSE) from R package mvtnorm to calculate the probability of x under the multivariate normal distribution with mean equal to mean and covariance matrix sigma. I become the following Error in

Re: [R] converting a vector of bytes to a PNG/JPEG image

2009-11-19 Thread Uwe Ligges
Can't be too hard to convert it in pixmap representation, but you need to tell us hoe the representation looks like, probbaly with a short code that makes it reproducible (e.g. generate some data that we can use in R). Uwe Ligges Rajarshi Guha wrote: Hi, I have some code that uses rJava.

Re: [R] ddply function nesting problems

2009-11-19 Thread Jason Rupert
Awesome!  Thanks a ton!  I guess I had overlooked how it was really working. I will still have to reflect on why it was working running it straight through, but not being nested. That is kind of a mystery.  Oh well... Thanks again.    - Original Message From: baptiste auguie

Re: [R] person-level to person-period xfm

2009-11-19 Thread David Winsemius
On Nov 19, 2009, at 12:23 PM, Richard M. Heiberger wrote: Robert Terwilliger wrote: Dear R experts, I have a so-called person-level data frame that I need to transform into a person-period data frame. If the lingo is unclear, the data have one row for each subject, with repeated measures

Re: [R] Problem with zoo and BootPR packages

2009-11-19 Thread Gabor Grothendieck
Actually it may be that the documentation is at fault more than the code. The help page for Plot.Fore says that the first argument is a time series data set and although that seems to suggest that it should be a ts object the code seems to be written assuming a plain numeric vector; therefore,

[R] How do I specify a partially completed survival analysis model?

2009-11-19 Thread RWilliam
Hello, I just started using R to do epidemiologic simulation research using the Cox proportional hazard model. I have 2 covariates X1 and X2 which I want to model as h(t,X)=h0(t)*exp(b1*X1+b2*X2). I assume independence of X from t. After I simulate Time and Censor data vectors denoting the

[R] Reading in a series of files using a for loop

2009-11-19 Thread Thomas Jackson
Dear R Users, I am trying to read in a series of csv files which vary by the letter on the end of he file name. When I input what seems to be a logical for loop I get an error message that doesn't make sense to me. for(i in 1:12){ paste(GP, LETTERS[i],sep='')

Re: [R] How to install older version of R?

2009-11-19 Thread 潘家群
Dear Paul, Thank you for your response. I am investigating a differentially expressed gene detected method that cannot be applied on R 2.9.1 that currently I used. The reason about why it cannot work due to the changes of new version of R. And, recently, I study a lots of methods that developed

Re: [R] Reading in a series of files using a for loop

2009-11-19 Thread Jorge Ivan Velez
Hi Thomas, Here are two suggestions which storage all the files in a list: # Parameters route - /Users/thomasjackson/Data/GEPCO/GEPCO lts - LETTERS[1:12] toread - paste(route,lts,/HPLC_,lts,12.csv,sep='') # Reading in the data -- option 1 myfiles - list() for(i in 1:12) myfiles[i] -

Re: [R] Reading in a series of files using a for loop

2009-11-19 Thread Greg Snow
Try something like: gp - lapply( paste((/Users/thomasjackson/Data/GEPCO/GEPCO,LETTERS[1:12],/HPLC_,LETTERS[1:12],12.csv,sep=''), read.csv, header=TRUE, sep=',' ) names(gp) - paste(GandP, LETTERS[1:12], sep='') Now gp (or whatever you want to call it) will be a list with your 12

Re: [R] plot filled.contour over continent map

2009-11-19 Thread Ray Brownrigg
On Fri, 20 Nov 2009, Matthias Demuzere wrote: Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary =

Re: [R] Reading in a series of files using a for loop

2009-11-19 Thread David Winsemius
You need to look at the assign function: http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-turn-a-string-into-a-variable_003f (... and in addition to R FAQ 7.21 also perhaps read the rest of the R FAQ.) There are a ton of very similar questions in the r-help archives, so you could

Re: [R] plot filled.contour over continent map

2009-11-19 Thread Greg Snow
You probably want to put the map boundaries on top of the contour plot rather than the other way around. Here is one example of doing that (find the correct asp may be the hardest part): library(maps) x - seq( -124.7, -67, length.out=25 ) y - seq( 25, 49, length.out=25 ) z - outer(x,y,'+')

[R] Lme

2009-11-19 Thread Patrick Bedard
Hi In my experiment I have 2 Groups (control patient) and 4 blocks of trials per subject, n=12 for Control and n=11 for patient The variances of the patients is very much higher than that of the Control. 1- I tried an ANOVA with Groups as between factor and Blocks as within factor (repeated

[R] Calling R (GAMM) from Fortran

2009-11-19 Thread Paul Warren Simonin
Hello,   I am currently working on a modeling project using Fortran to run large repetitive loops (many DO loops). As part of this process I would like to use a model fit in R and currently stored as an R object. This is a rather complex model, a GAMM, and I am curious whether

[R] Partial derivatives of the multivariate cumulative distribution

2009-11-19 Thread Stephane LUCHINI
I'm currently using the mvtnorm package to model unobserved heterogeneity in a structural model and using optim to estimate the model. I have got good clues that convergence is not really a problem but the hessian matrix estimate is very bad. To overcome this problem, I'm constructing an

Re: [R] Calling R (GAMM) from Fortran

2009-11-19 Thread Charlie Sharpsteen
On Thu, Nov 19, 2009 at 12:14 PM, Paul Warren Simonin paul.simo...@uvm.edu wrote:   Hello,     I am currently working on a modeling project using Fortran to run large repetitive loops (many DO loops). As part of this process I would like to use a model fit in R and currently stored as an R

Re: [R] How to wait for a user response in Rscript?

2009-11-19 Thread Peng Yu
Why I don't see 'par(ask=T)' in /usr/lib/R/library/graphics/R-ex/barplot.R? Is 'par(ask=T)' implicitly called by example()? On Mon, Sep 7, 2009 at 10:55 AM, RIOS,ALFREDO ARTURO ar...@ufl.edu wrote: Hi Peng I think this is what you are looking for par(ask=T) Alfredo On Sun Sep 06

[R] Performance of 'by' and 'ddply' on a large data frame

2009-11-19 Thread Tahir Butt
I've only recently started using R. One of the problems I come up against is after having extracted a large dataset (5M rows) out of database, I realize I need another variable. In this case I have data frame with dates. I want to find the minimum date for each value of x1 and add that minimum

[R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Peng Yu
There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but just silently apply a function to the list argument. For example, the plot function is applied to each element in 'alist'. It is redundant to return anything from apply.

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Marc Schwartz
On Nov 20, 2009, at 10:21 AM, Peng Yu wrote: There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but just silently apply a function to the list argument. For example, the plot function is applied to each element in 'alist'.

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Charlie Sharpsteen
On Thu, Nov 19, 2009 at 2:21 PM, Peng Yu pengyu...@gmail.com wrote: There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but just silently apply a function to the list argument. For example, the plot function is applied to

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread jim holtman
invisible(apply(...)) On Thu, Nov 19, 2009 at 5:21 PM, Peng Yu pengyu...@gmail.com wrote: There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but just silently apply a function to the list argument. For example, the plot

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Peng Yu
On Thu, Nov 19, 2009 at 4:27 PM, Marc Schwartz marc_schwa...@me.com wrote: On Nov 20, 2009, at 10:21 AM, Peng Yu wrote: There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but just silently apply a function to the list

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Peng Yu
I want to use both the name and the content. Although, I could do the following thing. for(x in names(List)) { do some thing with x do some thing with List[[x]] } However, I'd prefer something like the following if R offers such functionality. But it seems not. for(x in List) { do

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Gabor Grothendieck
On Thu, Nov 19, 2009 at 5:31 PM, Peng Yu pengyu...@gmail.com wrote: Is there a way to get the name of the list in the loop body? This was just discussed! See this thread: https://stat.ethz.ch/pipermail/r-help/2009-November/218919.html __

[R] Change color of single panel in xyplot

2009-11-19 Thread jimdare
Hi there, I have created a single page, multi-panel, xyplot using lattice. Each panel is a trial (total of 8) from an experiment with 3 variables on the x axis and the observed value on the y axis. I have added the mean of all trials so the entire plot looks like this (where M=mean panel and

Re: [R] Is there an variant of apply() that does not return anything?

2009-11-19 Thread Marc Schwartz
On Nov 19, 2009, at 4:31 PM, Peng Yu wrote: On Thu, Nov 19, 2009 at 4:27 PM, Marc Schwartz marc_schwa...@me.com wrote: On Nov 20, 2009, at 10:21 AM, Peng Yu wrote: There are a few version of apply() (e.g., lapply(), sapply()). I'm wondering if there is one that does not return anything but

Re: [R] How to wait for a user response in Rscript?

2009-11-19 Thread Peter Ehlers
You could try ?example and see what the help page says. -Peter Ehlers Peng Yu wrote: Why I don't see 'par(ask=T)' in /usr/lib/R/library/graphics/R-ex/barplot.R? Is 'par(ask=T)' implicitly called by example()? On Mon, Sep 7, 2009 at 10:55 AM, RIOS,ALFREDO ARTURO ar...@ufl.edu wrote: Hi

[R] loess smoothing

2009-11-19 Thread rkevinburton
Hello, In reading the loess description I see: span: the parameter alpha which controls the degree of smoothing. The default seems to be 0.75. Would it be possible to expand on this decription so I can avoid trail and error? Can I increase this pass 'span' 1? Qualitatively to what degree

[R] optim(.. ,SANN,..)

2009-11-19 Thread lloyd barcza
I have a problem using optim, so I am hoping someone can help me out with it: Suppose I have the list: list(D,R,P) [[1]] V1 V2 V3 V4 1 0 1 0 1 2 1 1 0 0 3 1 0 1 0 4 0 0 1 1 5 0 1 0 1 6 1 1 0 0 7 1 0 1 0 8 0 0 1 1 9 1 0 1 0 10 0 0 1 1 [[2]]

Re: [R] Change color of single panel in xyplot

2009-11-19 Thread David Winsemius
This was just answered by Deepayan earlier today: This code did not produce the plot you have linked to. The answer to your question depends on how you created the plot, so you have to tell us that. Changing the color in all panels is easy: histogram(rnorm(100), col = goldenrod) Different

[R] Passing filenames to the getopt package

2009-11-19 Thread EJ Nikelski
Hi all, I've finally started to use Rscript for my statistical scripting needs, and find I'm being blocked by what must be a very simple problem. Specifically, the command lines for my scripts usually contain: (1) the script name, (2) one or more options and their arguments, and finally, (3)

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread ychu066
what is that ? Hrishi Mittal wrote: Add col=i in the histogram call. -- View this message in context: http://old.nabble.com/How-do-I-change-the-colour-and-format-for-the-trelli-plot---tp26418382p26421453.html Sent from the R help mailing list archive at Nabble.com.

[R] sqldf

2009-11-19 Thread JoK LoQ
Hello, I would like some help with sqldf syntax. Suppose I have table 1 and table 2. What do I have to do to generate a table with columns 2,5,6 from table 1 (for example), and columns 3,4,5,9 from table 2, but only when values in column 2 from table 1 are equal to values in column

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread Hrishi Mittal
To have a different colour for each histogram. Wasn't that your question? Sorry if I misunderstood. Did you try it? ychu066 wrote: what is that ? Hrishi Mittal wrote: Add col=i in the histogram call. -- View this message in context:

Re: [R] Tinn-R related problem

2009-11-19 Thread AlisonT
As of last night I have this exact same problem, Enter, Backspace and the arrow keys do not work for me in TinnR! I am using windows vista prof on a 64bit laptop. I have uninstalled, re-installed many times and even cleaned the registry etc. Does anyone know of a reason and a fix for this

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread ychu066
Fore example my code is histogram(~data[,8]|data[,2], ylab = Frequency, xlab = Score, xlim = c(1,5), ylim = c(0,100),layout = c(3,1),col=data[,2] ) and i want the colour to be depended on the level of the factor in data[,2]. how do i do that ? ychu066 wrote: what is that ?

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread ychu066
I want to change the colour of the histogram in each panel by the levels of the conoditonal factor. I have 3 levels in the factor levels(data[,2])[1] levels(data[,2])[2] levels(data[,2])[3] what can i do , in order to chnage the colour ?? ychu066 wrote: Fore example my code is

[R] RE shaping large dataset

2009-11-19 Thread ychu066
I am doing a project which involve reshaping a large dataset, can any of you please sugguest me some good reading/websites/ examples can be in R and SAS Thanks everyone !!! -- View this message in context: http://old.nabble.com/REshaping-large-dataset-tp26421513p26421513.html Sent from

Re: [R] sqldf

2009-11-19 Thread Gabor Grothendieck
Google for sql join and see the examples in Example 4 on the sqldf home page: http://code.google.com/p/sqldf/#Example_4._Join On Thu, Nov 19, 2009 at 2:30 PM, JoK LoQ jok...@gmail.com wrote: Hello,   I would like some help with sqldf syntax.   Suppose I have table 1 and table 2.   What do

Re: [R] How do I change the colour and format for the trelli plot ?

2009-11-19 Thread David Winsemius
On Nov 19, 2009, at 3:10 PM, ychu066 wrote: I want to change the colour of the histogram in each panel by the levels of the conoditonal factor. I have 3 levels in the factor levels(data[,2])[1] levels(data[,2])[2] levels(data[,2])[3] what can i do , in order to chnage the colour ??

[R] how to plot mean+/-SD from externally calculated values

2009-11-19 Thread Dan Kortschak
Hi, I am trying to plot a set of means+/-SD calculated by an external program (an RDBMS) since the data set is too large to happily fit in R (740M x 100 values - which are summarised to 100 means/SD by the DB). I want to have a mean with whiskers at +/-1SD. Can anyone suggest a way to do this?

Re: [R] how to plot mean+/-SD from externally calculated values

2009-11-19 Thread David Winsemius
On Nov 19, 2009, at 6:33 PM, Dan Kortschak wrote: Hi, I am trying to plot a set of means+/-SD calculated by an external program (an RDBMS) since the data set is too large to happily fit in R (740M x 100 values - which are summarised to 100 means/SD by the DB). I want to have a mean with

[R] Finding replacing non-ASCII characters

2009-11-19 Thread Steven Kang
Hi guys, Are there any feasible methods in searching finding non-ASCII characters in R? For example, from the following object, x - mia. SzaÌmitaÌó The desired output is, x.out - mia. SzaImitaIA Your help in resolving this would be greatly appreciated. [[alternative

  1   2   >