Re: [R] plotCI (plotrix) problem

2009-04-20 Thread Dieter Menne
Derek Ogle DOgle at northland.edu writes: I am attempting to create a plot with intervals stretched in the x-direction using plotCI() in the plotrix package. The same data provides an appropriate set of intervals when stretched in the y-direction but I only get a lower interval when

[R] to extract data

2009-04-20 Thread Roslina Zakaria
Hi R-users, I have a set of data from 1958-2009, how do I extract the data from 1927 and 2007?   beechworth.dt Year Month  Rain 1    1858 3  21.8 2    1858 4  47.0 3    1858 5  70.1 4    1858 6  78.7 5    1858 7 101.6 6    1858 8 129.8 7    1858 9  80.8 8   

Re: [R] to extract data

2009-04-20 Thread Uwe Ligges
Roslina Zakaria wrote: Hi R-users, I have a set of data from 1958-2009, how do I extract the data from 1927 and 2007? beechworth.dt Year Month Rain 11858 3 21.8 21858 4 47.0 31858 5 70.1 41858 6 78.7 51858 7 101.6 61858 8 129.8 7

Re: [R] to extract data

2009-04-20 Thread ronggui
If Year is numeric data, then you can use: beechworth.dt.2 - subset(beechworth.dt, Year=1997 Year =2008) If it is character, then you can use: beechworth.dt.2 - subset(beechworth.dt, Year %in% as.character(1997:2008)) 2009/4/20 Roslina Zakaria zrosl...@yahoo.com: Hi R-users, I have a set of

Re: [R] to extract data

2009-04-20 Thread djhurio
Hi, Try this one: beechworth.dt.2 - beechworth.dt[beechworth.dt$Year=1927 beechworth.dt$Year=2007,] Martins On Apr 20, 9:59 am, Roslina Zakaria zrosl...@yahoo.com wrote: Hi R-users, I have a set of data from 1958-2009, how do I extract the data from 1927 and 2007?   beechworth.dt   

Re: [R] to extract data

2009-04-20 Thread Christian Ritz
Hi, maybe the following line works for you: beechworth.dt.2 - beechworth.dt[as.numeric(beechworth.dt$Year) %in% 1927:2007, ] (using as.numeric() to make sure that Year is numeric, maybe not needed?). Christian __ R-help@r-project.org mailing list

[R] how can I run to multinomial PLS regression in R?

2009-04-20 Thread 풍이
__ 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 provide commented, minimal, self-contained, reproducible code.

Re: [R] to extract data

2009-04-20 Thread Stefan Grosse
On Sun, 19 Apr 2009 23:59:43 -0700 (PDT) Roslina Zakaria zrosl...@yahoo.com wrote: RZ I have a set of data from 1958-2009, how do I extract the data from RZ 1927 and 2007? RZ beechworth.dt RZ Year Month  Rain how about: beech.cut-subset(beechworth.dt,(Year1926Year2008)) hth Stefan

Re: [R] savePlot error when type = eps or wmf

2009-04-20 Thread Paul Murrell
Hi jimdare wrote: Hello, When I use savePlot(filename=xy,type=eps) or savePlot(filename=xy,type=wmf) , I get the following error: Error in grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : Polygon edge not found (zero-width or zero-height?) This doesn't occur when I

Re: [R] (senza oggetto)

2009-04-20 Thread ONKELINX, Thierry
Giuseppe, Please enter a meaningfull subject line. Have you tried the examples from the variogram helpfile? They clearly demonstrate how to do it. If you still have problems, then you need to give us a reproducible example of your dataset and code. HTH, Thierry

Re: [R] Buglet in plotCI

2009-04-20 Thread Jim Lemon
Dieter Menne wrote: Hi, Jim, there is a typo at the bottom of plotCI: there is an y.to.in which should be x.to.in. See list for an example. http://article.gmane.org/gmane.comp.lang.r.general/147103 Should be: nz - (abs(ui - pmin(x + gap, ui)) * x.to.in) 0.001 Dieter Hi Dieter,

[R] Odp: (senza oggetto)

2009-04-20 Thread Petr PIKAL
Hi I found that variogram is a function in spatial library. I tried ??variogram library(spatial) ?variogram I used an example and I did not have any problem with it. data(topo, package=MASS) topo.kr - surf.ls(2, topo) variogram(topo.kr, 25) I do not have gstat available but from

