[R] recommended way to group function calls in Sweave

2012-04-25 Thread Liviu Andronic
Dear all When using Sweave, I'm always hitting the same bump: I want to group repetitive calls in a function, but I want both the results and the function calls in the printed output. Let me explain myself. Consider the following computation in an Sweave document: summary(iris[,1:2])

Re: [R] recommended way to group function calls in Sweave

2012-04-25 Thread Liviu Andronic
On Wed, Apr 25, 2012 at 3:41 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: I would use the last method, or if the calls were truly repetitive (i.e. always identical, not just the same pattern), use a named chunk. Labeled chunks are indeed what I was looking for [1]. As far as I understand,

Re: [R] Organizations where IT has approved the use of R software

2012-04-13 Thread Liviu Andronic
On Fri, Apr 13, 2012 at 8:56 AM, Michael Bibo michael_b...@health.qld.gov.au wrote: Have you seen this from the archives? http://www.mail-archive.com/r-help@r-project.org/msg45768.html In addition see [1], which highlights the likes of Google, Merck or Stanford University. Liviu [1]

[R] plyr: set '.progress' argument to default to text

2012-04-10 Thread Liviu Andronic
Dear all Is it possible to set globally the option .progress = text to all the apply functions in 'plyr'. For example, current default is daply(..., .progress = none). I would like to set it to daply(..., .progress = text), so as to avoid writing the argument every time I call such a function. I

Re: [R] plyr: set '.progress' argument to default to text

2012-04-10 Thread Liviu Andronic
On Tue, Apr 10, 2012 at 5:32 PM, R. Michael Weylandt michael.weyla...@gmail.com wrote: You might try the Defaults package. Thanks for the hint. Unfortunately library(Defaults) setDefaults('create_progress_bar', name = text) doesn't do the trick. But one can set defaults for each *ply()

Re: [R] R generated means are different from the boxplot!

2012-04-06 Thread Liviu Andronic
Hello On Fri, Apr 6, 2012 at 10:35 AM, Jhope jeanwaij...@gmail.com wrote: Hi R-listers, 1) I am having trouble understanding why the means I have calculated from Aeventexhumed (A, B, and C) are different from the means showing on the boxplot I generated (see attached).  I have added the

Re: [R] reclaiming lost memory in R

2012-04-06 Thread Liviu Andronic
On Fri, Apr 6, 2012 at 2:21 PM, Ramiro Barrantes ram...@precisionbioassay.com wrote: Please let me know if you have any other suggestions or clues. See this older post by Brian [1] and check ?Memory-limits. Otherwise, I remember someone suggesting that even if R releases the memory internally,

Re: [R] spaghetti plots in R

2012-04-04 Thread Liviu Andronic
On Wed, Apr 4, 2012 at 4:22 PM, uday uday_143...@hotmail.com wrote: I would like to plat some spaghetti plots from my data , ma data is as See: require(sos) findFn('spaghetti') Liviu follows ak[1:3,]          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7] [,8]      [,9]

Re: [R] spaghetti plots in R

2012-04-04 Thread Liviu Andronic
On Wed, Apr 4, 2012 at 5:04 PM, uday uday_143...@hotmail.com wrote: Hi Liviu , thanks for post , but I could not find findFn('spaghetti') , I can see the following functions in sos package Extract.findFn findFn After installing 'sos', use the 'findFn()' function. For example, run

Re: [R] How to use Latex code in R loop?

2012-04-01 Thread Liviu Andronic
On Sun, Apr 1, 2012 at 4:48 PM, Manish Gupta mandecent.gu...@gmail.com wrote: Hi, I am newbie in Latex and R. I am working on one report in which i need to read file and display content of file by formatting (adding color boxes and colorful text for each record). For this i need to use latex

[R] proper order of calls when estimating nested models

