Re: [R] MySQL + R as a Replacement for SAS Proc SQL + Various Stat Procs

Hi, there, it seems size of your data is very manageable in R. i've been using proc sql in SAS a lot and found sqldf package is very similar to proc sql. check it out. one more thought. instead of using R as scripting lang. to process data, it might be worthy to take a look at python. On Jan 18, 2

Re: [R] MySQL + R as a Replacement for SAS Proc SQL + Various Stat Procs

On Fri, 18 Jan 2008, [EMAIL PROTECTED] wrote: > For cost reasons, I'd like to replace SAS on my PC under Win XP Pro. > > Nearly all my work involves medium-size datasets (100k-10M) records which I > cleanup, relate, fliter and get into shape for analysis using SAS/SQL > followed by standard statis

Re: [R] How do we get two-tailed p-values for rlm?

On Fri, 18 Jan 2008, Darren Weber wrote: > How do we get 2-tailed p-values for the rlm summary? > > I'm using the following: > >> fit <- rlm(oatRT ~ oatoacData$erp, psi=psi.bisquare, maxit=100, > na.action='na.omit') >> fitsum <- summary(fit, cor=F) >> print(fitsum) > > Call: rlm(formula = o

[R] MySQL + R as a Replacement for SAS Proc SQL + Various Stat Procs

For cost reasons, I'd like to replace SAS on my PC under Win XP Pro. Nearly all my work involves medium-size datasets (100k-10M) records which I cleanup, relate, fliter and get into shape for analysis using SAS/SQL followed by standard statistical procedures, e.g. regression using SAS proc reg.

Re: [R] how to use different variable to store values with different length in a loop

I would recommend that you use a "list" instead of creating the variables. Here is how you might do a list: # preallocate the list myList <- vector('list', 100) for (i in 1:100){ .you calculations. myList[[i]] <- result } If you really want to create 100 variables, then use 'assi

[R] how to use different variable to store values with different length in a loop

Hi, List I am trying to use variables named A1, A2, ..., A100 to store some values, each variable could store some values with different length, how can I achieve this? Thanks, Jack [[alternative HTML version deleted]] __ R-help@r-project.org

[R] How do we get two-tailed p-values for rlm?

How do we get 2-tailed p-values for the rlm summary? I'm using the following: > fit <- rlm(oatRT ~ oatoacData$erp, psi=psi.bisquare, maxit=100, na.action='na.omit') > fitsum <- summary(fit, cor=F) > print(fitsum) Call: rlm(formula = oatRT ~ oatoacData$erp, psi = psi.bisquare, maxit = 100,

[R] Serez -vous notre prochain gagnant?

-- Pour vous d�sabonner � cette liste, visitez http://jbala.ws/lists/lt.php?id=MkxTBQYBCFVIAU8AVV1XAQ%3D%3D Pour mettre vos pr�f�rences � jour, visitez http://jbala.ws/lists/lt.php?id=MkxTBQYBCFRIAU8AVV1XAQ%3D%3D Faire suivre un message a quelqun http://jbala.ws/lists/lt.php?id=MkxTBQYBCFdIA

Newsletter: Vielen dank f�r ihre Anmeldung!

Hallo, Ihre E-Mailadresse "[EMAIL PROTECTED]" wurde in den Newsletterverteiler hinzugefügt. URL: http://www.fueralles.de/cgi-bin/nsl_pro/newsletter.cgi?id=steinfeuerfred Mit freundlichen Grüssen Ihr Newsletter-Team __ R-help@r-project.org mailing l

Re: [R] Adding a table to a plot

