[R] rterm user command line input

In order to get user input when running a batch file with Rterm in Windows XP and Windows Server 2003, I have been using the R script below with R-2.8.1 and earlier. This opens a command line window with the statement Please enter the ID: and feeds the users response into the variable 'id'.

Re: [R] Announcing a new R news site: R-bloggers.com

Hi Dirk, In order to reply to you in a way people will be able to follow, I am going to concentrate my answers (per point) here: 1) PlanetR doesn't encourage the adding of new feeds. No where on the site is there a contact page or a submit new feed form. I give the full credit for the site

[R] savePlot for Mac and / or Linux?

Hi all, In the package rtlu, I use the function savePlot. It is convenient since it let the user decide in which graphic format he wants his graph to be export. But when I run R CMD check, I get the following message : rtlu(V1,fileOutput=First.tex,textBefore=\\section{Variable 1 to

[R] inverse gamma truncated

Hi, all How can generate a sample from truncated inverse gamma distribution in R? thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] yeroon.net/ggplot2 web application v0.11

You seriously rock Jeroen, Way to go! Best, Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com/ (English)

Re: [R] Announce: edtdbg, integrating R's debug() with your text editor

Would it be possible to add this to notepad++ with NPPtoR ? (A tip for whom it may concern, I personally use the debug package and am very happy with it) Cheers, Tal Contact Details:--- Contact me: tal.gal...@gmail.com |

Re: [R] Label individual points in lattice by and ID

If you want to plot each point as the text of its ID, use a panel function: panel = function(x, y, groups, subscripts, ...) panel.text(x, y, groups[subscripts]) If you want to add labels interactively to selected points: xyplot(y~x | f1 + f2, groups=ID,dat) ## repeat for each panel:

Re: [R] savePlot for Mac and / or Linux?

How about pressing ``Print Screen'' to capture image. However, the simple/smart solution for you is to output data and draw the same thing by other softwares. This can save you much time under any operating systems.

Re: [R] yeroon.net/ggplot2 web application v0.11

Jeroen, I just wrote about your update on: http://www.r-statistics.com/2009/12/a-web-application-of-rs-ggplot2/ Also, does your blog at: http://www.stat.ucla.edu/~jeroen/ Has an RSS? (I would love to add it to www.r-bloggers.com) If it doesn't will you (finally) consider to update your blog to

[R] column statistics

Hi everybody, I would like to compute the mean for 1 variable between the rows with the same levels. For example, with the dataset below: Factor1 Factor2 Value A X 1 A X 2 A Y 3 A Y

Re: [R] yeroon.net/ggplot2 web application v0.11

Wow, this is really cool Jeroen! Jeroen Ooms-3 wrote: A new version of the ggplot2 web application is available at http://www.yeroon.net/ggplot2. New features include 1D geom’s (histogram, density, freqpoly), syntax mode (by clicking the tiny arrow at the bottom), and some additional

Re: [R] column statistics

You have several options in R. 1) cast from the reshape package cast(Factor1 + Factor2 ~ . , data = your.data.frame, value = Value, fun = mean) 2) ddply from the plyr package ddply(your.data.frame, c(Factor1, Factor2), function(x){mean(x$Value)}) HTH, Thierry

[R] A GGPLOT ploting question

library(ggplot2) m - ggplot(movies, aes(x=rating)) m + geom_histogram() Now in x-axis instead of numbers like 2,4,6,8,10..., I want to write like 2%,4%,6%,8%,10%... Is there any way to do that through GGPLOT ? Your help will be highly appreciated. Thanks, -- View this message in context:

Re: [R] A GGPLOT ploting question

Hi, I think you have two options: 1- use a specific formatter in scale_x_continuous, e.g. last_plot + scale_x_continuous(formatter=percent) 2- specify breaks and labels, last_plot + scale_x_continuous(breaks=c(2,6), labels=paste(c(2,6),%)) HTH, baptiste 2009/12/7 Megh

[R] install.packages error

Dear listeRs, I am trying to install a package from the command line, using install.packages(). Specifying all the parameters, I simply get this error: r - getOption(repos) install.packages(rJava, dependencies=T, repos=r, lib=/home/adi/myRlib) Error: subscript out of bounds Additional info

Re: [R] Starting estimates for nls Exponential Fit

Thanks everybody, This has been quite helpful, the problem remains tricky but at least now I've got a version of my script that handles all my reactions without error. The DEoptim solution produced good starting values for a lot of reactions, but sadly not for all. I now use scaled parameters