2012-03-24 Thread Liviu Andronic
Dear all How should one proceed when estimating nested models containing missing data. What I would like to do is to first estimate the model with the control variables only, and then estimate the model containing also the variables of interest. For example, summary(reg.a - lm(IMC ~ STYLE + SEXE

Re: [R] Recommendations regarding textbooks

2012-03-22 Thread Liviu Andronic
On Thu, Mar 22, 2012 at 1:17 PM, richard willey richard.wil...@gmail.com wrote: I was hoping to get some advice regarding teaching R in an academic environment. thinks regarding “An R Companion to Applied Regression”. CAR features an excellent introduction to R programming. Liviu

Re: [R] vif function using lm object

2011-12-23 Thread Liviu Andronic
On Wed, Dec 21, 2011 at 9:28 AM, arunkumar akpbond...@gmail.com wrote: Hi,   can anyone please explain why the vif should have more than 2 terms. *vif.lm(lmobj) : model contains fewer than 2 terms* why it is throwng error if it is one variable. The _VIF_ is a measure of

Re: [R] Multicollinearty in logistic regression models

2011-12-16 Thread Liviu Andronic
On Fri, Dec 16, 2011 at 1:47 PM, David Winsemius dwinsem...@comcast.net wrote: Harrell's rms package has a vif function that is intended for use with fits from his logistic regression model function, lrm. This uses the variance Also see vif() in 'car'. Liviu

Re: [R] nice report generator?

2011-12-07 Thread Liviu Andronic
On Wed, Dec 7, 2011 at 7:43 PM, Janko Thyson janko.thyson.rst...@googlemail.com wrote: I like knitr. IMHO Yihui really came up with a killer package there: http://yihui.github.com/knitr/ If we're talking about nice, I'll chip in for LyX. It has support for Sweave, and will soon support knitr.

Re: [R] how to view/edit large matrix/array in R?

2011-12-05 Thread Liviu Andronic
On Tue, Dec 6, 2011 at 4:42 AM, Michael comtech@gmail.com wrote: I didn't mean the speed of saving the file is slow... I meant the manual procedures of exporting and then opening Excel, etc. is slow and inconvenient and unproductive... You may want to use RExcel then. Liviu On Mon, Dec

Re: [R] Help! Big problem when using browser() to do R debugging?

2011-12-03 Thread Liviu Andronic
On Sat, Dec 3, 2011 at 2:38 AM, Michael comtech@gmail.com wrote: Hi all, Could you please help me? I am having the following weird problem when debugging R programs using browser(): In my function, I've inserted a browser() in front of Step 1. My function has 3 steps and at the end of

Re: [R] detecting autocorrelation structure in panel data

2011-11-28 Thread Liviu Andronic
On Mon, Nov 28, 2011 at 3:54 AM, Matteo Richiardi matteo.richia...@gmail.com wrote: Hello, I'm a newby in R. I have created a data.frame holding panel data, with Take a look at the 'plm' package and its vignette. Regards Liviu the following columns: id,time,y, say: periods = 100 numcases

Re: [R] First read (was: Re: Looping and paste)

2011-11-24 Thread Liviu Andronic
On Thu, Nov 24, 2011 at 9:27 AM, Patrick Burns pbu...@pburns.seanet.com wrote: It's very seldom that I disagree with Bert, but here is one time. I don't think An Introduction to R is a suitable first read for people with little computational experience. I must agree with Patrick here. The

[R] understanding all.equal() output: Mean relative difference

2011-11-24 Thread Liviu Andronic
Dear all How should one parse all.equal() output? I'm specifically referring to the 'mean relative difference' messages. For example, all.equal(pi, 355/113) [1] Mean relative difference: 8.491368e-08 But I'm not sure how to understand these messages. When they're close to 0 (or 1xe-16), then

[R] proper work-flow with 'formula' objects and lm()

2011-11-24 Thread Liviu Andronic
Dear all I have a work-flow issue with lm(). When I use lm(y1~x1, anscombe) Call: lm(formula = y1 ~ x1, data = anscombe) Coefficients: (Intercept) x1 3.0001 0.5001 I get as expected the formula, y1 ~ x1, in the print()ed results or summary(). However, if I pass through a

Re: [R] proper work-flow with 'formula' objects and lm()

2011-11-24 Thread Liviu Andronic
On Thu, Nov 24, 2011 at 9:14 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: It is retained.  terms(fit) will give it to you, if fit is an lm object. Thank you. The following works nicely. (form - formula(y1~x1)) y1 ~ x1 x - lm(form, anscombe) formula(terms(x)) y1 ~ x1 However, I was

Re: [R] proper work-flow with 'formula' objects and lm()

2011-11-24 Thread Liviu Andronic
On Thu, Nov 24, 2011 at 10:25 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: Yes.  That's a job for substitute (the second time today). form - formula(y1~x1) x - eval(substitute(lm(f, anscombe), list(f = form))) summary(x) Call: lm(formula = y1 ~ x1, data = anscombe) That's what I

Re: [R] proper work-flow with 'formula' objects and lm()

2011-11-24 Thread Liviu Andronic
)    Call$formula - as.formula(terms(formula))    eval(Call) } From limited testing this seems to do exactly what I wanted! I thank all those who came up with suggestions. Regards Liviu On Thu, 24 Nov 2011, Liviu Andronic wrote: On Thu, Nov 24, 2011 at 10:25 PM, Prof Brian Ripley rip

Re: [R] [OT] 1 vs 2-way anova technical question

2011-11-22 Thread Liviu Andronic
On Tue, Nov 22, 2011 at 2:09 PM, Giovanni Azua brave...@gmail.com wrote: Mr. Gunter did not read/understand my problem, and there were no useful tips but only ad hominem attacks. By your side-taking I suspect you are in the same party club if you want to defend him maybe you should start by

Re: [R] changelog for MASS?

2011-11-17 Thread Liviu Andronic
On Thu, Nov 17, 2011 at 7:33 AM, R. Michael Weylandt michael.weyla...@gmail.com wrote: Hmmm...sorry -- the only thing I can suggest is maybe striking some sort of deal that you change when it gets however many months out of date:  if you look here

Re: [R] Running R with browser without installing anything