Re: [R] problem with new version

2009-04-20 Thread Jim Lemon
Roslina Zakaria wrote: Hi R-users, I just change to the new version of R. I just wonder why everytime I run my function I will get this message: source(.trPaths[4], echo=TRUE, max.deparse.length=1) Thank you. Hi Roslina, If your function consists of the line:

Re: [R] data$ID - I always get a NULL

2009-04-20 Thread Grześ
This is my result: class(data) [1] data.frame str(data) 'data.frame': 2193 obs. of 83 variables: $ X.ID. : Factor w/ 2193 levels '18201',..: 1 2 3 4 5 6 7 8 9 10 ... $ X.kod. : Factor w/ 20 levels '01','02',..: 1 1 1 1 1 1 1 1 1 1 ...

[R] problem with new version

2009-04-20 Thread Roslina Zakaria
Hi R-users, I just change to the new version of R.  I just wonder why everytime I run my function I will get this message: source(.trPaths[4], echo=TRUE, max.deparse.length=1) Thank you. __ R-help@r-project.org mailing list

Re: [R] Sweave: Changing the background color, adding a border

2009-04-20 Thread Dieter Menne
Christophe Genolini cgenolin at u-paris10.fr writes: ... Thanks for your answer. I finally succeed. I used the listings package to define an environment that put code in a grey box. The command is \lstnewenvironment{Sinput}[1][]{ ... This is too nice code to be lost, but it suffers from

[R] [R-pkgs] xterm256

2009-04-20 Thread Romain Francois
This is the first release of the xterm256 package. xterm256 is a small package that takes advantage of the 256 color mode of xterm, enabling use of foreground and background colors in the R console. The package exposes one function style that takes three arguments: - (x) the text to style -

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Uwe Ligges
jjh21 wrote: Hello, The journal I am publishing in requires MS Word files. What is my best option for getting a high quality image of a graph done in R into Word? JPEG? Postscript? Windows metafile, if you are under Windows anyway. Uwe Ligges Thanks.

[R] R graph into MS Word: which format to use?

2009-04-20 Thread jjh21
Hello, The journal I am publishing in requires MS Word files. What is my best option for getting a high quality image of a graph done in R into Word? JPEG? Postscript? Thanks. -- View this message in context:

[R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread justin bem
See ade4 or mva package.  Justin BEM BP 1917 Yaoundé De : nikolay12 nikola...@gmail.com À : r-help@r-project.org Envoyé le : Lundi, 20 Avril 2009, 4h37mn 41s Objet : [R] PCA and automatic determination of the number of components Hi all, I have relatively

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Simon Pickett
I actually get superior results from creating a pdf, opening it in adobe acrobat, adjust the image size so it is big as possible with the screen, then copying it into word (by using the little square capture icon). HTH, Simon. - Original Message - From: jjh21 jjhar...@gmail.com To:

Re: [R] PCA and automatic determination of the number of components

2009-04-20 Thread Bruno Falissard
You can also use parallel analysis using the scree.plot function of the psy package. Regards, Bruno Bruno Falissard INSERM U669, PSIGIAM Paris Sud Innovation Group in Adolescent Mental Health Maison de Solenn 97

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Stefan Grosse
On Mon, 20 Apr 2009 03:01:42 -0700 (PDT) jjh21 jjhar...@gmail.com wrote: J The journal I am publishing in requires MS Word files. What is my J best option for getting a high quality image of a graph done in R J into Word? JPEG? Postscript? I use png or eps. The latter has the better output but

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Gabor Grothendieck
Use wmf format. That's a vector format rather than bitmapped so it will give better quality. Also you will be able to edit the image right in Word, e.g. change axis labels. On Mon, Apr 20, 2009 at 6:01 AM, jjh21 jjhar...@gmail.com wrote: Hello, The journal I am publishing in requires MS

[R] R-Squared with biglm?

2009-04-20 Thread Bryan Lim
I've been working with a rather large data set (~10M rows), and while biglm works beautifully for generating coefficients, it does not report an r-squared. It does report RSS. Any idea on how one could coax an R-squared out of biglm? Thanks in advance for any help with this! Bryan Lim

[R] Two or more dimensional root (Zero) finding

2009-04-20 Thread enrico.fosco...@libero.it
Good morning to all, I should find the zero of a specific function with respect to a vector of arguments. Does it exist something similar in R? Thank you very much, Enrico Foscolo __ R-help@r-project.org mailing list

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Duncan Mackay
Word can be memory hungry when using graphics. Another alternative is the metafile which is very compact and can give good results eg from the example win.metafile(Rplot%02d.wmf, pointsize = 10) Regards Duncan Mackay Department of Agronomy Soil Science University of New England ARMIDALE NSW

Re: [R] data$ID - I always get a NULL

2009-04-20 Thread David Winsemius
On Apr 20, 2009, at 4:33 AM, Grześ wrote: This is my result: class(data) [1] data.frame str(data) 'data.frame': 2193 obs. of 83 variables: $ X.ID. : Factor w/ 2193 levels '18201',..: 1 2 3 4 5 6 7 8 9 10 ... $ X.kod. : Factor w/

Re: [R] Buglet in plotCI

2009-04-20 Thread Derek Ogle
Jim and Dieter, Thank you for the quick reply and precise fix. My plot now works as expected. Thanks again. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of Jim Lemon Sent: Monday, April 20, 2009 4:21 AM To:

Re: [R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread Jari Oksanen
justin bem justin_bem at yahoo.fr writes: See ade4 or mva package.  Justin BEM BP 1917 Yaoundé I guess the problem was not to find PCA (which is easy to find), but finding an automatic method of selecting (determining sounds like that selection would be correct in some objective sense)

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Frank E Harrell Jr
jjh21 wrote: Hello, The journal I am publishing in requires MS Word files. What is my best option for getting a high quality image of a graph done in R into Word? JPEG? Postscript? Thanks. Double check the instructions for authors. Most journals take pdf and eps as separate files. Frank

Re: [R] help with this code

2009-04-20 Thread Tal Galili
Did you try using:debug(f2) and then running the function, and see what is your variables condition before the error ? On Sun, Apr 19, 2009 at 12:58 PM, li li hannah@gmail.com wrote: Hi, can anyone help me with the following code? Thanks! library(mvtnorm) f2 - function(n, rho) { var -

Re: [R] Two or more dimensional root (Zero) finding

2009-04-20 Thread Ben Bolker
enrico.fosco...@libero.it wrote: Good morning to all, I should find the zero of a specific function with respect to a vector of arguments. Does it exist something similar in R? Someone else may pipe up with a better answer, but this is generally a difficult problem (see _Numerical

[R] system is exactly singular

2009-04-20 Thread onyourmark
Hi. I have a csv file. I imported it with mydata-read.table(C:/dataForR/radiology/WordFrequency.csv, header=TRUE, sep=,) dim(mydata) [1] 982 925 The first column had the doc numbers like doc1, doc2, etc. so I did mydataNum-mydata[,-1] dim(mydataNum) [1] 982 924 The second to last column

[R] Paired test for 2-way ANOVA

2009-04-20 Thread A Ezhil
Hi, I have an experimental data with 2 factors: visit and treatment. Each subject has 2 visits and on each visit they get a treatment. We take samples before and after treatment. I can easily to do the 2-way ANOVA in R with visit and treatment as factors. anova( lm(data ~ visit*treatment) )

[R] bug when subtracting decimals?

2009-04-20 Thread wolfgang.siewert
Try this: 0.7-0.3==0.4 (We get FALSE) 0.7-0.30.4 (We get TRUE) but 0.8-0.3==0.5 (TRUE) 0.8-0.30.5 (FALSE) Funny, he? There is a way around: round(0.7-0.3,1)==0.4 (TRUE) Obviously there is a problem with some combinations of decimal subtractions, that - we have the feeling - shouldt be

Re: [R] problem with new version

2009-04-20 Thread Ben Bolker
Roslina Zakaria wrote: Hi R-users, I just change to the new version of R.  I just wonder why everytime I run my function I will get this message: source(.trPaths[4], echo=TRUE, max.deparse.length=1) Are you by any chance using Tinn-R? Could be a glitch in the interaction

[R] doing zero inflated glmm for count data with fmr

2009-04-20 Thread levyofi
Hello R users, Doing My PhD I collected count data which I believe is zero inflated. I have run a statistical model with lmer and family=poisson and got summary(model)@sigma=1 so I believe there is no overdispertion. I would like to use the fmr function from the 'gnlm' library but I just cannot

Re: [R] bug when subtracting decimals?

2009-04-20 Thread Henrique Dallazuanna
This is a FAQ. Try this: all.equal(0.7-0.3, 0.4) On Mon, Apr 20, 2009 at 10:07 AM, wolfgang.siewert wolfgang.siew...@gmail.com wrote: Try this: 0.7-0.3==0.4 (We get FALSE) 0.7-0.30.4 (We get TRUE) but 0.8-0.3==0.5 (TRUE) 0.8-0.30.5 (FALSE) Funny, he? There is a way around:

Re: [R] bug when subtracting decimals?

2009-04-20 Thread Dimitris Rizopoulos
this is a (very) Frequently Asked Question; check: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f Best, Dimitris wolfgang.siewert wrote: Try this: 0.7-0.3==0.4 (We get FALSE) 0.7-0.30.4 (We get TRUE) but 0.8-0.3==0.5 (TRUE) 0.8-0.30.5

Re: [R] Two or more dimensional root (Zero) finding

2009-04-20 Thread Ravi Varadhan
Ben, That is not a good approach, i.e. finding the zero, x*, of F(x), such that F(x*) = 0, as a minimum of ||F(x)|| is NOT a good approach. Any root of F(x) is indeed a global minimum of ||F(x)||, or for that matter, the global minimum of any f(F(x)), where f(.) is a mapping from R^p to R, such

Re: [R] Paired test for 2-way ANOVA

2009-04-20 Thread Mike Lawrence
Paired (aka. Repeated measures, aka. within-Ss) tests can be achieved by using aov() and specifying the within-Ss effect in the error term: my_aov = aov( dependent_variable~between_Ss_variable*within_Ss_variable + Error(Ss_id/(within_Ss_variable))) summary(my_aov) On Mon, Apr 20, 2009 at 10:25

[R] how to access file backed big matrix (package bigmemory)

2009-04-20 Thread utkarshsinghal
__ 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 provide commented, minimal, self-contained, reproducible code.

Re: [R] explicit documentation

2009-04-20 Thread Duncan Murdoch
ronggui wrote: It is always unfair to complain about volunteer work, and what you should do is to make contributions. I only half agree with this: I think that it's fair to complain, as long as you make contributions. With documentation, if you don't think it's clear, as part of your

[R] How estimate VAR(p)-model robustly?

2009-04-20 Thread Irene Schreiber
Hello, Does anyone know about robust estimation of vector autoregressive models (VAR(p)) in R? Or in Matlab? Currently I am using the function ar(). The problem is, that the variances of my data change a lot with time, and we also have some outliers in the data. That is why, I presume, that

[R] setting levels in contourplot()

2009-04-20 Thread Cable, Samuel B Civ USAF AFMC AFRL/RVBXI
Does the contourplot() routine have an argument analogous to the levels argument in the contour() routine? More generally, is there a way for the user to fix the contour levels in contourplot()? Thanks. --Sam Cable __ R-help@r-project.org mailing

Re: [R] Two or more dimensional root (Zero) finding

2009-04-20 Thread Ben Bolker
Ravi Varadhan wrote: That is not a good approach, i.e. finding the zero, x*, of F(x), such that F(x*) = 0, as a minimum of ||F(x)|| is NOT a good approach. Any root of F(x) is indeed a global minimum of ||F(x)||, or for that matter, the global minimum of any f(F(x)), where f(.) is a mapping

Re: [R] setting levels in contourplot()

2009-04-20 Thread David Winsemius
On Apr 20, 2009, at 10:30 AM, Cable, Samuel B Civ USAF AFMC AFRL/RVBXI wrote: Does the contourplot() routine have an argument analogous to the levels argument in the contour() routine? More generally, is there a way for the user to fix the contour levels in contourplot()? Thanks.

Re: [R] doing zero inflated glmm for count data with fmr

2009-04-20 Thread Ben Bolker
levyofi wrote: Hello R users, Doing My PhD I collected count data which I believe is zero inflated. I have run a statistical model with lmer and family=poisson and got summary(model)@sigma=1 so I believe there is no overdispertion. You have been misled. sigma is set to 1 by

Re: [R] setting levels in contourplot()

2009-04-20 Thread Dieter Menne
Cable, Samuel B Samuel.Cable at hanscom.af.mil writes: Does the contourplot() routine have an argument analogous to the levels argument in the contour() routine? More generally, is there a way for the user to fix the contour levels in contourplot()? Thanks. Assuming lattice: cuts Dieter

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Dieter Menne
Uwe Ligges ligges at statistik.tu-dortmund.de writes: The journal I am publishing in requires MS Word files. What is my best option for getting a high quality image of a graph done in R into Word? JPEG? Postscript? Windows metafile, if you are under Windows anyway. But make sure that

Re: [R] setting levels in contourplot()

2009-04-20 Thread David Winsemius
On Apr 20, 2009, at 11:22 AM, Dieter Menne wrote: Cable, Samuel B Samuel.Cable at hanscom.af.mil writes: Does the contourplot() routine have an argument analogous to the levels argument in the contour() routine? More generally, is there a way for the user to fix the contour levels in

[R] importing spreadsheet data - linera regression - panel data

2009-04-20 Thread Millo Giovanni
Dear Cecilia, just adding some examples to Stefan's post, which says everything already. I've recently gone mad with reshaping, so I assume it is a little tricky. Or maybe what I tell you is obvious, then just skip it. **import** Your files are spreadsheets, so the best way to import is to save

[R] what is R best for; what should one learn in addition to R

2009-04-20 Thread Juliet Hannah
Hi, I've been working with R for a couple of years, and I've been able to get most of the things done that I needed (sometimes in a roundabout way). A few experienced statisticians told me that R is best for interactive data analysis, but for large-scale computations, one needs something else. I

[R] Constraining equality between parameters of multinomial logit regression using VGAM)

2009-04-20 Thread Rémi Kazma
Hello, I am using the VGAM library to perform a multinomial logistic regression with 3 outcome categories and a binary independant factor. How can I contraint the beta [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Constraining equality between parameters of multinomial logit regression using VGAM)

2009-04-20 Thread Rémi Kazma
Hello, I am using the VGAM library to perform a multinomial logistic regression with 3 outcome categories and a binary independant factor. How can I constraint the 2 beta parameters to be equal. I intend to do that in order to obtain a Log-Likelihood to test against the model with 2 different beta

Re: [R] data$ID - I always get a NULL

2009-04-20 Thread Grześ
Now everything is clear for me. Thanks David! :) David Winsemius wrote: On Apr 20, 2009, at 4:33 AM, Grześ wrote: This is my result: class(data) [1] data.frame str(data) 'data.frame':2193 obs. of 83 variables: $ X.ID. : Factor w/ 2193

[R] How to force axis to have the same range

2009-04-20 Thread Sebastien Bihorel
Dear R-users, I am trying to produce (standard and trellis) scatterplots which use the same range of the x and y axes. In addition, I would like the plots to be physically square. Is there one or more specific argument(s) to plot and xyplot what would do that? I have looked at various

Re: [R] How to force axis to have the same range

2009-04-20 Thread Dieter Menne
Sebastien Bihorel Sebastien.Bihorel at cognigencorp.com writes: I am trying to produce (standard and trellis) scatterplots which use the same range of the x and y axes. In addition, I would like the plots to be physically square. Is there one or more specific argument(s) to plot and xyplot

Re: [R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread William Revelle
At 12:08 PM + 4/20/09, Jari Oksanen wrote: justin bem justin_bem at yahoo.fr writes: See ade4 or mva package. Justin BEM BP 1917 Yaoundé I guess the problem was not to find PCA (which is easy to find), but finding an automatic method of selecting (determining sounds like that

Re: [R] doing zero inflated glmm for count data with fmr

2009-04-20 Thread levyofi
Thank you very much Ben. I took you advice and run the lmer model with family=quasipoisson and got sigma~10. I guess this means that the data has overdispertion but not too high (15) for me to must use zero inflated model. Am I right? I will also post to the r-sig-mixed-models specialty list...

Re: [R] How to force axis to have the same range

2009-04-20 Thread David Winsemius
On Apr 20, 2009, at 12:18 PM, Dieter Menne wrote: Sebastien Bihorel Sebastien.Bihorel at cognigencorp.com writes: I am trying to produce (standard and trellis) scatterplots which use the same range of the x and y axes. In addition, I would like the plots to be physically square. Is there

Re: [R] bug when subtracting decimals?

2009-04-20 Thread Dieter Menne
wolfgang.siewert wolfgang.siewert at gmail.com writes: There is a way around: round(0.7-0.3,1)==0.4 (TRUE) Obviously there is a problem with some combinations of decimal subtractions, that - we have the feeling - shouldt be solved. Oh no, not that one again! This was lecture two in my

Re: [R] drawing data structures

2009-04-20 Thread Greg Snow
I am not sure that I fully understand what you want, but look at the TkListView function in the TeachingDemos package for one possible solution. If that is not what you meant, then a better problem description will help us to help you. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center

Re: [R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread Stéphane Dray
ade4 has the 'testdim' function which implements a recent method for estimating the number of dimension for PCA on correlation matrix. Paper describing the approach is available at http://pbil.univ-lyon1.fr/members/dray/articles/SD805.php William Revelle wrote: At 12:08 PM + 4/20/09,

[R] automatic exploration of all possible loglinear models?

2009-04-20 Thread Christopher W. Ryan
Is there a way to automate fitting and assessing loglinear models for several nominal variables . . . something akin to step or drop1 or add1 for linear or logistic regression? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch

Re: [R] how can I run to multinomial PLS regression in R?

2009-04-20 Thread Max Kuhn
There are functions in the caret package to do classification with pls and spls, if that is what you mean. They are plsda and splsda. -- Max __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] what is R best for; what should one learn in addition to R

2009-04-20 Thread Stefan Grosse
On Mon, 20 Apr 2009 11:49:29 -0400 Juliet Hannah juliet.han...@gmail.com wrote: JH I've been working with R for a couple of years, and I've JH been able to get most of the things done that I needed (sometimes in JH a roundabout way). A few experienced statisticians told me that JH R is best for

Re: [R] Random Forests: Question about R^2

2009-04-20 Thread Dimitri Liakhovitski
I would like to summarize. Would you please confirm that my summary is correct? Thank you very much! Determining R^2 in Random Forests (for a Regression Forest): 1. For each individual case, record a mean prediction on the dependent variable y across all trees for which the case is OOB

[R] Sweave and executive summaries - WORK AROUND

2009-04-20 Thread Karen_Byron
__ 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 provide commented, minimal, self-contained, reproducible code.

Re: [R] Re : PCA and automatic determination of the number of components

2009-04-20 Thread nikolay12
Thanks to all for the suggestions. Are you aware of a convenient implementation of AIC/BIC for the problem of selecting the number of factors? Nick William Revelle wrote: At 12:08 PM + 4/20/09, Jari Oksanen wrote: justin bem justin_bem at yahoo.fr writes: See ade4 or mva

[R] generic genotype calling algorithm?

2009-04-20 Thread Nigel Birney
Hi, I have Agilent and Illumina SNP data. That's the only thing I know about the files - there seem to be no version specification in any of them. Is there a generic genotype calling algorithm that I could try to use on these datasets? thanks, N. -- View this message in context:

[R] SEM package

2009-04-20 Thread Tijana Gonja
Hi I tried to install the sem package and it was succsessful but I keep getting the following message Paket 'sem' erfolgreich ausgepackt und MD5 Summen abgeglichen Die heruntergeladenen Pakete sind in C:\Users\Dean\AppData\Local\Temp\RtmpVRxlLZ\downloaded_packages aktualisiere HTML

[R] RES: How to force axis to have the same range

2009-04-20 Thread Rodrigo Aluizio
Hi Sebastien, take a look at the par(pty='s') argument. Maybe its can solve your issue. Best wishes. Rodrigo. -Mensagem original- De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em nome de Sebastien Bihorel Enviada em: segunda-feira, 20 de abril de 2009 12:26

Re: [R] what is R best for; what should one learn in addition to R

2009-04-20 Thread Josh Stumpf
You might also try looking in to Python. It has lots of mathematical and computational libraries (numpy and scipy are particularly useful), as well some to deal with concurrency (e.g. Parallel Python), and can interact with R (rpy and rpy2). Josh Stumpf On Mon, Apr 20, 2009 at 8:49 AM, Juliet

Re: [R] what is R best for; what should one learn in addition to R

2009-04-20 Thread David M Smith
On Mon, Apr 20, 2009 at 10:37 AM, Stefan Grosse singularit...@gmx.net wrote: I know of people doing optimization stuff which needs a lot of computational power. They use Matlab since it is easy for them to use multiple processors (+multiple pc's). R at the moment only uses one processor and

Re: [R] SEM package

2009-04-20 Thread Dimitri Liakhovitski
Nina, did you try to install this package using some other mirror? Dimitri On Mon, Apr 20, 2009 at 1:35 PM, Tijana Gonja tijana.go...@gmail.com wrote: Hi I tried to install the sem package and it was succsessful but I keep getting the following message Paket 'sem' erfolgreich ausgepackt und

Re: [R] R graph into MS Word: which format to use?

2009-04-20 Thread Mark Wardle
They may stipulate word files for the manuscript, but during submission, journals usually request EPS or PDF formats. In fact, all I've dealt with stipulate NEVER include graphcis into a Word document. May be worth checking again! Mark 2009/4/20 jjh21 jjhar...@gmail.com: Hello, The journal

[R] Random Forests: Predictor importance for Regression Trees

2009-04-20 Thread Dimitri Liakhovitski
Hello! I think I am relatively clear on how predictor importance (the first one) is calculated by Random Forests for a Classification tree: Importance of predictor P1 when the response variable is categorical: 1. For out-of-bag (oob) cases, randomly permute their values on predictor P1 and then

[R] R Golf?

2009-04-20 Thread Lauri Nikkinen
Hello, This is probably off-topic but have the R community ever organized R Golf contents similar to Perl Golf: http://perlgolf.sourceforge.net/ ? It would be nice to see how R gurus solve problems in many different ways (although you see it at this list every day :-)). Regards, -L

Re: [R] Calling objects in a loop

2009-04-20 Thread Duncan Murdoch
On 4/20/2009 2:35 PM, Brendan Morse wrote: Hi everyone, I am trying to calculate a particular variable (vector) from some previously defined variables in a loop but I am having trouble figuring out how to get the loop to recognize that it should index for the previously defined objects. Here is

Re: [R] Calling objects in a loop

2009-04-20 Thread Dimitri Liakhovitski
Brendan, I think you should create objects outside of the for loop. You can't create objects instide the loop. You can try this: metalist1-list() for(i in 1:10) {metalist1[[i]]-assign(paste(theta1_,i,sep=),data.frame(scale(rnorm(250} lapply(metalist1,function(x){print(dim(x))}) # Checking

[R] graph with 15 combinations

2009-04-20 Thread Penner, Johannes
Dear R helpers, I have a data set with 4 types (W, C, E S). Now I have values for all types plus all possible combinations (the order is unimportant): W, C, WC, E, WE, CE, WCE, S, WS, CS, WCS, ES, WES, CES WCES. Ideally I would like to represent everything in one graph and as concise as

[R] Matrix package,solve() errors and crashes

2009-04-20 Thread Surendar Swaminathan
Hello All, I am working on graph object using IGRAPH package wanted to do Bonacich Power. This is my graph object. The file 'Graph.RData' (4.2 MB) is available for download at http://dropbox.unl.edu/uploads/20090424/cfe4fcb854bb17f2/Graph.RData Graph size Vertices: 20984 Edges: 326033

Re: [R] Calling objects in a loop

2009-04-20 Thread Jorge Ivan Velez
Hi Brendan, If you really, really want to work with the for() loop, then for(i in 1:10){ assign(paste(theta1_,i,sep=),data.frame(scale(rnorm(250 assign(paste(theta2_,i,sep=),data.frame(scale(rnorm(250 assign(paste(theta3_,i,sep=),

[R] Calling objects in a loop

2009-04-20 Thread Brendan Morse
Hi everyone, I am trying to calculate a particular variable (vector) from some previously defined variables in a loop but I am having trouble figuring out how to get the loop to recognize that it should index for the previously defined objects. Here is a simplified version of what I am trying to

Re: [R] bug when subtracting decimals?

2009-04-20 Thread Gavin Simpson
Dieter Menne wrote: wolfgang.siewert wolfgang.siewert at gmail.com writes: There is a way around: round(0.7-0.3,1)==0.4 (TRUE) Obviously there is a problem with some combinations of decimal subtractions, that - we have the feeling - shouldt be solved. Oh no, not that one again! This was

Re: [R] what is R best for; what should one learn in addition to R

2009-04-20 Thread Greg Snow
What things you should learn in addition to R depends on what types of things you want to do. But here are some of the programs that I recommend people take a look at: gnuplot ggobi Imagemagick LaTeX xfig Make Emacs (or other programming editor) SQL Perl (or other scripting language) Tk

Re: [R] SEM package

2009-04-20 Thread John Fox
Dear Tijana, I'm afraid that I can't entirely read the messages, but the problem doesn't look like it's peculiar to the sem package. It appears as if the package may have been properly installed but you didn't have permission to update HTML help links. If that's the case, then you can still load

[R] Is latest R faster?

2009-04-20 Thread Ted Harding
Hi Folks, I just upgraded R to the latest R version 2.9.0 (2009-04-17) (from the Debian Etch repository). It seems to me that it at least starts up much faster than it used to and, though I don't have comparative timing data for long jobs, also that it does its work faster. Is that just an

Re: [R] graph with 15 combinations

2009-04-20 Thread Greg Snow
http://en.wikipedia.org/wiki/Venn_Diagram#Extensions_to_higher_numbers_of_sets shows a couple of solutions, not in R, but the ideas could be implemented in R. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111

Re: [R] bug when subtracting decimals?

2009-04-20 Thread Stephan Kolassa
Hi, Gavin Simpson wrote: I bemoan the apparent inability of those asking such questions to use the resources provided to solve these problems for themselves... Looking at all the people who quite obviously do NOT read the posting guide and provide commented, minimal, self-contained,

Re: [R] generic genotype calling algorithm?

2009-04-20 Thread Martin Morgan
Hi Nigel -- ask on the Bioconductor mailing list. http://bioconductor.org/docs/mailList.html Look at their package descriptions. http://bioconductor.org/packages/release/Software.html Martin Nigel Birney wrote: Hi, I have Agilent and Illumina SNP data. That's the only thing I know

Re: [R] graph with 15 combinations

2009-04-20 Thread Penner, Johannes
Thanks a lot! That is exactly what I was looking for! Best wishes Johannes -- Project Coordinator BIOTA West Amphibians Museum of Natural History Dep. of Research (Herpetology) Invalidenstrasse 43 D-10115 Berlin Tel: +49 (0)30 2093 8708 Fax: +49 (0)30 2093 8565 http://www.biota-africa.org

Re: [R] R-Squared with biglm?

2009-04-20 Thread Thomas Lumley
On Mon, 20 Apr 2009, Bryan Lim wrote: I've been working with a rather large data set (~10M rows), and while biglm works beautifully for generating coefficients, it does not report an r-squared. It does report RSS. Any idea on how one could coax an R-squared out of biglm? Hmm. I don't ever

[R] R 2.9.0 MASS package

2009-04-20 Thread Tom La Bone
I can't seem to find MASS for the latest version of R. Is it coming or has the name of the package changed? Tom -- View this message in context: http://www.nabble.com/R-2.9.0-MASS-package-tp23144198p23144198.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] R 2.9.0 MASS package

2009-04-20 Thread Peter Dalgaard
Tom La Bone wrote: I can't seem to find MASS for the latest version of R. Is it coming or has the name of the package changed? Tom Where did you look? It is in the sources (as part of VR), and also in my (SUSE) test builds. -- O__ Peter Dalgaard Øster Farimagsgade 5,

Re: [R] R 2.9.0 MASS package

2009-04-20 Thread Tom La Bone
In Windows Xp Pro: R2.8.1 USA(CA1) repository markerSearchPower MASS(VR) MasterBayes R2.9.0 USA(CA1) repository markerSearchPower MasterBayes MASS is not where it used to be. I checked a couple of other repositories in the US and got similar results. Tom Peter Dalgaard wrote: Tom

[R] Cross-Correlation function (CCF) issues

2009-04-20 Thread manta
Dear all, I have two series of returns and I want to find the cross-correlations between these two series. I know of the ccf, but it does not work as I'd like if i type ccf(x,y,lag.max=20,type=correlation,plot=FALSE) i got the error message Error in na.fail.default(ts.intersect(as.ts(x),

Re: [R] graph with 15 combinations

2009-04-20 Thread Emmanuel Charpentier
Le lundi 20 avril 2009 à 21:04 +0200, Penner, Johannes a écrit : Dear R helpers, I have a data set with 4 types (W, C, E S). Now I have values for all types plus all possible combinations (the order is unimportant): W, C, WC, E, WE, CE, WCE, S, WS, CS, WCS, ES, WES, CES WCES. Ideally I

  1   2   >