Re: [R] column statistics

Use aggregate: aggregate(x['Value'], x[c('Factor1', 'Factor2')], mean) On Mon, Dec 7, 2009 at 8:25 AM, Ivan Calandra ivan.calan...@uni-hamburg.de wrote: Hi everybody, I would like to compute the mean for 1 variable between the rows with the same levels. For example, with the dataset

Re: [R] Announce: edtdbg, integrating R's debug() with your text editor

On Sun, Dec 06, 2009 at 11:39:06PM -0800, Norm Matloff wrote: I've just developed edtdbg, a small package that integrates R's debug() with one's text editor. Congratulations! I was planning to add this functionality to vim-r-plugin, but I had no idea how to do it. I downloaded edtdbg files from

Re: [R] Sweave, cacheSweave, and data frame

It's been a while since this thread appeared, but I had the same problem and don't want to keep the solution for myself. A small, reproducible example: cache=true= a=1 @ = a+1 @ Compiling with cacheSweave gives an error: library(cacheSweave) Lade nötiges Paket: filehash filehash: Simple

[R] Filtering a zoo object based on index of another object

Hello everybody, I have two datasets, observed and predicted. Since my observed dataset is not in regular intervals, I need to filter my predicted dataset based on the measurement date of my observed data. Here, is an example similar to what I have library(chron);library(zoo) DATE-

Re: [R] Announce: edtdbg, integrating R's debug() with your text editor

On 07/12/2009 2:39 AM, Norm Matloff wrote: I've just developed edtdbg, a small package that integrates R's debug() with one's text editor. Excerpt from the README file: Goal The debug() function in R is primitive. My goal was to make it more usable by integrating it with one's

Re: [R] Filtering a zoo object based on index of another object

On Mon, 7 Dec 2009, Subodh Acharya wrote: Hello everybody, I have two datasets, observed and predicted. Since my observed dataset is not in regular intervals, I need to filter my predicted dataset based on the measurement date of my observed data. Here, is an example similar to what I have

Re: [R] Filtering a zoo object based on index of another object

Since DATE2 is a subset of time(myzoo) either of these should work window(myzoo, DATE2) myzoo[DATE2] If DATE2 were not a subset of time(myzoo) then you could do this: myzoo2 - na.locf(merge(myzoo, zoo(, DATE2))) and now use either of the first two expressions with myzoo2 in place

Re: [R] regression tests still fail for R version 2.10.0

Do you have R_LIBS set to you library of packages? I got this error without fail when trying to R CMD check on my hand compiled R 2.10-patched. R_LIBS was set to my own library and it was full of 2.9.x packages. When I looked in ./tests/reg-tests-1.Rout I got the error about out of date packages

[R] data.frame to stacked frame

Hi, have anybody a hint how i could avoid the cumbersome way (..especially the value part) change the data representation from n column data.frame to an stacked representation. many thanks Christian #example data dfw - as.data.frame(matrix(runif(10*10),ncol=10)) dfw$group -

Re: [R] data.frame to stacked frame

Try this: cbind(stack(dfw, select = -group), dfw$group) On Mon, Dec 7, 2009 at 11:20 AM, Christian Schulz chsch...@email.de wrote: Hi, have anybody a hint how i could avoid the cumbersome way (..especially the value part) change the data representation from  n column data.frame to an

Re: [R] Installing RandomForest on SuSe Linux - warnings

Those are the same warnings I get when I test the package (before submitting to CRAN) and have been that way for a long time. They stemmed from conditional allocation of arrays in C. gcc -wall seems to always pick on that. As far as I know, they are harmless. Andy -Original

Re: [R] RandomForest - getTree status code

Is that the entire tree? If so there's a problem. The node status is defined as follows in rf.h of the source code: #define NODE_TERMINAL -1 #define NODE_TOSPLIT -2 #define NODE_INTERIOR -3 i.e., -3 means non-terminal node. Andy -Original Message- From:

[R] outputting functions in lapply

How come the k is 3 in all of this output? I expected it to be equal to r. tmp3 - lapply(1:3, function(k) function(r) print(paste(r, - r | k -, k))) for (i in 1:3) { tmp3[[i]](i) } [1] 1 - r | k - 3 [1] 2 - r | k - 3 [1] 3 - r | k - 3 -- Med venlig hilsen Rune Schjellerup Philosof