2011-10-20 Thread Liviu Andronic
On Thu, Oct 20, 2011 at 7:29 PM, Bogaso Christofer bogaso.christo...@gmail.com wrote: Dear all, the company I work for has Matlab installed for statistical/mathematical calculations and really not ready to go with R (even installing exe file for R). Therefore I was wondering is it possible to

Re: [R] R and Forex

2011-10-11 Thread Liviu Andronic
On Wed, Oct 12, 2011 at 3:29 AM, Yves S. Garret yoursurrogate...@gmail.com wrote: Hi all,   I recently started learning about Forex and found this O'Reilly book in Barnes Nobles about R.  I bought it out of pure curiosity.  I like what I see.  However, I have a question.  Has anyone tried to

Re: [R] Understanding the workflow between sweave, R and Latex

2011-10-03 Thread Liviu Andronic
On Fri, Sep 30, 2011 at 2:41 PM, Duncan Murdoch murdoch.dun...@gmail.com wrote: As an aside, I don't recommend the workflow you describe:  it's very slow and cumbersome.  It's much better to tell your text editor how to run both Sweave and Latex in one command.  In the upcoming release of R

[R] extracting p-values in scientific notation

2011-10-03 Thread Liviu Andronic
Dear all How does print.htest display the p-value in scientific notation? (x - cor.test(iris[[1]], iris[[3]])) Pearson's product-moment correlation data: iris[[1]] and iris[[3]] t = 21.65, df = 148, p-value 2.2e-16 alternative hypothesis: true correlation is not equal to 0 95 percent

Re: [R] extracting p-values in scientific notation

2011-10-03 Thread Liviu Andronic
Thanks all for your pointers. The following does trick: base::format.pval(x$p.value) ##Hmisc also has such a function [1] 2e-16 On Mon, Oct 3, 2011 at 10:25 AM, Rolf Turner rolf.tur...@xtra.co.nz wrote: Isn't it true that 0 2.2e-16? Yes, but it doesn't mean that the p-value actually hits

Re: [R] Supporting R/Membership

2011-09-25 Thread Liviu Andronic
On Sat, Sep 24, 2011 at 5:43 PM, Joshua Wiley jwiley.ps...@gmail.com wrote: I printed the form to mail off, but anecdotally I had planned on this for well over a year and kept putting it off.  If this is something that useRs and the R Foundation would be interested in, I would be happy to help

Re: [R] Reading R Code aloud

2011-09-14 Thread Liviu Andronic
On Wed, Sep 14, 2011 at 4:50 AM, cwdillon cwdil...@gmail.com wrote: Well, I've bee reading the R Manual, Nathan Yau's book on Flowing Data, and a few others like this to my kids to get them to go to sleep. I started off doing it the *wrong* way, just to keep it boring but slipped into doing it

Re: [R] Solve your R problems

2011-09-12 Thread Liviu Andronic
On Mon, Sep 12, 2011 at 9:41 AM, Patrick Burns pbu...@pburns.seanet.com wrote: R-help is all about solving R problems. So here ya go: http://www.portfolioprobe.com/2011/09/12/solve-your-r-problems/ Sweet. :) May I suggest a font change: anything but the default CM should do the trick. For one

Re: [R] Latex + R + sweave

2011-09-12 Thread Liviu Andronic
On Mon, Sep 12, 2011 at 9:21 AM, Twaha Mlwilo uddessy2...@hotmail.com wrote: Hello all, Good day, I have problem on how to remove the source code from the pdf output.Here I mean this.  code in sweave Rnw files  =  x-c(1,2,3,4,5,6) x  mean(x) sd(x) @ then would like it appear as

Re: [R] problem with plm package

2011-09-03 Thread Liviu Andronic
(please reply to the list, too, as it increases your chances of getting a helpful answer) On Sat, Sep 3, 2011 at 2:51 AM, loth lorien lothlorie...@hotmail.com wrote: Hi Liviu, Thank you for your suggestions, but they didn't seem to work. I got R to estimate a fixed effects model for me:

Re: [R] problem with plm package

2011-09-01 Thread Liviu Andronic
On Thu, Sep 1, 2011 at 4:25 AM, Ash lothlorie...@hotmail.com wrote: Hi, I am trying to complete a very simple panel analysis on some bank data. Call: Formula-plm(RoE~RoA+CAR+Inc.Dep+Cash.TL+NPL.Loans, data=Banks1, model=random, index=c(Bank.I.D.,Year)) summary(Formula) I get the

Re: [R] !!!function to do the knn!!!

2011-09-01 Thread Liviu Andronic
I would nominate the following fortune: On Thu, Sep 1, 2011 at 12:16 PM, Petr PIKAL petr.pi...@precheza.cz wrote: Try to call 112 if you are in Europe. David Winsemius: Thank your for your entry in the Poorly Capitalized and Inadequately Searched Posting Contest. mark: help, help ,help!!!

Re: [R] Treat an Unquoted Character String as a Data Frame