Here is the brute force method. It would probably not be very difficult to turn in into a function. In fact I am pretty sure that someone has but I cannot find what package it is in. x <- "Time SD Mean 12.5 30 2a.3 40" xx <- read.table(textConnection(x), header=TRUE, as.is=T

Re: [R] Fortran level mismatch between package and R

Hi, I recompiled R using gfortran in my personal laptop. But I think it may help if you build an agent C program which calls the FORTRAN program, and then call the C program from R. I'm also curious on how to solve this problem without too much work. Please let me know if you find a good solution

[R] Newsletter - Aktivierungslink

Sehr geehrte(r) Abonnent(in), Ihre E-Mailadresse "[EMAIL PROTECTED]" wurde in unseren Newsletter eingetragen. Um die Anmeldung abzuschliessen, folgen Sie dem Link unterhalb. Dieses Verfahren schützt Sie vor ungewollten Mails von unserer Seite aus. Wir bitten Sie daher um Verständins. Sollten Si

[R] modifing barplot

Hi, I would like to modify how barplot plots. For reasons only my supervisor comprehends i need to do the crisscross lines in the bar thicker. If i use lwd=2 the y-axis is thicker, (i use the option beside = TRUE) but nothing else is affected by this parameter - it seems. I looked at th

Re: [R] things that are difficult/impossible to do in SAS or SPSSbut simple in R

hadley wickham wrote: > On Jan 18, 2008 1:19 PM, Jeffrey J. Hallman <[EMAIL PROTECTED]> wrote: >> Frank E Harrell Jr <[EMAIL PROTECTED]> writes: >> >>> Rob Robinson wrote: I wonder if those who complain about SAS as a programming environment have discovered SAS/IML which provides a progra

Re: [R] Regex magic anyone?

Great! Joh On Friday 18 January 2008 15:22:29 Gabor Grothendieck wrote: > strapply in the gsubfn package can split strings based on > content rather than delimiters. > > > library(gsubfn) > > strapply(s, ".[(].[)]|.")[[1]] > > [1] "A""B""C""D""E""F""G""T(P)" "H""I

Re: [R] html search not working

On 1/18/2008 3:07 PM, [EMAIL PROTECTED] wrote: > Sorry the first email had no subject. > I've searched around and havent had much luck trying to explain why my > HTML search no longer works. If I use > ?cor.test > it is unable to display the webpage. Is this a java problem? I am > running window

[R] How to write a Greek letter in a graph?

Hi, Sorry for bothering you again. I have a small question about the graph and the following is my code: xyplot(power1+power2 ~ pi0, data=resultb, ylab="Average Power", type="l", pch=LETTERS[1:2], lty=c(1:2), main="(a)") The name of the X-axis will be “pi0” in the graph. What should I do if

Re: [R] things that are difficult/impossible to do in SAS or SPSSbut simple in R

On Jan 18, 2008 1:19 PM, Jeffrey J. Hallman <[EMAIL PROTECTED]> wrote: > Frank E Harrell Jr <[EMAIL PROTECTED]> writes: > > > Rob Robinson wrote: > >> I wonder if those who complain about SAS as a programming environment have > >> discovered SAS/IML which provides a programming environment akin to

[R] html search not working

Sorry the first email had no subject. I've searched around and havent had much luck trying to explain why my HTML search no longer works. If I use ?cor.test it is unable to display the webpage. Is this a java problem? I am running windows xp, and use Internet explorer for my default browser. Tha

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

Jeffrey J. Hallman wrote: > Frank E Harrell Jr <[EMAIL PROTECTED]> writes: > >> Rob Robinson wrote: >>> I wonder if those who complain about SAS as a programming environment have >>> discovered SAS/IML which provides a programming environment akin to Matlab >>> which is more than capable (at least

Re: [R] things that are difficult/impossible to do in SAS or SPSSbut simple in R

Frank E Harrell Jr <[EMAIL PROTECTED]> writes: > Rob Robinson wrote: >> I wonder if those who complain about SAS as a programming environment have >> discovered SAS/IML which provides a programming environment akin to Matlab >> which is more than capable (at least for those problems which can be t

Re: [R] Adding a table to a plot

Judith Flores wrote: > Hi R-experts, > >Is it possible to add a table (text and lines) to > the plot area? > > I have a data frame that looks something like this: > >> mydata > > Time SD Mean >12.5 30 >2a.3 40 > > > I am using R 2.6.1, in Windows. > > Thank y

Re: [R] re siduals from pcaiv

Hi Nikolaos, >> How can I extract the residuals from a pcaiv/rda in ade4? >From what you can remember, you can get them as follows (ade4): ## If your call to pcaiv was as follows: pcaiv(dudi$obj, df=FactorIV, scannf=F, nf=4) ## then get your residuals as follows: ortho <- pcaivortho(dudi$obj,

Re: [R] Selecting rows conditionally between 2 data.frames

Looks like a job for ?merge HTH, Chuck On Fri, 18 Jan 2008, jiho wrote: > Hello everyone, > > I have two data.frames that look like > > calib: > place zoomscale > left 0.658 > left 0.805.6 > left 1.203 > right

[R] Adding a table to a plot

Hi R-experts, Is it possible to add a table (text and lines) to the plot area? I have a data frame that looks something like this: >mydata Time SD Mean 12.5 30 2a.3 40 I am using R 2.6.1, in Windows. Thank you, Judith Flores _

[R] forming a linear discriminant function from the output of lda()

Hello all- I am a relatively new user of R and am working through a graduate course in Statistics that uses Minitab, SAS and some Matlab. I like using R but am having some trouble lining up the output from lda() to that of the other programs' results. The dataset below is a modified set of wine

[R] (no subject)

I've searched around and havent had much luck trying to explain why my HTML search no longer works. If I use ?cor.test it is unable to display the webpage. Is this a java problem? I am running windows xp. Thanks -- David Kowalske Jr Graduate Rese

Re: [R] plotting other axes for PCA

Hi Silvia, >> What I need is exactly what I get using biplot (pca.object) but for other >> axes. You need to look at ?biplot.prcomp (argument: choices=) ## Try biplot(prcomp(USArrests), choices=c(1,2)) ## plots ax1 and ax2 biplot(prcomp(USArrests), choices=c(1,3)) ## plots ax1 and ax3

[R] more help needed Re: communicate from Rterm

I think the message from Prof. Ripley refers to calling shell and system from R. There is no "system" call in Windows cmd or VBA and I am calling R from VBA. Now I want to use tryCatch to output the error to a text file and exit gracefully. My expression will be source("file.r") and that will

Re: [R] image/area plot

Yes, that is it, a square pie chart :-) I did not knew the name... sorry... Does anyone knows about it? Thank you very much, Best wishes, Marta hadley wickham wrote: > Do you have an example graphic that shows what you're trying to > create? I can't figure out if you want something like a square

Re: [R] plotting other axes for PCA

On Fri, 18 Jan 2008, Silvia Lomascolo wrote: > > Hi R-community, > I am doing a PCA and I need plots for different combinations of axes (e.g., > PC1 vs PC3, and PC2 vs PC3) with the arrows indicating the loadings of each > variables. What I need is exactly what I get using biplot (pca.object) but

[R] Amusez vous et choisissez vos cadeaux!

-- Pour vous d�sabonner � cette liste, visitez http://jbala.ws/lists/lt.php?id=MkxaBQUFAB0FTgEDXVVW Pour mettre vos pr�f�rences � jour, visitez http://jbala.ws/lists/lt.php?id=MkxaBQUFAR0FTgEDXVVW Faire suivre un message a quelqun http:

[R] Fortran level mismatch between package and R

We have a user here who is attempting to compile a package which requires f90 (gfortran will do). Our version of R was built with f77 (g77). When the user attempts to build the package, f77 is always invoked and fails due to the presence of f90 features. Is there a simple method or some docu

Re: [R] image/area plot

Do you have an example graphic that shows what you're trying to create? I can't figure out if you want something like a square pie chart (aka waffle chart), a stacked barchart, a levelplot, or something else. Hadley On Jan 18, 2008 6:06 AM, Marta Rufino <[EMAIL PROTECTED]> wrote: > Dear R users,

Re: [R] histogram with NAs

> "Typical plots with vertical bars are not histograms. Consider barplot > or plot(*, type = "h") for such bar plots. ". But no worry, I've mixed > them up myself a number of times. Or you can use ggplot2, which will do the right thing regardless of whether you have continuous or categorical data:

Re: [R] constrOptim with method SANN

Hi, thanks again for the clarification. However, then the documentation of constrOptim is somehow misleading, because it says "Any optim method that permits infinite values for the objective function may be used (currently all but "L-BFGS-B")." But I assume that this (documentation issues) is n

[R] plotting other axes for PCA

Hi R-community, I am doing a PCA and I need plots for different combinations of axes (e.g., PC1 vs PC3, and PC2 vs PC3) with the arrows indicating the loadings of each variables. What I need is exactly what I get using biplot (pca.object) but for other axes. I have plotted PC2 and 3 using the sc

Re: [R] gboxplot (JMP Diamond plot in R?)

On Fri, 18 Jan 2008, Richard Yanicky wrote: > I am looking to replicate the functionality of the JMP diamond plot in > R. I have found the following enhancement and no other references(see > below). I have tried to contact the author via email and I have had no > luck. > > Does any know if this

[R] Selecting rows conditionally between 2 data.frames

Hello everyone, I have two data.frames that look like calib: place zoomscale left0.658 left0.805.6 left1.203 right 0.658.4 right 0.806 right 1.202.9 X:

Re: [R] constrOptim with method SANN

On Fri, 18 Jan 2008, Jonas Rumpf wrote: > Hi, > > thank you for your help, Thomas, > (and my apologies to everyone for the double-post earlier), > however I'm not sure if I understand you correctly: > > Are you basically saying that it doesn't really make sense > to use SANN with constrOptim and I

Re: [R] histogram with NAs

On Jan 18, 2008 4:49 PM, <[EMAIL PROTECTED]> wrote: > Dear list, > > I have a categorical variable in a data.frame that I would like to > plot using a histogram to show number of events. Values are 0, 1 and > some NAs. I can´t make the hist() function to > 1) include a column with the number of NA

Re: [R] constrOptim with method SANN

Hi, thank you for your help, Thomas, (and my apologies to everyone for the double-post earlier), however I'm not sure if I understand you correctly: Are you basically saying that it doesn't really make sense to use SANN with constrOptim and I should use another algorithm (like CG etc.) ? Best re

[R] gboxplot (JMP Diamond plot in R?)

I am looking to replicate the functionality of the JMP diamond plot in R. I have found the following enhancement and no other references(see below). I have tried to contact the author via email and I have had no luck. Does any know if this functionality exists in R? If it does not does anyone

Re: [R] constrOptim with method SANN

On Fri, 18 Jan 2008, Jonas Rumpf wrote: > Hi Everyone, > > I'm trying to minimize a function using constrOptim with > the simulated annealing method SANN. > > If I understand constrOptim well, it basically passes most > of its arguments to optim while somehow enforcing the constraints. > > My prob

Re: [R] Trouble receiving messages from the mailing list

On Friday 18 January 2008, Ted Harding wrote: > [...] > > This policy would carry a risk that spam mail "From:" a > subscribed address (probably with a forged "From:", but > possibly also from a subscriber's compromised computer) > will reach the list. But I dare say the rest of us can > live with

Re: [R] Converting plots to ggplot2

Dear David, Try this solution. The data are stored is a dataframe instead of a list but they should be the same. Once you have this dataframe, thinks get rather easy. library(ggplot2) a <- seq(0, 360, 5)*pi/180 s <- seq(5, 45, 10)*pi/180 dataset <- expand.grid(a = a, s = s) d

[R] histogram with NAs

Dear list, I have a categorical variable in a data.frame that I would like to plot using a histogram to show number of events. Values are 0, 1 and some NAs. I can´t make the hist() function to 1) include a column with the number of NAs 2) have the x axis to be categorical, I always get 0, 0.2,

