Re: [R] vegdist Error en double(N * (N - 1)/2) : tama?o del vector especificado es muy grande

2013-02-09 Thread Prof Brian Ripley
Suppose N = 138037 (you haven't really told us). A dissimilarity half-matrix would have 9 billlion elements. The maximum size of a vector in current versions of R is 2 billion. You will be able to get further in R-devel (3.0.0-to-be) with a 64-bit version of R, although as you appear to be

Re: [R] help with double looping

2013-02-09 Thread christel lacaze
Hi Rui, thanks very much for your help - i've finally cracked it... the correct code was iscores-list()for (j in seq(1,166,11)){for (i in j:(j+10)) {iscores[[i]=(i-j)*data[[i]] }} instead of iscores[[i-j+1]]=(i-j)*data[[i-j+1]] I probably wasn't explaining myself very well if you're

[R] Addressing Columns in a Data Frame

2013-02-09 Thread Lorenzo Isella
Dear All, Probably a one liner, but I am banging my head against the floor. Consider the following DF - data.frame( x=1:10, y=10:1, z=rep(5,10), a=11:20 ) mn-names(DF) but then I cannot retrieve a column by doing e.g, DF$mn[2] I tried to play with the quotes and so on, but so

Re: [R] Addressing Columns in a Data Frame

2013-02-09 Thread Michael Weylandt
On Feb 9, 2013, at 1:40 PM, Lorenzo Isella lorenzo.ise...@gmail.com wrote: Dear All, Probably a one liner, but I am banging my head against the floor. Consider the following DF - data.frame( x=1:10, y=10:1, z=rep(5,10), a=11:20 ) mn-names(DF) but then I cannot retrieve

Re: [R] Addressing Columns in a Data Frame

2013-02-09 Thread Denis Francisci
I don't know if I understood your problem, but maybe you can retrieve your columns simply in this way: DF$y [1] 10 9 8 7 6 5 4 3 2 1 or (for looking at column) DF[2] y 1 10 2 9 3 8 4 7 5 6 6 5 7 4 8 3 9 2 10 1 Bye D. 2013/2/9 Lorenzo Isella

[R] help on plotting series with different x-axis values on a graph with one x-axis

2013-02-09 Thread Eike Marie Thaysen
Hi, I have run into this problem a couple of times now and hope you can help! If I want to plot mulitiple series with differing x-axis values (however, all in the same range) in 1 plot with one common axis R obstruses the plots. E.g. if I plot water content against time and I start with the

Re: [R] help on plotting series with different x-axis values on a graph with one x-axis

2013-02-09 Thread David Winsemius
On Feb 9, 2013, at 4:52 AM, Eike Marie Thaysen wrote: Hi, I have run into this problem a couple of times now and hope you can help! If I want to plot mulitiple series with differing x-axis values (however, all in the same range) in 1 plot with one common axis R obstruses the plots.

[R] Swaeve, Beamer and \alt

2013-02-09 Thread Jarrod Hadfield
Hi, I am having trouble getting \alt (or \altenv) to work with Schunk/Sinput and was wondering if anybody had had success? With the slide \begin{frame}[fragile]\frametitle{Basic R} \alt2{ echo=TRUE= 2+2 @ }{ echo=TRUE, eval=FALSE= 2+2 @ } \end{frame} I get the error message: ! FancyVerb

Re: [R] Troubleshooting underidentification issues in structural equation modelling (SEM)

2013-02-09 Thread John Fox
Dear Ruijie, Your model is underidentified by virtue of two of the factors having only one observed indicator each. No SEM software can magically estimate this model as it stands. Beyond that, I won't comment on the wisdom of what you're doing, such as computing covariances between ordinal

Re: [R] Addressing Columns in a Data Frame

2013-02-09 Thread Jim Holtman
Try DF [[mn [2]]] Sent from my Verizon Wireless 4G LTE Smartphone Original message From: Lorenzo Isella lorenzo.ise...@gmail.com Date: 02/09/2013 8:40 AM (GMT-05:00) To: r-h...@stat.math.ethz.ch Subject: [R] Addressing Columns in a Data Frame Dear All, Probably a one

Re: [R] qcc package

2013-02-09 Thread Uwe Ligges
On 08.02.2013 23:37, Steven LeBlanc wrote: Greets, My data looks like: p3.18 s xbar subgroup 1 0.84 12.21 2 1.64 11.22 3 2.07 10.63 4 2.49 12.24 5 0.84 11.25 ... Using the command

Re: [R] Swaeve, Beamer and \alt

2013-02-09 Thread Uwe Ligges
On 09.02.2013 17:00, Jarrod Hadfield wrote: Hi, I am having trouble getting \alt (or \altenv) to work with Schunk/Sinput and was wondering if anybody had had success? With the slide \begin{frame}[fragile]\frametitle{Basic R} \alt2{ echo=TRUE= 2+2 @ }{ echo=TRUE, eval=FALSE= 2+2 @ }

Re: [R] Swaeve, Beamer and \alt

2013-02-09 Thread Jarrod Hadfield
Hi, Unfortunately not, nor \begin{frame}[fragile], which does work with the uncover environment. Cheers, Jarrod Quoting Uwe Ligges lig...@statistik.tu-dortmund.de on Sat, 09 Feb 2013 18:45:39 +0100: On 09.02.2013 17:00, Jarrod Hadfield wrote: Hi, I am having trouble getting \alt

Re: [R] funcion ff

2013-02-09 Thread Uwe Ligges
Two comments: This list is in English. The package maintainer may know better answers than this list. Best, Uwe Ligges On 05.02.2013 15:21, Tania Patiño wrote: Hola R, tengo las siguientes preguntas: Pregunta 1: Cargar las tablas de los datos de peliculas en R usando `ff`. Cómo se

Re: [R] plot significant spearman correlation

2013-02-09 Thread Uwe Ligges
On 05.02.2013 17:08, Ozgul Inceoglu wrote: Hello, I calculate the correlation between two matrices cor(x,y, method=spearman) and I am wondering if it possible to see only the significant correlations. I can do that for single OTUs with cor.test() command but I would like to have an output

Re: [R] maximum amount of lines of code to be copy pasted to R console

2013-02-09 Thread Uwe Ligges
On 05.02.2013 20:30, Jannis wrote: Dear list members, is there any maximum amount of lines of code that can be copy pasted into the R console? Or some maximum amount of levels of nested parantheses? Whenever I copy paste some big chunk of code into the R console I get errors of missing

Re: [R] qcc package

2013-02-09 Thread Uwe Ligges
On 09.02.2013 19:10, Steven LeBlanc wrote: Hi Uwe, Thanks much. I don't think this is the issue. My entire data set is: p3.18 s xbar subgroup 1 0.84 12.21 2 1.64 11.22 3 2.07 10.63 4 2.49 12.24 5 0.84 11.25 6 1.82 12.66 7 2.95

Re: [R] character strings with embedded commands: perl /gee ?

2013-02-09 Thread ivo welch
joshua wiley told me how to do this. so, here is something I find quite useful: estring - function (e2) { syntax - regmatches(e2, gregexpr((?=\\{\\{).*?(?=\\}\\}), e2, perl=TRUE))[[1]] res - lapply(syntax, function(p) eval(parse(text=p))) res - lapply(1:length(res), function(i) { e2

Re: [R] impossible to invert a spam-object, but possible when it's a matrix-object

2013-02-09 Thread Uwe Ligges
On 05.02.2013 13:29, Camarda, Carlo Giovanni wrote: Dear R-users, a question concerning sparse matrices in package spam (spam_0.29-2). On one hand I have a spam object (n X n) from which I cannot compute the inverse. On the other hand, if I convert this object in a plain matrix, I can find

Re: [R] order function

2013-02-09 Thread Uwe Ligges
On 01.02.2013 17:52, Kripa R wrote: Hi I'm having a simple issue with the order function. When I use the following code my data is not ordered correctly output[order(output[,3]),] Namebeta pval 881 9.09303277751237 0.000100253350483199 7402

Re: [R] Error: package/namespace load failed for ‘arm’

2013-02-09 Thread Uwe Ligges
On 01.02.2013 12:29, Pernille Heelsberg wrote: Hi! I get the following error messages when trying to run the package 'arm'. library(arm) Error : Functions found when exporting methods from the namespace ‘arm’ which are not S4 generic: ‘fixef’, ‘ranef’ Error: package/namespace load failed

Re: [R] Saving history across sessions

2013-02-09 Thread Uwe Ligges
On 07.02.2013 19:21, James Jong wrote: I have an environment variable `$R_HISTFILE` pointing to `/home/my.username/.RHistory` and the following code in my `.Rprofile` in my home directory: .Last - function() { if (!any(commandArgs()=='--no-readline') interactive()){

Re: [R] question about reproducibility/consistency of principal component and lda directions in R

2013-02-09 Thread Uwe Ligges
On 08.02.2013 20:14, David Romano wrote: Hi everyone, I'm not exactly sure how to ask this question most clearly, but I hope that giving the context in which it occurs for me will help: I'm trying to compare the brain images of two patient populations; each image is composed of voxels (the

Re: [R] Feature selection for kmeans

2013-02-09 Thread Uwe Ligges
On 07.02.2013 19:02, SpaceSeller wrote: I know that within sum of squares, DB, sillhouette and cophenetic are indicators of clustering quality, but what indicators I need to observe when I choose attributes for kmeans? There are lots of papers about ways to measure the quality. For example,

Re: [R] help with creating new variables using a loop

2013-02-09 Thread Uwe Ligges
On 07.02.2013 17:16, christel lacaze wrote: Hi there, I've got a set of 10 numeric variables called Mood1 to Mood10 in a dataset called mood. I'm trying to create a set of 10 new variables called m1 to m10 so that m1=Mood1*1, m2=Mood2*2, etc to m10=Mood10*10 Trawling through the

Re: [R] Empty cluster / segfault using vanilla kmeans with version 2.15.2

2013-02-09 Thread Uwe Ligges
We need a reproducible example. Uwe Ligges On 03.02.2013 15:03, Luca Nanetti wrote: Dear experts, I am encountering a version-dependent issue. My laptop runs Ubuntu 12.04 LTS 64-bit, R 2.14.1; the issue explained below never occurred with this version of R My desktop runs Ubuntu 11.10

[R] Optimization Problem in R

2013-02-09 Thread Axel Urbiz
Dear List, I'm new in R. I'm trying to solve a simple constrained optimization problem. Essentially, let's say I have a matrix as in the object 'mm' inside the function below. My objective function should have a matrix of parameters, one parameter for each element 'mm' (4 in this case). The

[R] R maps Help

2013-02-09 Thread Frank Rodgers
I am fairly new to R and am plotting flight data on a map. Everything is working well except the size of the map is really too small to show the data effectively and I can't seem to figure out how to make the output map larger. Do I need to change the device characteristics or is it a map.???

Re: [R] help with creating new variables using a loop

2013-02-09 Thread arun
Hi,  mood- data.frame(Mood1=1:5,Mood2=6:10,Mood3=11:15) lapply(seq_along(mood),function(i) i*mood[i]) A.K. - Original Message - From: christel lacaze christellac...@hotmail.co.uk To: r-help@r-project.org Cc: Sent: Thursday, February 7, 2013 11:16 AM Subject: [R] help with creating new

Re: [R] Optimization Problem in R

2013-02-09 Thread Berend Hasselman
On 09-02-2013, at 21:08, Axel Urbiz axel.ur...@gmail.com wrote: Dear List, I'm new in R. I'm trying to solve a simple constrained optimization problem. Essentially, let's say I have a matrix as in the object 'mm' inside the function below. My objective function should have a matrix of

Re: [R] confidence / prediction ellipse

2013-02-09 Thread Michael Friendly
Hi Giuseppe You've posted a series of questions on this topic, never with any code you've tried or data, and usually with some undefined references to 'a' and 'b'. In spite of this, a variety of people have tried to give you helpful replies, intuiting what it is you might have meant. This is

[R] Questions about caret - Best place to ask?

2013-02-09 Thread James Jong
Hi there. I am relatively new to R and to the community. I have previously asked my R questions at R-help (through Nabble, not sure if there is a better way other than Nabble or regular email), and StackOverflow with the [R] tag. I have some questions that are somewhat specific about the

[R] Correct value of R_HOME

2013-02-09 Thread James Jong
What is the purpose of the variable R_HOME? What value should I set it to? Thanks, James [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] rpy2 (Python -- R)

2013-02-09 Thread James Jong
Hi - I am having trouble installing rpy2. I have already compiled R as a shared library, but I do not have admin privileges on so I am trying to install rpy2 with: pip install -user rpy2 but I am getting the following error:

Re: [R] First R package, advice

2013-02-09 Thread ivo welch
I am still researching my best choices. IMHO, every programming and every typesetting language ever invented should have had a compatible basic escape sequence built in. simple. just compatible with one another. something like a specific escape sequence on the first column, like #%@@==R

Re: [R] R maps Help

2013-02-09 Thread Ray Brownrigg
I don't understand what you mean by too small, and you haven't provided a reproducible example. The map will be as large as the default graphics window on your screen will allow (leaving room for axes and labels as well). If you are not happy with this, try something like: x11(width=30,

Re: [R] help with creating new variables using a loop

2013-02-09 Thread arun
Hi, Just to add: If the columns are  not in the order (Mood1:Mood10) in the dataset.  mood- data.frame(Mood2=6:10,Mood1=1:5,Mood10=11:15,Mood3=16:20)  do.call(cbind,lapply(names(mood),function(i) mood[i]*as.numeric(gsub(\\D+,,i #  Mood2 Mood1 Mood10 Mood3 #1    12     1    110    48 #2