2011-09-01 Thread Liviu Andronic
On Thu, Sep 1, 2011 at 5:48 PM, Jean V Adams jvad...@usgs.gov wrote: Try this data - eval(parse(text=paste(study, level, ., population, sep=))) ..or this: data - get(paste(study, level, ., population, sep=)) Liviu Jean - dbateman wrote on 08/31/2011 17:44:44: I have several

Re: [R] how to split a data frame by two variables

2011-09-01 Thread Liviu Andronic
On Thu, Sep 1, 2011 at 7:53 PM, Changbin Du changb...@gmail.com wrote: HI, Dear R community, I want to split a data frame by using two variables: let and g It's not clear what you want to do, but investigate the following: require(plyr) Loading required package: plyr ddply(x, .(let, g),

Re: [R] Grouping variables in a data frame

2011-08-27 Thread Liviu Andronic
On Sat, Aug 27, 2011 at 7:26 AM, Andra Isan andra_i...@yahoo.com wrote: Hi All, I have a data frame as follow: user_id time age location gender . and I learn a logistic regression to learn the weights (glm with family= (link = logit))), my response value is either zero or one. I

Re: [R] obtaining p-values for lm.ridge() coefficients (package 'MASS')

2011-08-24 Thread Liviu Andronic
Dear Michael Thank you for your pointers. On Tue, Aug 23, 2011 at 4:05 PM, Michael Friendly frien...@yorku.ca wrote: First, you should be using rms::ols, as Design is old. Good to know. I've always wondered why Design and rms, in many cases, were providing similar functions and (upon cursory

Re: [R] obtaining p-values for lm.ridge() coefficients (package 'MASS')

2011-08-24 Thread Liviu Andronic
The attachment seems to have been dropped, so I'm pasting the code below. Sorry for that Liviu On Wed, Aug 24, 2011 at 1:44 PM, Liviu Andronic landronim...@gmail.com wrote: Second, penalty in ols() is not the same as the ridge constant in lm.ridge, but rather a penalty on the log likelihood

Re: [R] debugging functions in R

2011-08-24 Thread Liviu Andronic
On Wed, Aug 24, 2011 at 4:20 PM, Eran Eidinger e...@taykey.com wrote: Hi, I am not sure if this is the right list to ask this question (though I did not find a more appropriate one). I've started using R a month ago, and small scripts work fine. However, when I start writing more complex

[R] obtaining p-values for lm.ridge() coefficients (package 'MASS')

2011-08-23 Thread Liviu Andronic
Dear all I'm familiarising myself with Ridge Regressions in R and the following is bugging me: How does one get p-values for the coefficients obtained from MASS::lm.ridge() output (for a given lambda)? Consider the example below (adapted from PRA [1]): require(MASS) data(longley) gr -

Re: [R] How to navigate (zoom, pan) in a plot/graph

2011-08-21 Thread Liviu Andronic
On Sun, Aug 21, 2011 at 3:04 PM, Eran Eidinger e...@taykey.com wrote: Thank you Uwe, That solves the second question. Still looking for some solution to zooming in and out dynamically. Try package 'playwith'. Liviu __ R-help@r-project.org mailing

[R] Hmisc::rcorr on a 'data.frame'?

2011-08-19 Thread Liviu Andronic
Dear all ?Hmisc::rcorr states that it takes as main argument a numeric matrix. But is it normal that it fails in such an ugly way on a data frame? (See below.) If the function didn't attempt any conversion to a matrix, I would have expected it to state that in the error message that it didn't

[R] display only the top-right half of a correlation matrix?

2011-08-19 Thread Liviu Andronic
Dear all Is there an easy way to display only one half (top-right or bottom-left) of a correlation matrix? require(Hmisc) rcorr(as.matrix(mtcars[ , 1:4])) mpg cyl disphp mpg 1.00 -0.85 -0.85 -0.78 cyl -0.85 1.00 0.90 0.83 disp -0.85 0.90 1.00 0.79 hp -0.78 0.83 0.79

Re: [R] display only the top-right half of a correlation matrix?

2011-08-19 Thread Liviu Andronic
On Fri, Aug 19, 2011 at 9:02 PM, Peter Langfelder peter.langfel...@gmail.com wrote: Use as.dist: here's an example. Seems promising, but for one issue: I would like to keep the diagonal and thus specify 'diag=T', but then as.dist() replaces the diagonal values with zero. (See below.) Is there a

Re: [R] display only the top-right half of a correlation matrix?

2011-08-19 Thread Liviu Andronic
On Fri, Aug 19, 2011 at 9:38 PM, Peter Langfelder peter.langfel...@gmail.com wrote: if as.dist doesn't work, use brute force: x = matrix(rnorm(5*100), 100, 5) mat = signif(cor(x), 2); mat[lower.tri(mat)] = data.frame(mat) Yes, brute force works. This isn't quite how I wanted to do this,

Re: [R] round() a data frame containing 'character' variables?

