Re: [R] Transforming relational data

2011-02-18 Thread mathijsdevaan
OK, for the last step I have tried this (among other things): library(data.table) DT = data.table(read.table(textConnection(A B C 1 1 a 1999 2 1 b 1999 3 1 c 1999 4 1 d 1999 5 2 c 2001 6 2 d 2001 7 3 a 2004 8 3 b 2004 9 3 d 2004),head=TRUE,stringsAsFactors=FALSE))

Re: [R] Find peaks in dataset(x,y) and area for each peak

2011-02-18 Thread Laurent Gatto
Dear Ramya On 18 February 2011 05:08, Ramya ramya.vict...@gmail.com wrote: http://r.789695.n4.nabble.com/file/n3312061/x_and_y_values.txt x_and_y_values.txt I have the absorbance values form HPLC Chromatogram. I need to find the peaks in datapoints and area under each peak. I am not sure

Re: [R] color.scale error

2011-02-18 Thread Jim Lemon
On 02/18/2011 04:31 AM, Alaios wrote: Dear all when I call color.scale like this: require('plotrix') colcolor-color.scale(c(range_sr,sr),extremes=c(red,blue)) Error in if (min(reds) 0 || max(reds) 1) reds- rescale(reds, c(0, : missing value where TRUE/FALSE needed range_sr [1]

Re: [R] RGtk2 on Debian Testing

2011-02-18 Thread Jan van der Laan
It has been a while back, but I believe I had to install libgtk2.0-dev (that was on Ubuntu) You could also try to install the r-cran-rgtk2 debian-package using dpkg, aptitude, or whatever you use as package manager. This makes rgtk available for all users. HTH, Jan Quoting Lorenzo

Re: [R] Variable length datafile import problem

2011-02-18 Thread Ingo Reinhold
Hi John, seems there is no easy way. I'll just precondition it with AWK as described here http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg53401.html There are some remarks in the thread that R is not supposed to read too large files for political reasons. Maybe that's it. Many thanks

Re: [R] RGtk2 on Debian Testing

2011-02-18 Thread Jan van der Laan
It has been a while back, but I believe I had to install libgtk2.0-dev (that was on Ubuntu) You could also try to install the r-cran-rgtk2 debian-package using dpkg, aptitude, or whatever you use as package manager. This makes rgtk available for all users. HTH, Jan Quoting Lorenzo

Re: [R] How to change dataframe to tables

2011-02-18 Thread Dimitris Rizopoulos
say, 'Dat' is your data frame, then one way to do it is: with(Dat, tapply(freq, list(hair, eye, sex), c)) I hope it helps. Best, Dimitris On 2/18/2011 8:50 AM, Lao Meng wrote: The data is in the attachment. What I wanna get is: , , Sex = Male Eye HairBrown Blue Hazel Green

Re: [R] RGtk2 on Debian Testing

2011-02-18 Thread Lorenzo Isella
Thanks Jan, I am CCing to the list just to let everybody know that your suggestion pretty much fixes my problem. Installing RGtk2 was the last step to get rattle up and running. Somehow I had problems in running rattle with the RGtk2 package from the standard debian testing repositories (I

Re: [R] Dependency on R-Forge Package

2011-02-18 Thread Jari Oksanen
Damian Betebenner dbetebenner at nciea.org writes: In building a package, is it possible to make the package depend upon another package that is only available on R-Forge (not CRAN). For example, by doing something in the DESCRIPTION file I'd like to add a dependency to my package such that

Re: [R] need tick marks by 10s or 20s up to 200 on x axis

2011-02-18 Thread Jim Lemon
On 02/18/2011 07:56 AM, juliemcwh wrote: Sorry, I have read other posts on tick marks but nothing is working. Here is my R code. I was able to turn off the x axis at one point (I removed that code), but never got the tick marks on by 20s. It always goes to the default of 50, 100, 150, 200.

Re: [R] color.scale error

2011-02-18 Thread Alaios
Dear Jim, this is just to confirm that it your solution worked :) Regards Alex --- On Fri, 2/18/11, Jim Lemon j...@bitwrit.com.au wrote: From: Jim Lemon j...@bitwrit.com.au Subject: Re: [R] color.scale error To: Alaios ala...@yahoo.com Cc: R-help@r-project.org Date: Friday, February 18,

Re: [R] xyplot formula

2011-02-18 Thread Dennis Murphy
Hi: Here's another approach using the reshape package: library(reshape) # id defines the grouping variables; the others are stacked, creating # new variables 'variable' (containing the variable names as factor levels) and # 'value' (which holds the corresponding values) dfm - melt(df, id =

Re: [R] xyplot formula

2011-02-18 Thread Geoff Russell
On Fri, Feb 18, 2011 at 8:35 PM, Dennis Murphy djmu...@gmail.com wrote: Hi: Here's another approach using the reshape package: library(reshape) # id defines the grouping variables; the others are stacked, creating # new variables 'variable' (containing the variable names as factor levels)

Re: [R] How to change dataframe to tables

2011-02-18 Thread Dennis Murphy
This is a built-in dataset in R - see ?HairEyeColor and str() it. I smell homework... Dennis On Thu, Feb 17, 2011 at 11:50 PM, Lao Meng laomen...@gmail.com wrote: The data is in the attachment. What I wanna get is: , , Sex = Male Eye HairBrown Blue Hazel Green Black32

[R] lm without intercept

2011-02-18 Thread Jan
Hi, I am not a statistics expert, so I have this question. A linear model gives me the following summary: Call: lm(formula = N ~ N_alt) Residuals: Min 1Q Median 3Q Max -110.30 -35.80 -22.77 38.07 122.76 Coefficients: Estimate Std. Error t value Pr(|t|)

[R] line type in plots of the ade4 package

2011-02-18 Thread Duarte Viana
Hello all, In package ade4, does anyone know whether there is a simple way of setting the line type for the ellipses in an s.class plot? The argument lty does not work. Thanks, Duarte __ R-help@r-project.org mailing list

Re: [R] lm without intercept

2011-02-18 Thread Achim Zeileis
On Fri, 18 Feb 2011, Jan wrote: Hi, I am not a statistics expert, so I have this question. A linear model gives me the following summary: Call: lm(formula = N ~ N_alt) Residuals: Min 1Q Median 3Q Max -110.30 -35.80 -22.77 38.07 122.76 Coefficients:

Re: [R] How to change dataframe to tables

2011-02-18 Thread Henrique Dallazuanna
Try this: xtabs(Freq ~ Hair + Eye + Sex, Dat) Using Dimitri's Dat example. On Fri, Feb 18, 2011 at 5:50 AM, Lao Meng laomen...@gmail.com wrote: The data is in the attachment. What I wanna get is: , , Sex = Male Eye HairBrown Blue Hazel Green Black32 1110 3

Re: [R] censoring symbols on survfit plot

2011-02-18 Thread Terry Therneau
--begin included message Hi, when ploting Kaplan-Meier estimate curves as below, the censoring symbols (crosses) to not change thickness along the lines plot(survfit(surv ~ I(x=cut.off) ),lty=c(1,2), lwd=2) is there any strightforward way to make it happen? thanks -- End inclusion ---

Re: [R] lm without intercept

2011-02-18 Thread Dennis Murphy
Hi: On Fri, Feb 18, 2011 at 2:49 AM, Jan jrheinlaen...@gmx.de wrote: Hi, I am not a statistics expert, so I have this question. A linear model gives me the following summary: Call: lm(formula = N ~ N_alt) Residuals: Min 1Q Median 3Q Max -110.30 -35.80 -22.77

[R] multiple regression plane

2011-02-18 Thread Rosario Garcia Gil
Hello I have a multiple linear regression with two cofactors, I would like to represent a plane but I could not find any help which worked out. Any suggestions. Regards and thanks in advance. Rosario __ R-help@r-project.org mailing list

[R] Shared nearest neighbor (SNN) clustering algorithm implementation?

2011-02-18 Thread Jay
Hello, is there an implementation available for a shared nearest neighbor (SNN) clustering algorithm? //Jay __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] reading a matlab file

2011-02-18 Thread Krijnen WP, Wim
Dear R-people, I am trying to read a matlab file, unfortunately without success. I enclosed the zipped file. Do you have an idea how to read this file by ReadMat package into the R system? Thanks for your time! Wim Krijnen __

Re: [R] missing values in party::ctree

2011-02-18 Thread Torsten Hothorn
On Thu, 17 Feb 2011, Andrew Ziem wrote: After ctree builds a tree, how would I determine the direction missing values follow by examining the BinaryTree-class object? For instance in the example below Bare.nuclei has 16 missing values and is used for the first split, but the missing

Re: [R] error in optim, within polr(): initial value in 'vmmin' is not finite

2011-02-18 Thread Ben Haller
An update for the benefit of the list/posterity: I resolved this issue by switching over to using the lrm() function of package rms. It seems to pick better starts, or something; in any case, it has been able to converge on a solution for every model I've tried, although for the most complex

[R] Boundaries of R

2011-02-18 Thread Michael Holt
Hello Everyone, I'm pretty new to R and I'm trying to get some idea of the capabilities of the language. I work with some pretty large data sets and the consensus seems to be that R does not work well with big data. I've started talking to the guys at Revolution, but I need to get some outside

Re: [R] lm without intercept

2011-02-18 Thread Mike Marchywka
Date: Fri, 18 Feb 2011 12:25:36 +0100 From: achim.zeil...@uibk.ac.at To: jrheinlaen...@gmx.de CC: r-help@r-project.org Subject: Re: [R] lm without intercept On Fri, 18 Feb 2011, Jan wrote: Hi, I am not a statistics expert, so I have this

[R] how to incorporate prior base probabilities into binomial glmm

2011-02-18 Thread sprocket
Hi, I am modeling a gen linear mixed model with binomial repeated measures responses (y=cbind(correct/not correct)) of some individuals (id). I have a 2X2 design (effort and costs) and another continuous independent variable (scr). This will lead to something like this:

[R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread Maas James Dr (MED)
I'm attempting to sort a 3 dimensional array that looks like this x , , 1 [,1] [,2] [1,]99 [2,]79 , , 2 [,1] [,2] [1,]65 [2,]46 , , 3 [,1] [,2] [1,]21 [2,]32 Such that it ends up like this y , , 1 [,1] [,2] [1,]21

Re: [R] lm without intercept

2011-02-18 Thread Jay Emerson
No, this is a cute problem, though: the definition of R^2 changes without the intercept, because the empty model used for calculating the total sums of squares is always predicting 0 (so the total sums of squares are sums of squares of the observations themselves, without centering around the

Re: [R] reading a matlab file

2011-02-18 Thread Ken Knoblauch
Krijnen WP, Wim w.p.krijnen at pl.hanze.nl writes: Dear R-people, I am trying to read a matlab file, unfortunately without success. I enclosed the zipped file. Do you have an idea how to read this file by ReadMat package into the R system? Thanks for your time! Wim Krijnen

Re: [R] how to incorporate prior base probabilities into binomial glmm

2011-02-18 Thread Ken Knoblauch
sprocket toelch at gmail.com writes: Hi, I am modeling a gen linear mixed model with binomial repeated measures responses (y=cbind(correct/not correct)) of some individuals (id). I have a 2X2 design (effort and costs) and another continuous independent variable (scr). This will lead

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread Joshua Wiley
Dear Jim, This works for your small example. I was not exactly sure what you were sorting on. In my example, I sum all values at each level of the third dimension and sort by that sum. Of course, if there are many levels or the 2 dimensions being summed are quite large, this will be quite

Re: [R] lm without intercept

2011-02-18 Thread Jan
Hello Achim, Not quite. Consult your statistics textbook for the correct interpretation of p-values. Under the null hypothesis of a true intercept of zero, it is very likely to observe an intercept as large as 13.52 or larger. thank you for that help. I suppose the net doesn't have a

Re: [R] lm without intercept

2011-02-18 Thread Jan
Hi, thanks for your help. I'm beginning to understand things better. If you plotted your data, you would realize that whether you fit the 'best' least squares model or one with a zero intercept, the fit is not going to be very good Do the data cluster tightly around the dashed line? No, and

Re: [R] Boundaries of R

2011-02-18 Thread Duncan Murdoch
On 18/02/2011 5:44 AM, Michael Holt wrote: Hello Everyone, I'm pretty new to R and I'm trying to get some idea of the capabilities of the language. I work with some pretty large data sets and the consensus seems to be that R does not work well with big data. I've started talking to the guys at

[R] R script HELP!

2011-02-18 Thread whizevans
The following is my R script which I am struggling with to assess ICESat data..perhaps it is the ID_min or ID_max that is wrong? I don't know, any help would be greatly appreciated :( # OPTIONS - CHANGE THESE VARIABLES IF NEEDED\par

Re: [R] Boundaries of R

2011-02-18 Thread Mike Marchywka
Date: Fri, 18 Feb 2011 08:39:05 -0500 From: murdoch.dun...@gmail.com To: mich...@aers.ca CC: r-help@r-project.org Subject: Re: [R] Boundaries of R On 18/02/2011 5:44 AM, Michael Holt wrote: Hello Everyone, I'm pretty new to R and

Re: [R] R script HELP!

2011-02-18 Thread Jonathan P Daily
It looks to me like when you assign ID_min and ID_max that you want them to be vectors. In R, this is accomplished by ?c. Example: ID_min - c(976834501, 290018099.3145, 1000) -- Jonathan P. Daily Technician - USGS Leetown Science Center 11649 Leetown Road

Re: [R] Matrix in R

2011-02-18 Thread Dieter Menne
danielepippo wrote: I'm building a matrix in R with a cycle for like this: pp_ris2=matrix(NA,6,6) for(i in 0:6){ . R is not like c, indexing starts with 1. Dieter -- View this message in context: http://r.789695.n4.nabble.com/Matrix-in-R-tp3312748p3312764.html Sent from the

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread Claudia Beleites
Dear James, this is what I understood your sorting along the third dimension to be: x - array(c(9, 9, 7, 9, 6, 5, 4, 6, 2, 1, 3, 2), dim = list(2, 2, 3)) y - apply (x, 1:2, sort) y , , 1 [,1] [,2] [1,]21 [2,]65 [3,]99 , , 2 [,1] [,2] [1,]32 [2,]

Re: [R] Boundaries of R

2011-02-18 Thread Ben Bolker
Mike Marchywka marchywka at hotmail.com writes: Duncan Murdoch: Vectors are limited to about 2 billion entries (2^31 - 1). Matrices are vectors, so that limit applies to the total count of entries. Dataframes are lists of vectors, so that limit applies separately to the numbers of rows

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread Maas James Dr (MED)
Hi Rex, Thanks, this is exactly what I want but have to do it with many big arrays ... thus if there were a way to do it with a vectorized function would it not be a lot more efficient? Much appreciated! J Subject: RE: sort a 3 dimensional array across third dimension ? Although I suggested

Re: [R] Matrix in R

2011-02-18 Thread danielepippo
but if in my function pp_ris2[i,j]=myfunction} must be the indexes 0-0,0-1,0-2,0-3, ? -- View this message in context: http://r.789695.n4.nabble.com/Matrix-in-R-tp3312748p3312780.html Sent from the R help mailing list archive at Nabble.com. __

[R] Matrix in R

2011-02-18 Thread danielepippo
Hi everyone, I'm building a matrix in R with a cycle for like this: pp_ris2=matrix(NA,6,6) for(i in 0:6){ for(j in 0:6){ if(ij){ pp_ris2[i,j]=myfunction} else if(i==j){ print(c(i,j)) pp_ris2[i,j]=myfunction}

Re: [R] problem running scripts

2011-02-18 Thread jose Bartolomei
Dear Uwe, Thanks for your time. Definitely something was bugging my script (my mistake). Since your replay I have been trying to find the bug, to a simple script :( , without success. Preparing a reply to post with data sample and the script, everything start to work normally.

[R] ERGM + Anova

2011-02-18 Thread Roger Gomes
Hi All, I would like to know where to find technical papers to how calculate ANOVA using ERGM. Which parameters and so on. Thanks. -- View this message in context: http://r.789695.n4.nabble.com/ERGM-Anova-tp3312792p3312792.html Sent from the R help mailing list archive at Nabble.com.

[R] variance of model delta glm

2011-02-18 Thread Melen
hi everybody, i know it a quite complicate subject but someone might have the solution. I am doing a delta model coupling a binomial glm and a lognormal one. Using the Laurent correction I can predict mean values and I would like to know if you know how to predict the variance? Do you know the

[R] plot3d, color points by group

2011-02-18 Thread Kari
Hello I've created a 3d scatterplot, and had no problems labeling the points. However, I've been really struggling to change the color of the points based on a factor (see 'group' below). Is such a thing possible? My data look like this: food group x y z apple fruit 0.216

[R] Weights in bagged regression trees

2011-02-18 Thread Simon Gillings
Has anyone any experience of applying observational weights in bagging? I am performing regression trees (continuous data on bird abundance) and need to account for sampling intensity. In a single tree, i.e. a call of rpart, I can specify weights either by having a separate vector called

Re: [R] speed up the code

2011-02-18 Thread rex.dwyer
Yes, remove the call to intersect, and rely on the results of match to tell you whether there is an overlap. If there are any matches, all(is.na(index)) will be false. Read help for match. ?match -Original Message- From: r-help-boun...@r-project.org

[R] calculating means

2011-02-18 Thread Mike Gibson
I apologize if you have already seen my question. I am new to the mailing list and I did not get any responses the first time I posted my question. However, I am not sure my post went through. Here is my situation: I have a spreadsheet with columns of fish species (text) and length

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread rex.dwyer
Although I suggested to someone else that for-loops be avoided, they are not in the inner loop in this code, and it's probably easier to understand than some sort of apply: a = array(round(100*runif(60)),dim=c(3,4,5)) a for (i in 1:dim(a)[1]) for (j in 1:dim(a)[2]) a[i,j,] = sort(a[i,j,]) a

Re: [R] {Spam?} Re: sort a 3 dimensional array across third dimension ?

2011-02-18 Thread Claudia Beleites
On 02/18/2011 04:11 PM, Maas James Dr (MED) wrote: Hi Claudia, It does help a lot, but not quite there yet ... I'm sure you are correct and is much appreciated, I need some sort of generalized form, actual arrays in my case are 3x3x1000. Do you suspect it could be done in one step with sapply?

Re: [R] R script HELP!

2011-02-18 Thread whizevans
OK thanks Jonathon, I will have a go. Is the plotdir file technically my output file? Or what is this otherwise? Also, having never used R script before, is it not very good at reading data, may I have to alter my input file in excel first? Thanks -- View this message in context:

Re: [R] calculating means

2011-02-18 Thread Bert Gunter
Have you read an Intro to R or, indeed, any R documentation? Perhaps you received no reply because you do not appear to have made an honest effort to answer the question for yourself. Maybe you should consider accessing R thru R Commander or one of the other R GUI's -- search CRAN for details. --

[R] Variogram (nlme) of a lme object - corSpatial element question

2011-02-18 Thread Arnaud Mosnier
Dear Users, From previous analysis (semi-variograms using package gstat), I found spatial autocorrelation in my dataset. The best fitted model to this spatial correlation structure is the Gaussian model (Spherical, Exponential, Linear tested and comparison done by Sum of Square errors). So I used

Re: [R] Matrix in R

2011-02-18 Thread Philipp Pagel
On Fri, Feb 18, 2011 at 06:32:01AM -0800, danielepippo wrote: but if in my function pp_ris2[i,j]=myfunction} must be the indexes 0-0,0-1,0-2,0-3, ? You'll have to take care of that yourself with a bit of index arithmetics. It's the same you encounter in C, if you are modelling

Re: [R] ERGM + Anova

2011-02-18 Thread Bert Gunter
What does this have to do with R? -- Bert On Fri, Feb 18, 2011 at 6:37 AM, Roger Gomes roger.mestr...@gmail.com wrote: Hi All, I would like to know where to find technical papers to how calculate ANOVA using ERGM. Which parameters and so on. Thanks. -- View this message in context:

Re: [R] Matrix in R

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 9:32 AM, danielepippo wrote: but if in my function pp_ris2[i,j]=myfunction} must be the indexes 0-0,0-1,0-2,0-3, ? From a search of RSiteSearch() that started with with terms: zero matrix indexing http://finzi.psych.upenn.edu/R/Rhelp02/archive/39031.html

Re: [R] Matrix in R

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 9:32 AM, danielepippo wrote: but if in my function pp_ris2[i,j]=myfunction} must be the indexes 0-0,0-1,0-2,0-3, ? I came across a posting in r-help that called this package blasphemy: http://cran.r-project.org/web/packages/Oarray/index.html -- David Winsemius,

[R] Bug in the function Variogram (package nlme) ?

2011-02-18 Thread Arnaud Mosnier
Dears UseRs, In Zuur's book (Mixed effects models and extensions in ecology with R), p 167, it is mentionned that the Variogram function (package nlme) scale the sill to 1. Here is the plot of the semi-variogram given by the Variogram function on my lme object. (http://imm.io/3OLe) Note that the

Re: [R] R script HELP!

2011-02-18 Thread Jonathan P Daily
Well, this script only shows the option settings, so I assumed that the string being used as your infile was formatted for R. As for how it will have to look for R to read it in, there are many ways to do this - I would start by reading R Data Import/Export that ships with every copy of R. As

Re: [R] Dependency on R-Forge Package

2011-02-18 Thread Spencer Graves
I echo Prof. Ripley's and Jari Oksanen's concerns. I see two options in case the R-Forge package maintainer does not submit the package to CRAN: (1) Copy the function(s) you want from that package into your own. In such cases, I recommend you change the name of the

[R] How to flag those iterations which yield a warning?

2011-02-18 Thread Ravi Varadhan
Hi, I am running a simulation study with the survival::coxph. Some of the simulations result in problematic fits due to flat partial likelihood. So, you get the warning message: Warning message: In fitter(X, Y, strats, offset, init, control, weights = weights, ... : Loglik converged before

Re: [R] How to flag those iterations which yield a warning?

2011-02-18 Thread Ravi Varadhan
May be I could do: for (i in 1:nsim) { last.warning - NULL # do model fitting if(!is.null(last.warning)) # discard simulation result } I think this might work. Any other ideas? Ravi. Ravi Varadhan, Ph.D. Assistant

[R] heatmap cluster

2011-02-18 Thread Duke
Hi all, I am new to heatmap function in stats package, so any advice will be greatly appreciated. So I have a dataset that I want to draw a heatmap, for example: dput(a) structure(c(-8.256023, -14.764429, -10.145779, -5.123104, -11.974729, -16.329029, -16.329029, -16.329029, -2.764114,

Re: [R] How to flag those iterations which yield a warning?

2011-02-18 Thread Ravi Varadhan
Nope - that does not work. The value of last.warning is not reset after the initial NULL. Ravi. Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins

Re: [R] How to flag those iterations which yield a warning?

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 11:35 AM, Ravi Varadhan wrote: May be I could do: for (i in 1:nsim) { last.warning - NULL # do model fitting if(!is.null(last.warning)) # discard simulation result You can raise warnings to behave as errors with options: options(warn=2) And you can define where R

Re: [R] How to flag those iterations which yield a warning?

2011-02-18 Thread William Dunlap
If you have a function that may produce warnings you can use withCallingHandlers() to make a new function that does what the original does but adds some contextual information to the usual warning report. E.g., addWarningReporting - function(FUN) { force(FUN) function(args) {

Re: [R] calculating means

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 9:45 AM, Mike Gibson wrote: I apologize if you have already seen my question. I am new to the mailing list and I did not get any responses the first time I posted my question. However, I am not sure my post went through. It didn't. Here is my situation: I have a

[R] fit a distribution

2011-02-18 Thread Simone Gabbriellini
hello list, I have this array of distribution of frequencies: freq-c(45,18,10,8,13,5,9,3,4,1,2,2,2,3,4,4,1,1,2,1,2,1,1,1,1,1,2,3,1,1,1,1,1,1,1,1,2,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1) for this array of values:

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread jamaas
Thanks Claudia (and all others!0 This works beautifully! The instructions for aperm are a bit vauge! J -- View this message in context: http://r.789695.n4.nabble.com/sort-a-3-dimensional-array-across-third-dimension-tp3312612p3313126.html Sent from the R help mailing list archive at

[R] hdf5 library install issue

2011-02-18 Thread Jerome
Dear all I'm trying to install the hdf5 library into R. The HDF5 package is installed in a special directory, distributed accross my cluster: /share/apps/HDF5 So i put the args option to the configure command as i read in previous post in the list: jerome]# R CMD INSTALL

Re: [R] sort a 3 dimensional array across third dimension ?

2011-02-18 Thread rex.dwyer
I was going to say: The problem with for-loops (as best I understand it) is that the R code gets interpreted over and over; what you normally want to do is design the computation so that you jump into the internals of R and stay there. But the inner loop is in the R internals of the sort in

Re: [R] ERGM + Anova

2011-02-18 Thread Roger Gomes
Hi Bert, thanks about your answer. Well, I would like to calculate anova.ergm in R. Wich values (edges, vertices) the library use to calculate F values or p-value with graphos. The calculus is like a Pandora Box... I put there a network (ergm) and then, values show me like a magic. Do you

[R] gaoshan

2011-02-18 Thread wang peter
i have got a svm MODEL from e1071 by svmModel-svm(x, y, kernel = linear, cost = gamma, scaled=FALSE) who can tell me how can i get the alpha of decision functions i just remembered rho is the bias of decision functions i think the manual of e1071 is a piece of usely garbage thank u gao shan

[R] Is the Diebold Mariano Test in forecast package adjusted?

2011-02-18 Thread Paka yag
Hello I would like to know if the Diebold Mariano Test in the forecast Package is adjusted to small samples (as Harvey, Leybourne, Newbold suggest) If not, how can I do that manually? Paka __ R-help@r-project.org mailing list

Re: [R] plot3d, color points by group

2011-02-18 Thread Dieter Menne
Kari-4 wrote: I've created a 3d scatterplot, and had no problems labeling the points. However, I've been really struggling to change the color of the points based on a factor (see 'group' below). Is such a thing possible? Its always polite to make the example self-running by adding

[R] How to calculate the perimeter and common border of polygons?

2011-02-18 Thread Leonardo Monasterio
Dear R-users, Is there any way of calculating the perimeter of a polygon in a shapefile object? Furthermore, how to calculate the length of the common border of two polygons? I've searched the code of spded, but could not find a hint on how to do it. Thank you very much, Leo.

[R] question of elimination drawn lines on competing risk graph

2011-02-18 Thread Gerard Smits
Hi All, I am using the package, cmprisk, to plot competing risks. In my case, I have four lines showing risk of going on dialysis (by a lab test [fgf-23] in quartiles), where the 4 lower lines are for the competing risk of death. I am trying to edit the function to just plot the upper curves

[R] segfault during example(svm)

2011-02-18 Thread Juergen Rose
If do: library(e1071) example(svm) I get: svm data(iris) svm attach(iris) svm ## classification mode svm # default with factor response: svm model - svm(Species ~ ., data = iris) svm # alternatively the traditional interface: svm x - subset(iris, select = -Species) svm y - Species svm

[R] Bootstraps standard error

2011-02-18 Thread danielepippo
Hi everyone, if I have the bootsrap values of the parameters of my regression model, how can I calculate the standard error of them to refer to the parameters? Thank you -- View this message in context: http://r.789695.n4.nabble.com/Bootstraps-standard-error-tp3313322p3313322.html Sent

[R] Scaling Lattice Graphics for tikzDevice

2011-02-18 Thread Elliot Joel Bernstein
I'm trying to use lattice graphics to produce some small plots for inclusion in a LaTeX file. I want the LaTeX fonts to be used in the plots, but to be scaled down to match the size of the plot. I have written the following code to apply a scaling factor to all the cex and padding entries in

[R] non-ascii characters in R output

2011-02-18 Thread Matt Shotwell
All, I'd like to automatically output text from R to HTML. In doing this I've run into trouble with non-ascii characters, as my browser (and presumably others) does not render such characters correctly. For example, the 'fancy' single quotes associated with summary.lm are multi-byte

[R] clustered bar chart help

2011-02-18 Thread Li, Qinghong, ST. LOUIS, NRC St. Louis Petcare
Hi, Can anyone help me to plot a chart graph? I have a data set like this. I would like a bar chart graph which mouse1's two treatments are clustered together, so on. I tried with barplot, but couldn't get it right. Thanks Johnny Mouse1 Mouse2 Mouse3

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

2011-02-18 Thread Tom Hopper
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, and to interactively create call-outs of sections of the data. Much of the data that I deal with

Re: [R] Bootstraps standard error

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 1:36 PM, danielepippo wrote: Hi everyone, if I have the bootsrap values of the parameters of my regression model, how can I calculate the standard error of them to refer to the parameters? The same way you would any other vector? Thank you -- View this message

Re: [R] segfault during example(svm)

2011-02-18 Thread Peter Ehlers
On 2011-02-18 11:16, Juergen Rose wrote: If do: library(e1071) example(svm) I get: svm data(iris) svm attach(iris) svm ## classification mode svm # default with factor response: svm model- svm(Species ~ ., data = iris) svm # alternatively the traditional interface: svm x-

Re: [R] question of elimination drawn lines on competing risk graph

2011-02-18 Thread David Winsemius
On Feb 18, 2011, at 2:10 PM, Gerard Smits wrote: Hi All, I am using the package, cmprisk, to plot competing risks. In my case, I have four lines showing risk of going on dialysis (by a lab test [fgf-23] in quartiles), where the 4 lower lines are for the competing risk of death. I am

Re: [R] segfault during example(svm)

2011-02-18 Thread Juergen Rose
Am Freitag, den 18.02.2011, 11:53 -0800 schrieb Peter Ehlers: On 2011-02-18 11:16, Juergen Rose wrote: If do: library(e1071) example(svm) I get: svm data(iris) svm attach(iris) svm ## classification mode svm # default with factor response: svm model- svm(Species ~

Re: [R] segfault during example(svm)

2011-02-18 Thread Claudia Beleites
Dear Jürgen, did you update.packages (checkBuilt = TRUE) ? I recently had segfaults, too on 64bit linux (with rgl, though) and they disappeared only after updating with checkBuilt (including also the packages originally installed via Dirk's .deb packages. HTH, Claudia On 02/18/2011 09:32

Re: [R] segfault during example(svm)

2011-02-18 Thread Peter Ehlers
On 2011-02-18 12:32, Juergen Rose wrote: Am Freitag, den 18.02.2011, 11:53 -0800 schrieb Peter Ehlers: On 2011-02-18 11:16, Juergen Rose wrote: If do: library(e1071) example(svm) I get: svm data(iris) svm attach(iris) svm ## classification mode svm # default with factor

Re: [R] segfault during example(svm)

2011-02-18 Thread Juergen Rose
Am Freitag, den 18.02.2011, 21:42 +0100 schrieb Claudia Beleites: update.packages (checkBuilt = TRUE) Dear Claudia, I just tried update.packages (checkBuilt = TRUE), it did not solved the issue. Nevertheless thanks, Juergen __ R-help@r-project.org

[R] Means

2011-02-18 Thread Dmitry Berman
Listers, Is there a command/function to get the population standard deviation (N) and the sample standard deviation (n-1) Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Means

2011-02-18 Thread Ted Harding
On 18-Feb-11 20:51:40, Dmitry Berman wrote: Listers, Is there a command/function to get the population standard deviation (N) and the sample standard deviation (n-1) Thanks Using data (1:10), with N=10, as an example: c(SampSD=1, PopSD=sqrt(9/10))*sd((1:10)) # SampSDPopSD #

Re: [R] segfault during example(svm)

2011-02-18 Thread Claudia Beleites
Dear Jürgen, I just tried update.packages (checkBuilt = TRUE), it did not solved the issue. Was just a guess. Nevertheless thanks, Bitte, bitte. In case that helps tracking down the problem: it works with my normal R as well as with r-devel. In case it helps you, I also put the

[R] breaks/bins question

2011-02-18 Thread Julie McWhorter
Thank you Joshua and Jim--I got the tick marks on finally! One last question: I am supposed to make bins in 1, 2 and 5 mm increments for fish up to 170 mm. I assume these are the breaks? So breaks = 170 for 1 mm and 340 for 2 mm and 1550 for 5 mm? I couldn't find any specific examples in

Re: [R] clustered bar chart help

2011-02-18 Thread Li, Qinghong, ST. LOUIS, NRC St. Louis Petcare
Sorry the data were all messed up. Let me try again. Can anyone help me to plot a chart graph? I have a data set like this. I would like a bar chart graph which mouse1's two treatments are clustered together, so on. I tried with barplot, but couldn't get it right. Mouse1 Mouse2 Mouse3

[R] R2wd

2011-02-18 Thread rookie r user
Hi all, I'm currently trying to get accustomed to R and some of its packages, thus need some help for a package called R2wd. As I am trying to use the wdGET() function I'm getting the following message: Error in if (!(tmp[[ActiveDocument]][[Name]] == filename)) tmp$Open(paste(path, :

  1   2   >