Re: [R] optim with constraints

Without the data / script, I'm guessing that it is likely an attempt to evaluate the loss function at an inadmissible point e.g., at the constraint where there is a log(0). Different optimization tools handle things differently, and there are a couple of us working (very slowly due to other

Re: [R] Passing arguments to gpar

Hi, Yes, that is exactly it. Charlie Sharpsteen gave me the same solution last week but he probably just replied to me, so his email did not go to the list. Thanks for the reply, I appreciate it I always forget about this do.call function. Paul Murrell wrote: Hi Sebastien Bihorel

Re: [R] savePlot for Mac and / or Linux?

Hi Christophe, No, there isn't anything close to savePlot() on either Linux or MacOS. But the difference is that savePlot() is used to save something that has already been plotted. Since your function appears to be intended to give the user the ability to save something that has yet to be

Re: [R] outputting functions in lapply

On 07/12/2009 8:53 AM, Rune Schjellerup Philosof wrote: How come the k is 3 in all of this output? I expected it to be equal to r. tmp3 - lapply(1:3, function(k) function(r) print(paste(r, - r | k -, k))) for (i in 1:3) { tmp3[[i]](i) } [1] 1 - r | k - 3 [1] 2 - r | k - 3 [1] 3 - r | k - 3

Re: [R] optim with constraints