2011-08-11 Thread Liviu Andronic
On Wed, Aug 10, 2011 at 5:52 PM, Jean V Adams jvad...@usgs.gov wrote: As it says in ?format the digits argument specifies ... how many significant digits are to be used ... enough decimal places will be used so that the smallest (in magnitude) number has this many significant digits ...

Re: [R] round() a data frame containing 'character' variables?

2011-08-11 Thread Liviu Andronic
Hello On Thu, Aug 11, 2011 at 4:34 PM, Martin Maechler maech...@stat.math.ethz.ch wrote: Well, as you say it's  View() only, so why don't you just view the result of something like          print(data, digits = 2) Unfortunately print(..., digits = 2) suffers of the same deficiencies as

Re: [R] round() a data frame containing 'character' variables?

2011-08-11 Thread Liviu Andronic
On Thu, Aug 11, 2011 at 4:58 PM, David Winsemius dwinsem...@comcast.net wrote: numVars - sapply(iris, is.numeric) iris[numVars] - lapply(iris[numVars], round, digits = 2) head(iris) Then drop the assignment operation and just print: ( lapply(iris[numVars], round, digits = 2) ) That's a

Re: [R] round() a data frame containing 'character' variables?

2011-08-11 Thread Liviu Andronic
On Wed, Aug 10, 2011 at 5:52 PM, Jean V Adams jvad...@usgs.gov wrote: data.frame(lapply(x, function(y) if(is.numeric(y)) round(y, 2) else y)) Building on this, I came up with a 'data.frame' method for round(). I'm not sure how orthodox it is, but it seems to do the job as intended:

Re: [R] convert 'list' to 'vector'?

2011-08-11 Thread Liviu Andronic
On Wed, Aug 10, 2011 at 11:39 PM, David Winsemius dwinsem...@comcast.net wrote: is.vector(x) [1] TRUE You probably didn't realize that lists _are_ vectors before this. Nice to know. Hopefully this expected invitation to grilling will not result in severe burns. I note that one of the

[R] round() a data frame containing 'character' variables?