Re: [R] Assigning into each of a list of dataframes

The problem is that your function is returning x$Y rather than x so the result is a list of Y's rather than a list of the x's with Y added. You want function(x) { x$Y <- x$A * x$B; x } or a bit shorter: lapply(mylist, transform, Y = A * B) On Jan 18, 2008 10:37 AM, Chuck Cleland <[EMAIL PROTE

Re: [R] Converting plots to ggplot2

Hello Hadley, Thank you (and Domenico) for your help. I'm almost there. For convenience, I'll restate the question. With the following data: a <- seq(0, 360, 5)*pi/180 ; a ac <- sin(a + (45*pi/180)) + 1 ; ac s <- seq(5, 45, 10)*pi/180 ; s asc <- lapply(s, function(x) x*cos(ac) + x*

[R] Assigning into each of a list of dataframes

What is the right way to assign a new variable into each a of list of data frames? Here is my failed attempt: mylist <- list(df1 = data.frame(A = runif(5), B = runif(5)), df2 = data.frame(A = runif(5), B= runif(5))) lapply(mylist, function(x){x$Y <- x$A * x$B}) $df1 [1] 0.25

[R] residuals from pcaiv

Dear R users, How can I extract the residuals from a pcaiv/rda in ade4? In Vegan there is the residuals() function, giving the approximation of the original data from the unconstrained ordination Is there something similar in ade4? Nikos __ R-help@r-p

Re: [R] Question about making active a previously opened graphic device

Perhaps: dev.set(3) On 18/01/2008, Monica Pisica <[EMAIL PROTECTED]> wrote: > > Hi, > > I have opened 2 different graphic devices on screen with windows(3) and > Windows(4) (I already have had an other graphic on screen opened in device 2). > > Now device 3 is inactive and device 4 is active. Ho

[R] FW: Question about making active a previously opened graphic device - SOLVED

Hi again, It seems that dev.set() will solve my problem. It took me quite a bit of time, search and frustration until i stumbled on dev.set ;-) Thanks anyway, Monica From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: [R] Question about making active a previously opened graphic deviceDate: Fri,

[R] Res: Leyenda (revisited translation)

(just tranlating...) (i made a little mistake on last on / miltinho). Hi, I am trying to put the legend on the outside of plot region. How can I do that? My code looks like... df<-read.table(stdin(),head=T,sep=",") row,LENGTH,LAT 091639,10.002,42.26282 091640,30.808,42.26834 091641,21.591,42.3

[R] Question about making active a previously opened graphic device

Hi, I have opened 2 different graphic devices on screen with windows(3) and Windows(4) (I already have had an other graphic on screen opened in device 2). Now device 3 is inactive and device 4 is active. How i can make device 3 active to plot in device 3 and after make device 4 active to hav

[R] Res: Leyenda

(just tranlationg...) Hi, I am trying to put the legend on the outside of plot region. How can I do that? My code looks like... df<-read.table(stdin(),head=T,sep=",") row,LENGTH,LAT 091639,10.002,42.26282 091640,30.808,42.26834 091641,21.591,42.31689 091642,22.030,41.53246 091643,22.744,42.01954

Re: [R] Regex magic anyone?

strapply in the gsubfn package can split strings based on content rather than delimiters. > library(gsubfn) > strapply(s, ".[(].[)]|.")[[1]] [1] "A""B""C""D""E""F""G""T(P)" "H""I" [11] "J""K" See home page at: http://gsubfn.googlecode.com and vignette: library

[R] Regex magic anyone?

Hi again, how to elegantly split s <- "ABCDEFGT(P)HIJK" into "A" "B" "C" "D" "E" "F" "G" "T(P)" "H" "I" "J" "K" (independently of which letters 'T' or 'P' actually represent ...). Please jumstart my regexing, Joh __ R-help@r-project.org mailing li

Re: [R] Select a group of data from a file

Hi ? subset should do it xx < subset(x, x$Var2==1) --- "LA SPINA, MICHELANGELO" <[EMAIL PROTECTED]> wrote: > Hello everybody! > > I've a file with several data six variables, three > quantitative and three qualitative, I would like to > select a group of data from the file to analyze > then, i.

Re: [R] Leyenda

> Hola estoy intentando colocar la leyenda de un gráfico fuera de la plot > region, es posible hacer esto? Como? Pilar, En la linea de comandos yo escribí: RSiteSearch(""legend plot region") y el primer link es http://finzi.psych.upenn.edu/R/Rhelp02a/archive/43964.html de Uwe Ligges, y hay varios

[R] ApnaCircle Invite for Online Social and Professional Networking Site

[1]ApnaCircle.com Online Social & Professional Networking Site [2]Make Freinds Earn Rewards [3]ApnaCircle.com gives you a valuable platform to make friends & professional contacts for

Re: [R] Any tools for working with US 2000 census data? (Zembower, Kevin)

Hi Kevin, Not an R-specific solution, but by far the easiest way to extract data from the 2000 census is to use American FactFinder. Go to the main page: http://factfinder.census.gov/home/saff/main.html?_lang=en Hover over Data Sets, click Decennial Census. Select the file you want (SF-2 for yo

[R] image/area plot

Dear R users, I am trying to produce an image plot, that represents the proportions of a factor (z variable), so that the number of squares of each colour represents each factor level, with the respective label inside (sorry for the crap English). # Something like this: kk=data.frame(fact=let

[R] Leyenda

Hola estoy intentando colocar la leyenda de un gráfico fuera de la plot region, es posible hacer esto? Como? Mi código es el siguiente: > df LENGTH LAT 091639 10.002 42.26282 091640 30.808 42.26834 091641 21.591 42.31689 091642 22.030 41.53246 091643 22.744 42.01954 091644 12.702 42.67

Re: [R] Problems with lattice plot

Dear Francesc, This is rather easy to do with the ggplot2 package. library(ggplot2) dataset <- expand.grid(x = 1:10, y = 1:50) dataset$z <- runif(nrow(dataset)) ggplot(data = dataset, aes(x = x, y = y, fill = z)) + geom_tile() + scale_fill_gradient(low="white", high="black") + coord_equal() Hav

Re: [R] Problems with lattice plot

> I have some problems when I try to plot a lattice. The lattice has 200 > columns (X value) and 5 rows (Y value), with a total of 1000 squares. I > want to plot this lattice with the squares having different colours > according to a variable value. > > > > The variable I want to plot is cove

[R] Problems with lattice plot

Hello useRs, I have some problems when I try to plot a lattice. The lattice has 200 columns (X value) and 5 rows (Y value), with a total of 1000 squares. I want to plot this lattice with the squares having different colours according to a variable value. The variable I want to plot is co

Re: [R] Trouble receiving messages from the mailing list

On 18-Jan-08 08:32:02, Martin Maechler wrote: >> "GC" == Gang Chen <[EMAIL PROTECTED]> >> on Thu, 17 Jan 2008 12:08:40 -0500 writes: > > GC> For some unknown reason I stopped receiving any messages > GC> from the R- help mailing list. See if this test gets > GC> through. >

[R] constrOptim with method SANN

Hi Everyone, I'm trying to minimize a function using constrOptim with the simulated annealing method SANN. If I understand constrOptim well, it basically passes most of its arguments to optim while somehow enforcing the constraints. My problem is, that since SANN does not need gradients, when us

[R] constrOptim with SANN

Hi Everyone, I'm trying to minimize a function using constrOptim with the simulated annealing method SANN. If I understand constrOptim well, it basically passes most of its arguments to optim while somehow enforcing the constraints. My problem is, that since SANN does not need gradients, when us

Re: [R] Select a group of data from a file

> I've a file with several data six variables, three quantitative and > three qualitative, I would like to select a group of data from the > file to analyze then, i.e: > my file is like that (but with 6 variables): > > Var1 Var2 > 2 1 > 5 1 > 8 1 > 7 2 > 3 2 > 8

Re: [R] Select a group of data from a file

Try subset. help(subset) Subset returns a subsets of vectors, matrices or data frames which meet conditions. There are some excellent examples in the help. Here is an example similar to your email. Create some data: Var1<-sample(1:10,20,replace=TRUE) Var2<-sample(1:10,20,replace=TRUE) the.o

Re: [R] Trouble receiving messages from the mailing list

> "GC" == Gang Chen <[EMAIL PROTECTED]> > on Thu, 17 Jan 2008 12:08:40 -0500 writes: GC> For some unknown reason I stopped receiving any messages GC> from the R- help mailing list. See if this test gets GC> through. It does; but you don't see it -- nor my answer -- so may

[R] Select a group of data from a file

Hello everybody! I've a file with several data six variables, three quantitative and three qualitative, I would like to select a group of data from the file to analyze then, i.e: my file is like that (but with 6 variables): Var1 Var2 2 1 5 1 8 1 7 2 3 2 8 2

Re: [R] Trouble receiving messages from the mailing list

> "GC" == Gang Chen <[EMAIL PROTECTED]> > on Thu, 17 Jan 2008 12:08:40 -0500 writes: GC> For some unknown reason I stopped receiving any messages GC> from the R- help mailing list. See if this test gets GC> through. It does; but you don't see it -- nor my answer -- so may

Re: [R] "non-plot" plotting

Marc, you are kind of crazy ... I asked for a discussion of the path to take to get there, not to have it done ... but hey, what am I complaining about? ;0) Thank you very much! Joh Marc Schwartz wrote: > Johannes Graumann wrote: >> I really do not know ho to else title this ... I want to draw