Steven, You haven't given us much information to help you. Are you `maximizing' the log-likelihood or are you `minimizing' the negative log-likelihood? If it is the former, you have to specify the control argument, control$fnscale = -1. This is a common mistake when using optim for likelihood

Re: [R] outputting functions in lapply

Hi, I was about to send the exact same answer as you just received, so I'll add a note instead. Your problem looks a bit like Currying, Curry - # original from roxygen function (f, ..., .left=TRUE) { .orig = list(...) function(...){ if(.left) {args - c(.orig, list(...))} else

Re: [R] Forest Plot

If you just want a forest plot, then the forest() function. If you have the betas and corresponding variances, then you can create a forest plot with: forest(betas, varbetas) And yes, estimate +/- 1.96*sqrt(variance of estimate) would be an *approximate* 95% CI. Best, -- Wolfgang

[R] Subset sum problem.

Hi, I'm quite new to the R-project. I was suggested to look into it because I am trying to solve the Subset sum problem, which basically is: Given a set of integers and an integer s, does any non-empty subset sum to s? (See http://en.wikipedia.org/wiki/Subset_sum_problem) I have been searching

Re: [R] Subset sum problem.

The problem is NP-Complete and the real problem is how you can solve it. According to the wiki page, you can use the bottom algorithm(Polynomial time approximate algorithm) to solve your problems. If you had trouble with writing R code, you can read ``R-introduction''. regards

[R] How to apply five lines of code to ten dataframes?

Hello R-helpers, I have 10 dataframes (named data1, data2, ... data10) and I would like to add 5 new columns to each dataframe using the following code: data1$LogDepth-log10(data1[,2]/data1[,4]) data1$LogArea-log10(data1[,3]/data1[,5]) data1$p-2*data1[,6]/data1[,7]

Re: [R] How to apply five lines of code to ten dataframes?

Hi, You could define a function that does the calculations for a given data.frame and apply it to all your data.frames, d1 - data.frame(a=1:10, b=rnorm(10)) d2 - data.frame(a=-(1:10), b=rnorm(10)) calculations - function(d){ if(is.character(d)) d - get(d) transform(d, c =

Re: [R] How to apply five lines of code to ten dataframes?

For example, tt1 = 2:4 ind = 1 evalstr = paste(xx = tt, ind, sep = ) eval(parse(text=evalstr)) == x = 2, 3, 4 It's easy for you to write a for loop to solve your problems Guo-Hao Huang

[R] A Gamma-GLM with log link

Hi, I have a set of data (total number of record = 144,122), and I would like to use gamma-glm with log link to set up a model. IC is number of records IL is paid amount The table below shows that I have 30.578% of the data in the level of 1 - 1000 paid amount 20.320% of the data in the level

Re: [R] SAS datalines or cards statement equivalent in R?

Regarding the various methods people have suggested, what if a typical tab-delimited data line looks like: SMS11001 1990 M01 688.0 and the SAS INPUT statement is INPUT survey $ 1-2 seasonal $ 3 state $ 4-5 area $ 6-10 supersector $ 11-12 @13 industry $8. datatype $ 21-22

[R] Subset of time observations where timediff 60 secs

Dear list members I have a rather large vector (part of a data frame) giving the time (date + time, POSIXct) of observations. The times are irregular (with both small and large jumps) but increasing, and there are several millions of them. I now wish to reduce my data set, so that I only have

Re: [R] latex.table for table with character and numeric columns

Ha Charlie, It helped. It is just what i was looking for. Gerrit. Gerrit Draisma wrote: Hallo, I have a dataset with one or two columns with character data and the rest with numeric data. Using latex.table from the quantreg package produced a table, but I cannot set the decimals. For

[R] different model performance because of nlevels()???

Hello everybody, I came across a problem when building a randomForest model. Maybe someone can help me. I have a training- and a testdataset with a discrete response and ten predictors (numeric and factor variables). The two datasets are similar in terms of number of predictor, name of

[R] multiple plots using summary in rms package

Dear All, I wonder if someone can point me in the right direction here. I'm working with the rms library, R 2.9.2 under Windows XP. I'm trying to arrange two plots side by side for a colleague. mfrow or mfcol do not seem to work, however, so I am obviously missing something important. I know

Re: [R] [ggplot2] Wind rose orientation

Hi Hadley, Thank you very much! The idea of plotting a wind rose must be fairly common. I wonder if it would make sense to have a switch that would wrap data around the ends of a continuous scale? All the best, Tom On Dec 6, 2009, at 7:51 PM, hadley wickham wrote: Hi Tom, The problem

Re: [R] optim with constraints

Steven, A quick trick to deal with problems of this nature (i.e. forcing parameters to stay positive although the parameter search isn't bound to non-negatives or positives) without resorting to another package is to re-parametrize the density function. Notice, x \in (-\infty,\infty) for x \in

Re: [R] [ggplot2] Wind rose orientation

The idea of plotting a wind rose must be fairly common.  I wonder if it would make sense to have a switch that would wrap data around the ends of a continuous scale? Probably - but it requires a lot of work, because ggplot2 doesn't currently support circular scales, which is what you really

[R] anova/factor

Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with thanks in advance. Group SI Sif SG

Re: [R] anova/factor

On Dec 7, 2009, at 11:38 AM, oscar linares wrote: Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with

Re: [R] SAS datalines or cards statement equivalent in R?

On Mon, Dec 7, 2009 at 3:53 PM, Marshall Feldman ma...@uri.edu wrote: Regarding the various methods people have suggested, what if a typical tab-delimited data line looks like:     SMS11001 1990 M01 688.0 and the SAS INPUT statement is   INPUT survey $ 1-2 seasonal $ 3 state $

[R] [GGPLOT] Legends at different layers

Here I have following code : dat = rnorm(100) ggplot() + geom_histogram(aes(dat, fill=..count..)) + scale_fill_gradient(Count, low=green, high=red) + opts(legend.position=none) # Above is without any legend ## Now I want to place two points with legends dat1 - data.frame(c(0,0), c(1,0)); Label

Re: [R] SAS datalines or cards statement equivalent in R?

On Mon, Dec 7, 2009 at 11:43 AM, Barry Rowlingson b.rowling...@lancaster.ac.uk wrote:  This just seems like horrendously bad practice, which is one reason it's kludgy in R. If it was good practice, someone would surely have written a way to do it neatly.  Keep your data in data files, and

[R] Setting Repositories in the .Rprofile file

Is there any way to select (during startup) all repositories by modifying the Rprofile.site file? If so, could someone let me know the syntax to use (or direct me to where I can find it)? Thank you, Patrick This email message, including any attachments, is for th...{{dropped:9}}

Re: [R] Subset of time observations where timediff 60 secs

There is an example at the end of the Prices and Returns section of the zoo-quickref vignette in the zoo package. library(zoo) vignette(zoo-quickref) If speed is your main concern check this recent thread that was posted on R-sig-finance:

Re: [R] SAS datalines or cards statement equivalent in R?

I totally agree with Barry, although it's sometimes convenient to include data with analysis code for debugging and/or documentation purposes. However, the example actually applies equally to separate data files. In fact, the example is from the U.S. Bureau of Labor Statistics at

Re: [R] Setting Repositories in the .Rprofile file

On Dec 7, 2009, at 12:18 PM, Richardson, Patrick wrote: Is there any way to select (during startup) all repositories by modifying the Rprofile.site file? If so, could someone let me know the syntax to use (or direct me to where I can find it)? ?options options(repos=c(CRAN =

[R] Output multiple files

Hello, I want the ouput to take the name of the original file. I have the following code - but the problem is that the output overrides the original file which I desire not. So given an input file called File1.txt, I want the ouput to be called File1.out.txt, or something along this line.

Re: [R] SAS datalines or cards statement equivalent in R?

[stuff about SAS CARDS statement, I think] INPUT survey $ 1-2 seasonal $ 3 state $ 4-5 area $ 6-10 supersector $ 11-12 @13 industry $8. datatype $ 21-22 year period $ value footnote $ ; ... This just seems like horrendously bad practice, which is one reason it's kludgy in R. If it

Re: [R] SAS datalines or cards statement equivalent in R?

On Mon, Dec 7, 2009 at 5:37 PM, Marshall Feldman ma...@uri.edu wrote: I totally agree with Barry, although it's sometimes convenient to include data with analysis code for debugging and/or documentation purposes. However, the example actually applies equally to separate data files. In fact,

[R] Multiple regression script

Hi all, I'm doing Multiple linear regression for a data set. However, it takes a lot of time, as I would like to check every possible combination of factors, evalute the results based for instance on their p values, and then choose the best regression model. So, I wonder if anyone might have

Re: [R] SAS datalines or cards statement equivalent in R?

strapply in the gsubfn package can parse any input format that can be described using a regular expression and generally requires only one line of code to do the parsing. See: http://gsubfn.googlecode.com On Mon, Dec 7, 2009 at 12:37 PM, Marshall Feldman ma...@uri.edu wrote: I totally agree

Re: [R] SAS datalines or cards statement equivalent in R?

On Dec 7, 2009, at 12:37 PM, Marshall Feldman wrote: I totally agree with Barry, although it's sometimes convenient to include data with analysis code for debugging and/or documentation purposes. However, the example actually applies equally to separate data files. In fact, the example

Re: [R] Output multiple files

On Dec 7, 2009, at 12:42 PM, Muhammad Rahiz wrote: Hello, I want the ouput to take the name of the original file. I have the following code - but the problem is that the output overrides the original file which I desire not. So given an input file called File1.txt, I want the ouput to

Re: [R] [ggplot2] Wind rose orientation

On Dec 7, 2009, at 7:00 AM, hadley wickham wrote: The idea of plotting a wind rose must be fairly common. I wonder if it would make sense to have a switch that would wrap data around the ends of a continuous scale? Probably - but it requires a lot of work, because ggplot2 doesn't

[R] Very basic R workflow question for Windows users

Hi all, If you use the R Editor (not another text editor), please read on... So, my usual R workflow involves having two windows open (R Console, R Editor), writing a line of code in the R Editor and then typing Ctrl-R to run that line. Then, quite frequently, I want to run a new command to

Re: [R] [Ann] R brush for SyntaxHighlighter

Hi, On Dec 5, 2009, at 2:43 AM, Romain Francois wrote: There is also the R package highlight [1,2] serving the same purpose but based on information from the R parser, not a static list. Thanks for the pointer - I wasn't aware of that package. I'll make a mention of it on my web page.

Re: [R] yeroon.net/ggplot2 web application v0.11

Hi Jeroen: This is great for someone who is point and click oriented. Can this program be downloaded to be used offline or its just web based? My field crew usually take toughbooks to the field and since they don't know how to program it will be nice to make the graphics just by pointing and

[R] Sweave and license message when loading mclust package

When loading mclust, it shows a license agreement message. This message shows up in my document when I use Sweave. I did the following: echo=FALSE,include=FALSE= library(mclust) @ Is this a problem with mclust, Sweave or with me? How can it be fixed? Thanks for any suggestions! Titus

Re: [R] Sweave and license message when loading mclust package

On Mon, 7 Dec 2009, Titus von der Malsburg wrote: When loading mclust, it shows a license agreement message. This message shows up in my document when I use Sweave. I did the following: echo=FALSE,include=FALSE= library(mclust) @ Is this a problem with mclust, Sweave or with me? How can

Re: [R] Very basic R workflow question for Windows users

On 07/12/2009 1:31 PM, Mark Na wrote: Hi all, If you use the R Editor (not another text editor), please read on... So, my usual R workflow involves having two windows open (R Console, R Editor), writing a line of code in the R Editor and then typing Ctrl-R to run that line. Then, quite

Re: [R] SAS datalines or cards statement equivalent in R?

On 8/12/2009, at 5:43 AM, Barry Rowlingson wrote: On Mon, Dec 7, 2009 at 3:53 PM, Marshall Feldman ma...@uri.edu wrote: Regarding the various methods people have suggested, what if a typical tab-delimited data line looks like: SMS11001 1990 M01 688.0 and the SAS INPUT

[R] Convert a list of N dataframes to N dataframes

Hello, I have used the following command: datalist-list(data1,data2,data3,data4,data5,data6) to make a list of my dataframes, which I then manipulated with these commands: datalist-llply(datalist,LogDepth) datalist-llply(datalist,LogArea) datalist-llply(datalist,p) datalist-llply(datalist,Exp)

Re: [R] [GGPLOT] Legends at different layers

You mean: dat - data.frame(x = rnorm(100)) dat1 - data.frame( x = c(0,0), y = c(1,0), Label = c(Point1, Point2) ) ggplot(dat, aes(x)) + geom_histogram(aes(fill = ..count..)) + geom_point(aes(x, y, colour = Label), data = dat1, size = 4) + scale_fill_gradient(Count, low = green, high

[R] Are there free R webinar recordings somewhere ?

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

Re: [R] Are there free R webinar recordings somewhere ?

Tal, You might look here: http://www.fort.usgs.gov/brdscience/LearnR.htm I just googled r stats webinar Tom Tal Galili wrote: Hi all, A friend just sent me this: http://www.mathworks.com/company/events/webinars/index.html?id=language=en

[R] Reminder: ASA Stat Comp/Graph Student Paper Competition Deadline 12/15

Just a gentle reminder that the deadline for submission to the Student Paper Competition is fast approaching. All application materials must be received by 5:00 PM EST, Monday, December 14, 2009. For submission guidelines, please visit http://stat-computing.org/awards/student/announcement.html

[R] Reminder: ASA Stat Comp/Graph Student Paper Competition Deadline 12/14

Just a gentle reminder that the deadline for submission to the Student Paper Competition is fast approaching. All application materials must be received by 5:00 PM EST, Monday, December 14, 2009. For submission guidelines, please visit

Re: [R] SAS datalines or cards statement equivalent in R?

Barry Rowlingson wrote: I'd love to duplicate this functionality of SAS, however, I fear: http://www.sas.com/news/preleases/SASsuit.html Amazing, since input statements in SAS bear an uncanny resemblance to how PL/I handles input from text files.

Re: [R] MLE in R

Hello Liang, Besides looking at ?optim, you may also want to look into this nice working example www.mayin.org/ajayshah/KB/R/documents/mle/mle.html Regards, Francisco Francisco J. Zagmutt Vose Consulting 1643 Spruce St., Boulder Boulder, CO, 80302 USA www.voseconsulting.com Liang Wang

Re: [R] SAS datalines or cards statement equivalent in R?

On Mon, Dec 7, 2009 at 8:54 PM, Marshall Feldman ma...@uri.edu wrote: Barry Rowlingson wrote:  I'd love to duplicate this functionality of SAS, however, I fear: http://www.sas.com/news/preleases/SASsuit.html Amazing, since input statements in SAS bear an uncanny resemblance to how PL/I

Re: [R] Are there free R webinar recordings somewhere ?

Tal I have a series of basic introduction tutorials for new R users at this link. http://processtrends.com/Learn_R_Toolkit.htm http://processtrends.com/Learn_R_Toolkit.htm These videos are geared to Excel users who want to learn about R and see it in action. D Kelly O'Day

Re: [R] SAS datalines or cards statement equivalent in R?

On Dec 7, 2009, at 12:37 PM, Marshall Feldman wrote: I totally agree with Barry, although it's sometimes convenient to include data with analysis code for debugging and/or documentation purposes. However, the example actually applies equally to separate data files. In fact, the example

Re: [R] bootstrap help

If I read your question now, I think, that I didn't understand you correctly, but... My code generates 50 bootstrap samples (50 tables) and for each calculates mean of column x and mean of column y, so the final result is matrix with means. If you don't want to calculate means, code below

[R] Help required with NLME syntax

Hi there, I am attempting to analyse data from a crossover study using a mixed effect Emax model in R using nlme. I am having problems ascertaining how to fit the model. Subjects are randomly assigned to 1 of 4 treatments in a complete block design. The response variable of interest is

Re: [R] yeroon.net/ggplot2 web application v0.11

Hi Felipe, my apps are completely webbased, so you cannot use them offline, unless you are hosting a mirror locally or in your LAN (which can be arranged if you are interested). Jeroen 2009/12/7 Felipe Carrillo mazatlanmex...@yahoo.com: Hi Jeroen: This is great for someone who is point and

Re: [R] SAS

Marshall Feldman marsh at uri.edu writes: Regarding the various methods people have suggested, what if a typical tab-delimited data line looks like: SMS11001 1990 M01 688.0 and the SAS INPUT statement is INPUT survey $ 1-2 seasonal $ 3 state $ 4-5 area $ 6-10

[R] сЛyЖБа PAcCЫЛoK

E-MAIL РАССЫЛКИ Доступные цены - высокое качество Наш прайс-лист: Рассылка по всей Москве - 5 000 руб., 4 рассылки по Всей Москве - 10 000 руб. Рассылка по всей России - 10 000 руб., 4 рассылки по всей России - 20 000 руб. Рассылка по России и СНГ - 15 000 руб. , 4 рассылки по России и СНГ - 30

[R] zelig logit survey

When I use zelig to compute parameter estimates with model = logit.survey, I receive the following error: Nicholas Carnes. 2007. logt.surveyWarning message: In eval(expr, envir, enclose) : non-integer #successes in a binomial glm! I believe this is because the model is not using quasibinomial

[R] problem in labeling the nodes of tree drawn by rpart

Hi all, I used rpart to fit a model, where the covariates in the model are categorical variables. Then I plotted the tree (mytree) and used the command text to add labels to the tree. In the nodes of the tree, the values of the covariates are represented with a, b or c. Is there a way to show

Re: [R] Label individual points in lattice by and ID

Thanks This works, but I lose the ability to say plot a line by group connecting the numbered points, how do I recover this, or say color the groups by ID.. xyplot(y~x | f1 + f2, data= dat, groups=ID, type=l, panel= function(x, y,groups, subscripts, ...){ panel.text(x, y, groups[subscripts])

Re: [R] Announce: edtdbg, integrating R's debug() with your text editor

The files seem to be inaccessible, and I can't see why. Neither can my department Web technician. So, for the time being, simply go to http://heather.cs.ucdavis.edu/~matloff/R/edtdbg/ and then download the 3 files as directed. Meanwhile, I've thought of an interesting enhancement, and a

Re: [R] RandomForest - getTree status code

Not the entire tree - I just selected a few rows. On the first tree, the only status codes I saw were -1 and -3, which makes sense. I am guessing the final tree should never have a status of -2, which looks like an intermediate status, unless I am missing something. Thanks for the response.

[R] Regular expression help

Hi there I have a string like this i want to extract 9831019 from this string i used a regular expresion \d+ by which i can only make it to see 7 and returns. This type of number(9831019) appears in any part of the string and is definitely more than 5 digits all the time and i want to give

Re: [R] data.frame to stacked frame

One option is library(reshape) dfr - melt(dfw, id.vars=group) -Ista On Mon, Dec 7, 2009 at 8:20 AM, Christian Schulz chsch...@email.de wrote: Hi, have anybody a hint how i could avoid the cumbersome way (..especially the value part) change the data representation from  n column data.frame

Re: [R] Regular expression help

On Dec 7, 2009, at 5:04 PM, Ramya wrote: Hi there I have a string like this i want to extract 9831019 from this string i used a regular expresion \d+ by which i can only make it to see 7 and returns. This type of number(9831019) appears in any part of the string and is definitely more

Re: [R] Regular expression help

Ramya - Try strings = c('UV7C11-F9-E1 MCS#9831019','MCS Lot #9512516') sub('^.*?(\\d{5,}).*?$','\\1',strings,perl=TRUE) [1] 9831019 9512516 The regular expression finds the first string of five or more numbers in the strings. Since you said the numbers could occur anywhere in the string,

Re: [R] Regular expression help

If I understand correctly you wish to extract strings of digits more than 5 characters long: s - c(UV7C11-F9-E1 MCS#9831019, MCS Lot #9512516) library(gsubfn) strapply(s, \\d{6,}, c) Depending on what you want to get back you might wish to add the simplify=TRUE argument to strapply, as well.

  1   2   >