2011-08-10 Thread Liviu Andronic
Dear all It is difficult to use round(..., digits=2) on a data frame since one has to first take care to remove non-numeric variables such as 'character' or 'factor': head(round(iris, 2)) Error in Math.data.frame(list(Sepal.Length = c(5.1, 4.9, 4.7, 4.6, 5, : non-numeric variable in data

Re: [R] round() a data frame containing 'character' variables?

2011-08-10 Thread Liviu Andronic
Hello On Wed, Aug 10, 2011 at 11:41 AM, Dimitris Rizopoulos d.rizopou...@erasmusmc.nl wrote: One approach is the following: numVars - sapply(iris, is.numeric) iris[numVars] - lapply(iris[numVars], round, digits = 2) head(iris) That's interesting, but still doesn't do what I need. Since it's

Re: [R] round() a data frame containing 'character' variables?

2011-08-10 Thread Liviu Andronic
Hello On Wed, Aug 10, 2011 at 2:31 PM, Jean V Adams jvad...@usgs.gov wrote: The function format() might serve your needs. This looks very promising, but yields some strange results. See below: x - data.frame(a=rnorm(10), b=rnorm(10), c=rnorm(10), d=letters[1:10]) x a b

[R] convert 'list' to 'vector'?

2011-08-10 Thread Liviu Andronic
Dear all How does one convert a non-symmetric list to a vector? See below: x - list() x[[1]] - letters[1:5] x[[2]] - letters[6:10] x[[3]] - letters[11:12] x [[1]] [1] a b c d e [[2]] [1] f g h i j [[3]] [1] k l paste(x) [1] c(\a\, \b\, \c\, \d\, \e\) c(\f\, \g\, \h\, \i\, \j\) [3] c(\k\,

Re: [R] convert 'list' to 'vector'?

2011-08-10 Thread Liviu Andronic
On Wed, Aug 10, 2011 at 9:02 PM, R. Michael Weylandt michael.weyla...@gmail.com wrote: unlist() Thanks all! This is perfect, and very R-ish: never where a novice would expect it to be. Cheers Liviu __ R-help@r-project.org mailing list

Re: [R] convert 'list' to 'vector'?

2011-08-10 Thread Liviu Andronic
On Wed, Aug 10, 2011 at 9:32 PM, David Winsemius dwinsem...@comcast.net wrote: Thanks all! This is perfect, and very R-ish: never where a novice would expect it to be. Well, since `unlist` is linked in the See Also on the help page for `list`, I can only hope you meant that in complete jest.

Re: [R] Very silent R

2011-08-07 Thread Liviu Andronic
On Fri, Aug 5, 2011 at 3:54 PM, khadeeja ismail haj...@yahoo.com wrote: other attached packages: [1] Rcmdr_1.6-2 I bet that this is what generates the lack of error messages. Check Rcmdr documentation on how to turn it off. Liviu __

Re: [R] Life Cycle Assessment with R.

2011-07-26 Thread Liviu Andronic
On Mon, Jul 25, 2011 at 8:03 PM, Jose Bustos Melo jbustosm...@yahoo.es wrote: Who knows if there's people working in Life Cycle Assesment (carbon emition) with R? or If there's someone interested in doing a package about it, please let me know! Some of the places to check would be: rseek.org

[R] add label attribute to objects?

2011-07-21 Thread Liviu Andronic
Dear all I know that the R way of documenting things is to work on your project in package development mode, and document each object (such as data frames) in a *.Rd files. This should work for gurus. What about a simpler way to document things, geared for mere mortals? I was thinking of a

Re: [R] Strange graphs

2011-07-18 Thread Liviu Andronic
On Mon, Jul 18, 2011 at 5:04 PM, UnitRoot akhussa...@gmail.com wrote: Hello, I am facing problems with plotting graphs in R. It was working well since today, but I cannot get some graphs for unknown reasons. Initially it gave me diagonal straight lines. Now it is giving me very strange curves.

Re: [R] Low Pain Unicode Characters in pdf graph?

2011-07-04 Thread Liviu Andronic
On Sun, May 15, 2011 at 3:06 PM, ivo welch ivo.we...@gmail.com wrote: Dear R-experts---is there a relatively low-pain way to get unicode characters into a plot to a pdf device? Have you tried Cairo package or cairo_pdf()? Both are making use of Cairo, which uses UTF-8 and automatically embeds

Re: [R] GUI Experience - Box Plots

2011-06-03 Thread Liviu Andronic
On Fri, Jun 3, 2011 at 6:46 AM, Vikas Garud information4vi...@gmail.com wrote: Hi, This is to share my experience of using GUI's for using various Quality Management/Assurance tools. A few days back, I had I had written a mail about Box plots and received some very good suggestions,

Re: [R] Box Plot under GUI (R Commander/RKward)

2011-05-17 Thread Liviu Andronic
On Tue, May 17, 2011 at 2:15 PM, Jannis bt_jan...@yahoo.de wrote:  Otherwise I would stick to SPSS/Excel/Sigmaplot or browse for some R GUIs with some basics beeing implemented on a click basis like RStudio, Rkward or TinnR. All R GUI implementations are however limited to basic

Re: [R] Powerful PC to run R

2011-05-16 Thread Liviu Andronic
On Fri, May 13, 2011 at 5:25 PM, John C Frain fra...@gmail.com wrote: Given the price of the Thinkpad X201 I would be very upset if it overheated and would contact my supplier for a replacement. Consider getting an aluminium cooling pad, such as those provided by Zalman. It helps keeping my HP

Re: [R] What does the - operator mean?

2011-04-21 Thread Liviu Andronic
On Thu, Apr 21, 2011 at 9:59 PM, Thomas Lumley tlum...@uw.edu wrote: The Evil and Wrong use is to modify variables in the global environment. I'm a bit at a loss. Why is it so wrong and evil? In other words, how should one modify variables in the global environment? Through the use of return()?

Re: [R] CairoPDF, Fonts, and Windows 7

2011-04-16 Thread Liviu Andronic
On Sat, Apr 16, 2011 at 6:39 PM, Al Roark hrbuil...@hotmail.com wrote: Hi All: I have some basic questions about Cairo graphics engine. I'm trying to use the Cairo package to produce PDF output, mainly because I perceive it to be easy to use with a wide variety of fonts. You may also want

[R] trouble with \SweaveOpts{grdevice=...}

2011-04-14 Thread Liviu Andronic
Dear all I've just tried the brand new 'grdevice' option in Sweave but couldn't make it work. When I declare \SweaveOpts{grdevice=pdf} or \SweaveOpts{grdevice=cairo_pdf} trying to plot something simple fig=T, echo=T= plot(1:10,1:10,main='Some title') @ would result in an Sweave error:

Re: [R] Sweave and Slides (Beamer)

2011-04-14 Thread Liviu Andronic
Hello I sympathise with you, since Beamer can quickly drive one mad. On Fri, Apr 15, 2011 at 3:11 AM, Brett Presnell presn...@stat.ufl.edu wrote: So, is this useful to anyone besides me?  What trick(s) am I missing that would make it easier/better, or that would obviate altogether the need

Re: [R] By function

2011-04-09 Thread Liviu Andronic
Hello On Fri, Apr 8, 2011 at 12:38 AM, Raoni Rodrigues caciquesamu...@gmail.com wrote: Hello all! I have a data frame with nine variables and 293 cases. (attached goes the csv file). The CSV didn't get through so it's difficult to replicate your example. Please post the output of: str(cpue)

Re: [R] R licence

2011-04-07 Thread Liviu Andronic
On Thu, Apr 7, 2011 at 9:45 AM, Stanislav Bek stanislav.pavel@gmail.com wrote: Hi, is it possible to use some statistic computing by R in proprietary software? I don't know if this covers your case, but SAS and SPSS provide interfaces to R. Regards Liviu Our software is written in c#,

[R] plyr workaround to converting by() to a data frame

2011-04-07 Thread Liviu Andronic
Dear all Is there a clean plyr version of the following by() and do.call(rbind, ...) construct: df-data.frame(a=1:10,b=11:20,c=21:30,grp1=c(x,y),grp2=c(x,y),grp3=c(x,y)) dfsum-by(df[c(a,b,c)], df[c(grp1,grp2,grp3)], range) as.data.frame(dfsum) Error in as.data.frame.default(dfsum) : cannot

[R] Sweave Cairo driver?

2011-04-06 Thread Liviu Andronic
Dear all I would like to use Sweave with the Cairo() graphics device instead of pdf(), since the former supports Unicode, allows for easier font selection out of a greater range of available fonts, and automatically embeds fonts into the resulting PDF. Following this older discussion [1], has

Re: [R] Sweave Cairo driver?

2011-04-06 Thread Liviu Andronic
On Thu, Apr 7, 2011 at 6:39 AM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: This is simple in 2.13.0 RC, which allows you to add graphics devices for Sweave. Interesting. Which help files should I look into? On Wed, 6 Apr 2011, Liviu Andronic wrote: since the former supports Unicode

Re: [R] Sweave Cairo driver?

2011-04-06 Thread Liviu Andronic
Dear Brian On Thu, Apr 7, 2011 at 6:39 AM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: allows for easier font selection out of a greater range of available fonts, As does cairo_pdf(), builtin to R and more reliable in my experience. I've looked into the help pages and searched the lists,

Re: [R] Syntax coloring in R console

2011-04-01 Thread Liviu Andronic
Hello On Fri, Apr 1, 2011 at 12:32 PM, January Weiner january.wei...@mpiib-berlin.mpg.de wrote: I have found the packages xterm256 and highlight, but I was not able to figure out how to use it to highlight the syntax in console output. I've been in this position before, without finding a

Re: [R] Syntax coloring in R console

2011-04-01 Thread Liviu Andronic
On Fri, Apr 1, 2011 at 3:48 PM, January Weiner january.wei...@mpiib-berlin.mpg.de wrote: RStudio might be a fine program, but it does not feature syntax highlighting, which is the only thing I am missing from R Console (it only colors the commands typed). The idea is that you shouldn't use the

Re: [R] Syntax coloring in R console

2011-04-01 Thread Liviu Andronic
hit More Set as working dir. A bit cumbersome, but seems to work, and a small price to pay for all the available functionality. Regards Liviu Thanks nonetheless, j. On Fri, Apr 1, 2011 at 4:20 PM, Liviu Andronic landronim...@gmail.com wrote: On Fri, Apr 1, 2011 at 3:48 PM, January Weiner

Re: [R] Performance Difference? Windows vs. Linux

2011-03-18 Thread Liviu Andronic
On Fri, Mar 18, 2011 at 11:00 PM, Brigid Mooney bkmoo...@gmail.com wrote: I'm not trying to start a Windows vs. Linux debate, but I've been using R on a Windows machine for a while, and was recently wondering if R's performance would be faster on a Linux machine.  And similarly, if any

Re: [R] Teaching R: To quote, or not to quote?

2011-03-09 Thread Liviu Andronic
On Wed, Mar 9, 2011 at 9:57 AM, Ingmar Visser i.vis...@uva.nl wrote: FWIW, the style for JSS now enforces (via Achimitization) the use of quotes as you show above. An excellent candidate for the fortunes package! We should probably address Achim for that. :) Regards Liviu best, Ingmar  

Re: [R] sorting subsetting a data.frame

2011-03-07 Thread Liviu Andronic
On Mon, Mar 7, 2011 at 1:38 AM, David Winsemius dwinsem...@comcast.net wrote: subset(x[order(x$Species1), ],  Sepal.Length==6.7 ) Thank you all for the suggestions. Now I can do exactly what I wanted. Regards Liviu __ R-help@r-project.org mailing list

[R] sorting subsetting a data.frame

2011-03-06 Thread Liviu Andronic
Dear all This may be obvious, but I cannot get it working. I'm trying to subset sort a data frame in one go. x - iris x$Species1 - as.character(x$Species) ##subsetting alone works fine with(x, x[Sepal.Length==6.7,]) ##sorting alone works fine with(x, x[order(Sepal.Length, rev(sort(Species1))),])

Re: [R] sorting subsetting a data.frame

2011-03-06 Thread Liviu Andronic
On Sun, Mar 6, 2011 at 11:49 PM, Liviu Andronic landronim...@gmail.com wrote: Dear all This may be obvious, but I cannot get it working. I'm trying to subset sort a data frame in one go. x - iris x$Species1 - as.character(x$Species) ##subsetting alone works fine with(x, x[Sepal.Length==6.7

Re: [R] sorting subsetting a data.frame

2011-03-06 Thread Liviu Andronic
On Sun, Mar 6, 2011 at 11:53 PM, Liviu Andronic landronim...@gmail.com wrote: On Sun, Mar 6, 2011 at 11:49 PM, Liviu Andronic landronim...@gmail.com wrote: Dear all This may be obvious, but I cannot get it working. I'm trying to subset sort a data frame in one go. x - iris x$Species1

Re: [R] problems with playwith

2011-03-01 Thread Liviu Andronic
On Tue, Mar 1, 2011 at 5:58 PM, R Heberto Ghezzo, Dr heberto.ghe...@mcgill.ca wrote: hello, i tried to run playwith but : library(playwith) Loading required package: lattice Loading required package: cairoDevice Loading required package: gWidgetsRGtk2 Loading required package: gWidgets

Re: [R] Is there any Command showing correlation of all variables in a dataset?

2011-03-01 Thread Liviu Andronic
On Tue, Mar 1, 2011 at 11:41 AM, JoonGi joo...@hanmail.net wrote: Thanks in advance. I want to derive correlations of variables in a dataset Specifically library(Ecdat) data(Housing) attach(Housing) cor(lotsize, bathrooms) this code results only the correlationship between two

Re: [R] When is *interactive* data visualization useful to use?

2011-02-19 Thread Liviu Andronic
On Fri, Feb 18, 2011 at 8:00 PM, Tom Hopper tomhop...@gmail.com wrote: Tal, One interactive capability that I have repeatedly wished for (but never taken the time to develop with the existing R tools) is the ability to interactively zoom in on and out of a data set, I believe that you can do

Re: [R] R-/Text-editor for Windows?

2011-01-29 Thread Liviu Andronic
Hello On Fri, Jan 28, 2011 at 8:03 PM, Jonathan P Daily jda...@usgs.gov wrote: Geany also handles a ton of filetypes, but lacks a direct interface to R in the windows version - I wrote an AutoHotKey script that did this for me in about 5 minutes. Could you please expand on that? I prefer to

Re: [R] how to get coefficient and scores of Principal component analysis in R?

2011-01-28 Thread Liviu Andronic
Hello On Fri, Jan 28, 2011 at 8:16 PM, Zunqiu Chen che...@ohsu.edu wrote: It might be a simple question. But I could not find the answer from function “prcomp” or “princomp”. Does anyone know what are the codes to get coefficient and scores of Principal component analysis in R? Try these:

Re: [R] speed up in R apply

2011-01-06 Thread Liviu Andronic
On Wed, Jan 5, 2011 at 10:49 PM, Young Cho young.s...@gmail.com wrote: When introduced to R, I learned how to use *apply whenever I could to avoid for-loops and all. And, getting the habit, I think I somehow got the mis-conception that it is a magic source, always an optimal way of coding in

Re: [R] Cost-benefit/value for money analysis

2011-01-05 Thread Liviu Andronic
On Wed, Jan 5, 2011 at 12:29 PM, Graham Smith myotis...@gmail.com wrote: maximal choices would break the budget. This sounds like a homework problem and I don't see any student effort yet. Search terms include: decision analysis , cost-benefit analysis, or utility theory. Hopefully,  my

Re: [R] Sweave for big data analysis

2011-01-01 Thread Liviu Andronic
(slightly off-topic) Hello On Sat, Jan 1, 2011 at 12:20 AM, Yihui Xie x...@yihui.name wrote: I still recommend the pgfSweave package (as usual) -- you can cache both data objects (using cacheSweave) and graphics (using pgf). Do these packages re-implement Sweave, or just use it as a

Re: [R] latex() etc.: How to nicely format a matrix for a LaTeX document?

2010-12-30 Thread Liviu Andronic
On Thu, Dec 30, 2010 at 9:48 PM, Marius Hofert m_hof...@web.de wrote: Thanks, John, I did. Here is an update [not all problems are solved; can they be?]: Consider: library(xtable) mat - matrix(c(1,NA,3,100,10012.23423,4), ncol = 3, byrow = TRUE) print(xtable(mat, digits = 1), floating =

Re: [R] What is the best way to lag a time series?

2010-12-26 Thread Liviu Andronic
On Sun, Dec 26, 2010 at 8:49 AM, Christian Schoder schoc...@newschool.edu wrote: Dear R-users, I've been using R for a while and I am very satisfied! Unfortunately, I still have not figured out an efficient and general way to construct and use lags of time series, especially when I need to

Re: [R] matrix indexing in 'for' loop?

2010-12-22 Thread Liviu Andronic
On Wed, Dec 22, 2010 at 2:57 AM, Phil Spector spec...@stat.berkeley.edu wrote: To make your loop work, you need to learn about the get function. I'm not going to give you the details because there are better approaches available. First, let's make some data that will give values which can be

Re: [R] matrix indexing in 'for' loop?

2010-12-22 Thread Liviu Andronic
On Wed, Dec 22, 2010 at 6:39 PM, govin...@msu.edu wrote: Thank you both for your suggestions. I have another question - is there a specific way to access the individual elements of a 'list' variable? i.e. dmi = matrix(rnorm(20),4,5) soi = matrix(rnorm(20),4,5) pe = matrix(rnorm(20),4,5) y

<    1   2   3   4   5   6   >