[R] Distance used in lda()

Whenever using lda(...,CV=TRUE), or if using lda.predict(), a classification is returned. Which distance is used in order to allocate a given individual to a class, once the original data have been transformed into discriminant components? Thanks Agustin Lobo aloboa...@gmail.com

[R] Problem with system() and source on linux

ot found" from within system()? Any help much appreciated, Thanks -- Agustin Lobo aloboa...@gmail.com __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] How to create a list of trellis objects for grid.arrange()

ects p1 and p2 within one single object (as in practice I have many objects generated within a for() loop). The following used to work: ps <- c(p1,p2) grid.arrange(ps,ncol=1,nrow=2) but does not work any more. How should I combine p1 and p2 into one single object that would be accepted by grid.arran

[R] R-square prob is not calculated by randomization in lmPerm::lmp

9e+05 on 1 and 998 DF, p-value: < 2.2e-16 -- Agustin Lobo aloboa...@gmail.com __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/post

[R] Friedman test

46,0.85,0.48,0.78,0.87,0.22,0.51,0.03,0.39,0.44) midata <- data.frame(treatments,people,y) midata[1:10,] I would do as follows: Test among treatments friedman.test(y ~ treatments|people,data=midata) Test among subjects friedman.test(y ~ people|treatments,data=midata) Is this correct? Thanks -

[R] margins and device: different results with pdf() and svg()

labeling_0.2 [7] lattice_0.20-24munsell_0.4.2 proto_0.3-10 RColorBrewer_1.0-5 reshape2_1.2.2 scales_0.2.3 [13] stringr_0.6.2 tools_3.0.2 -- Agustin Lobo aloboa...@gmail.com __ R-help@r-project.org mailing list https

[R] Different output from lm() and lmPerm lmp() if categorical variables are included in the analysis

I've found a problem when using categorical variables in lmp() from package lmPerm According to help(lmp): This function will behave identically to lm() if the following parameters are set: perm=, seq=TRUE, center=FALSE.) But not in the case of including categorical variables: require(lmPerm)

[R] Problem with lubridate and decimal seconds

] stats graphics grDevices utils datasets methods base other attached packages: [1] lubridate_1.3.0 loaded via a namespace (and not attached): [1] digest_0.6.3 memoise_0.1 plyr_1.8 stringr_0.6.2 tools_3.0.0 -- Agustin Lobo aloboa...@gmail.com

[R] query labels in iplot() (or other interactive scatterplot)

) or with identify(x,y,labels), but I do not want to get the labels permanently plotted on the graphic) Thanks -- -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 e-mail agustin.l...@ictja.csic.es

[R] strptime() with format %OS does not print millisecs in MacOS

tools_2.15.1 -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 e-mail agustin.l...@ictja.csic.es https://sites.google.com/site/aloboaleu/ __ R-help@r

[R] How to deal with thousands of seconds in R?

If I do: a [1] 2012_10_01_14_13_32.445 a2 [1] 2012_10_01_14_13_32.500 strptime(a,format=%Y_%M_%d_%H_%M_%S)-strptime(a2,format=%Y_%M_%d_%H_%M_%S) Time difference of 0 secs Is there any time object in R that would deal with thousands of seconds? Thanks Agus -- -- Dr. Agustin Lobo Institut de

Re: [R] How to deal with thousands of seconds in R?

Sorry for the typo, I meant thousandths, and thanks for pointing out the %OS format, which I had overlooked Agus On Tue, Oct 9, 2012 at 2:51 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: On 09/10/2012 09:54, Agustin Lobo wrote: If I do: a [1] 2012_10_01_14_13_32.445 a2 [1

[R] day of the year for chron objects

? a - chron(20100506,format=ymd) a [1] 100506 years(a) [1] 2010 Levels: 2010 but instead have to cast to Date to apply julian: julian(a) Error in names(d) : 'd' is missing julian(as.Date(a)) [1] 14735 attr(,origin) [1] 1970-01-01 Thanks Agus -- -- Dr. Agustin Lobo Institut de Ciencies de la Terra

Re: [R] day of the year for chron objects

, are all valid; Thanks! Agus On Fri, Jun 8, 2012 at 12:54 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Fri, Jun 8, 2012 at 4:01 AM, Agustin Lobo agustin.l...@ictja.csic.es wrote: Hi! Is not there an standard R function to retrieve the day of the year (since 1st Jan of the same year

[R] problems with xlim in plot.zoo() and window()

://dl.dropbox.com/u/3180464/alyL32007z.rda http://dl.dropbox.com/u/3180464/alyL32007.rda alyL32007z - zoo(alyL32007[-1],alyL32007$time) Agus -- -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 e-mail

Re: [R] problems with xlim in plot.zoo() and window()

-04 01:30:00 2007-03-04 02:00:00 2007-03-04 02:30:00 ... Would you advice in favor of using Date objects instead of POSIXct ? Agus On Tue, May 29, 2012 at 7:13 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Tue, May 29, 2012 at 11:47 AM, Agustin Lobo agustin.l...@ictja.csic.es wrote

[R] zoo: variable gets modified at making zoo object

Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 e-mail agustin.l...@ictja.csic.es https://sites.google.com/site/aloboaleu/ __ R-help@r-project.org

Re: [R] zoo: variable gets modified at making zoo object

Thanks for your interest. I've put the dataframe alyL32007 in http://dl.dropbox.com/u/3180464/alyL32007.rda ready to be used with load() Agus On Mon, May 28, 2012 at 6:11 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Mon, May 28, 2012 at 5:35 AM, Agustin Lobo agustin.l

[R] Problem with zoo::window()

in the manual: R window(z, start = as.Date(2005-02-15), end = as.Date(2005-02-28)) Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris S/N 08028 Barcelona Spain Tel. +34 934095410 Fax. +34 934110012 email: agustin.l...@ictja.csic.es

[R] cutree() and rect.hclust(): different labelling of classes

it be possible that future versions provide the same labeling? rect.hclust() is useful to display the classes over the dendrogram but the output from cutree() is often more handy. Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris S/N 08028 Barcelona Spain

[R] as.numeric() and POSIXct format

function of as.POSIXct()? I've tried as.numeric and unclass() using a origin= argument, but this does not work. Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: agustin.l

[R] heatmap with values

Hi! I'm displaying a contingency table with heatmap(): svm.predPix.tabla svm.predPix CC DD LL NN NN2 CC 22 0 3 8 3 DD 0 27 0 1 0 LL 1 1 90 3 7 NN 2 0 1 11 4 NN2 0 0 5 1 20 heatmap(svm.predPix.tabla[5:1,], Rowv=NA,

[R] cairoDevice on Windows: succesful installation but does not load

Hi! I've tried to install playwith on Windows XP but, install.packages(cairoDevice) trying URL 'http://cran.at.r-project.org/bin/windows/contrib/2.13/cairoDevice_2.15.zip' Content type 'application/zip' length 52371 bytes (51 Kb) opened URL downloaded 51 Kb package 'cairoDevice' successfully

[R] Date_Time detected as Duplicated (but they are not!)

I have a problem with duplicated date_time stamps that I do not see as duplicated. I read a file with observations taken every 30 minutes: aur2009=read.csv(paste(datadir,AUR_ECPP_2009.csv,sep=/),sep=;,stringsAsFactors=F) aur2009[1:3,1:5] Date.Time E_filled E_filled_flag LE_filled

Re: [R] Date_Time detected as Duplicated (but they are not!)

timezone), but conversion by ‘as.POSIXct’) will do so. On Wed, May 18, 2011 at 3:53 PM, Agustin Lobo agustin.l...@ictja.csic.es wrote: I have a problem with duplicated date_time stamps that I do not see as duplicated. I read a file with observations taken every 30 minutes: aur2009

[R] Reshpaing for xyplot() or ggplot()

is at the melt() or cast() level. I've also tried reshape(). Any help, advice and/or pointer to an appropriate tutorial would be appreciated. Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax

Re: [R] Problem at subsetting matrix by using dimnames

), c(b, c, b)) or diag(`[`(mitest, i = c(a, b, c), j = c(b, c, b))) On Wed, Sep 2, 2009 at 6:57 AM, Agustin Lobo aloboa...@gmail.com mailto:aloboa...@gmail.com wrote: Given: mitest - matrix(1:16,ncol=4) dimnames(mitest)[[1]] - c(a,b,c,d) dimnames(mitest)[[2]] - c(a,b,c,d

[R] Problem at subsetting matrix by using dimnames

,c,b))] [1] NA NA NA NA NA NA despite the fact that mitest[a,b] [1] 5 works fine. How can I do the equivalent of mitest[cbind(c(1,2,3),c(2,3,2))] using row and column names? Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028

[R] Problem with area()

/eospansite/dummy/areaproblem.rda Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 e-mail agustin.l...@ija.csic.es http://www.ija.csic.es/gt/obster

[R] Function to eliminate blank space within strings?

Is there an specific function to eliminate blank space within strings? Thanks -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: agustin.l...@ija.csic.es http://www.ija.csic.es/gt/obster

[R] robustbase: cannot install

-lR /usr/bin/ld: cannot find -lRlapack collect2: ld returned 1 exit status Despite the fact that I do have Rlapack: ~$ R CMD config LAPACK_LIBS -L/usr/lib/R/lib -lRlapack Am I doing something wrong? Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole

[R] paste( / ) and paste( /' )

I wish to write using paste(), but paste(\,Hola,\,sep=) [1] \Hola\ while the same approach works with '' paste(\',Hola,\',sep=) [1] 'Hola' why this difference? how could I do it to get Hola ? Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis

[R] Solved: [Fwd: Matching failure in merge()]

:48 +0100 From: Agustin Lobo agustin.l...@ija.csic.es Reply-To: agustin.l...@ija.csic.es To: r-help@r-project.org Hi! I've done a merging between 2 dataframes using merge(): delme - merge(miDUNS50peqB,Bnomscodmunicipis,by.x=POBLACION,by.y=NOMMUNI,all.x=T,sort=F) After realizing some problems

[R] Matching failure in merge()

; And R is in English too) Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: agustin.l...@ija.csic.es http://www.ija.csic.es/gt/obster

[R] prcomp(X,center=F) ??

that the calculation should always be done using a centered data matrix. Finally, all the examples in the help page use centering (or scaling, which implies centering) Therefore, why the option center=F ? Agus Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis

[R] Equivalent to apply(x[,2:5],1,sum) for dataframe

What's the most correct way of doing the equivalent to apply(x[,2:5],1,sum) if x is dataframe in which the only numeric fields are in columns 2:5 ? (using apply returns a character vector) Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris

[R] Problem with data.frame()

(data.matrix(x2)) num [1:5284, 1:4] 79 215 84 58 127 134 15 122 101 99 ... - attr(*, dimnames)=List of 2 ..$ : chr [1:5284] 3 6 9 12 ... ..$ : chr [1:4] 0-4 5-9 10-14 15-19 Thanks for any clarification Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole

[R] Component merge in hclust objects

? Thanks Agus Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http://www.ija.csic.es/gt/obster __ R-help@r-project.org mailing

Re: [R] R on a eeepc 901 ?

the Xandros OS in the internal SSD. So I'd appreciate help to install and test R on Xandros to compare. Thanks Agus Agustin, which EEE PC 901 is it? The one for Windoze or Linux (Xandros)? el on 9/16/08 6:50 PM Agustin Lobo said the following: Hi! it seems that the directions in http://wiki.r

Re: [R] R on a eeepc 901 ?

-started:installation:eeepc refer to installing R on a eeepc 701 and actually Yes, it does :) fail for a 901 at the first step. What exactly failed? You are actually do not need the KDE desktop. Greetings Sigbert -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis

[R] R on a eeepc 901 ?

Hi! it seems that the directions in http://wiki.r-project.org/rwiki/doku.php?id=getting-started:installation:eeepc refer to installing R on a eeepc 701 and actually fail for a 901 at the first step. Does anyone have any experience on installing R on a eeepc 901? Thanks -- Dr. Agustin Lobo

Re: [R] Greyed text in the background of a plot

/ Regards, Yihui On Fri, Sep 12, 2008 at 8:35 PM, Agustin Lobo [EMAIL PROTECTED] wrote: Hi! Is there any way of having a greyed (ghosted) text (i.e, 2006) in the background of a plot? I'm making a dynamic plot and would like to show the year of each time step as a big greyed text in the background

[R] Greyed text in the background of a plot

=4237353244338529080sourceid=searchfeed ) Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http://www.ija.csic.es/gt/obster

[R] Converting from char to POSIX:

] [1] 2008-06-09 18:25:46 Where does the length 9 come from? How is it that I can get the 1065th element of delme3 if the length is just 9? How can I get a vector of POSIX elements with the same length as delme2 so that I can put it in a data.frame? Thanks! Agus -- Dr. Agustin Lobo Institut

Re: [R] Recoding

- Phil Spector Statistical Computing Facility Department of Statistics UC Berkeley [EMAIL PROTECTED] On Fri, 27 Jun 2008, Agustin Lobo wrote: Hi! Given a vector (or a factor within a df),i.e. v1 - c

[R] Recoding

origen - imagen imagen - aux } m - match(v, origen, 0) #print(match finished) if(resto==0) v - v*0 v[m 0] - imagen[m] v } Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028

[R] subscripting in data frames with NA

and levels(b2$p_desc_no3)[levels(b2$position_geo)==3] - NA which does not result into error but leaves b2$p_desc_no3 unchanged: table(b2$p_desc_no3) 0 1 2 42 44 66 what am i doing wrong? Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris

[R] reformatting R scripts for htmlize()

to write an script (could be done in R itself?) that would automatically modify the current script into a htmlize-friendly script? One problem is that while R accepts png() providing a default name for the png file, htmlize() fails unless a name is provided. Thanks! Agus -- Dr. Agustin Lobo Institut

[R] Subset in cast: compact form?

| variable== Ba,mean,na.rm=T) Is there a more compact form? something like select=c(ph, temperature, Ba) Thanks Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http

[R] Proportional axes

and ylim from 0 to 30, i.e., plot(runif(10,0,10),runif(10,10,30),xlim=c(0,30),ylim=c(0,30)) but this is a waste of space. Thanks -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL

Re: [R] Proportional axes

the information you want into a form you can incorporate elsewhere, Richard Agustin Lobo wrote: Hi! If I'm plotting 2 variables with the same units (i.e. wing length of 2 species), how can I force that X and Y axes have the same size for length unit? For example, if X ranges from 0 to 10 and Y ranges

Re: [R] Log or diary file

interface you will find for R. HTH Vincent Le lun. 19 mai à 03:24, Agustin Lobo a écrit : Hi! Is it possible to set a file to which both commands and output would get automatically saved? I've tried with sink(), but only get the output. I mean something like a combined history and sink, as you

Re: [R] Starting R from .RData in linux

. On 05/19/08 09:06, Agustin Lobo wrote: Hi! Is it possible to start R by clicking the .RData file in linux as in Windows? I've tried with ubuntu hardy using the right button and selecting R, but does not work. Is there any way to set it up? Marianne Promberger Wrote: You presumably need

Re: [R] Log or diary file

] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Agustin Lobo Sent: Monday, May 19, 2008 1:24 AM To: r-help@r-project.org Subject: [R] Log or diary file Hi! Is it possible to set a file to which both commands and output would get

[R] Startin R from .RData in linux

Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http://www.ija.csic.es/gt/obster __ R-help@r-project.org mailing list

[R] Log or diary file

, and you can state diary on and diary off to control what is being saved to the file. Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http

[R] grid() : start values?

another problem) would be making the axes coincident with the rectangular plot region, that is, if ploting plot(0:10,0:10) get (0,0) just on the bottom down corner of the rectangle. I've tried with several par() and could not find the way. Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de

[R] Summary of importing vector graphic formats written by R into Open Office

reported since so long ago that there is little hope that this problem will ever be solved within Open Office). Thank you all. Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL

[R] space between plots

I've set layout(matrix(c(1,3,2,4),ncol=2)) and par(mar=c(0, 4, 4, 0) + 0.0) and par(mar=c(0, 0, 4, 3) + 0.0) for the upper-left and upper-right plots, but still get an space between both plots, is there any way to completely eliminate that space? Thanks Agus -- Dr. Agustin Lobo Institut de

Re: [R] space between plots

this: layout(matrix(c(1,3,2,4),ncol=2)) par(mar=c(0, 4, 4, 0) + 0.0) plot(1) par(mar=c(0, 0, 4, 3) + 0.0) plot(2) plot(3) plot(4) the upper left and right are butted together. Exactly what 'space' are you talking about? On Fri, Apr 18, 2008 at 4:19 AM, Agustin Lobo [EMAIL PROTECTED

Re: [R] SVG format from any R graphic GUI?

Stefan, Stefan Grosse escribió: On Wed, 16 Apr 2008 18:51:30 +0200 Agustin Lobo wrote: AL 2. I need the figure displayed on the screen, I'm using AL Impress (the equivalent of ppt). Thus EPS is not an option, 1. So what please is wrong with png? I am using it with Word (If for some reasons

Re: [R] SVG format from any R graphic GUI?

Thanks, actually I had thought on Lynx as an option, so your opinion encourages me to try. But wanted to solve the problem within the OO environment as well. Agus Marc Schwartz escribió: See comments inline: Agustin Lobo wrote: Thanks for your answers, let me summarize the situation: 1

[R] bty

Dear list: Is there any way of getting the equivalent to what you get with bty=l for the right and bottom axes? and the equivalent to bty=7 for the upper and left axes? Thanks! Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona

[R] SVG format from any R graphic GUI?

! Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http://www.ija.csic.es/gt/obster __ R-help@r-project.org mailing list

Re: [R] SVG format from any R graphic GUI?

, and emf files are poorly imported by OO. Agus Stefan Grosse escribió: On Wed, 16 Apr 2008 13:20:35 +0200 Agustin Lobo wrote: ALI have problems with emf files because OpenOffice does ALa poor job at importing figures in this format. On the other ALhand, imported eps figures are not displayed

Re: [R] SVG format from any R graphic GUI?

seems to be to use EPS with a preview (which GSView can add for you). On Wed, 16 Apr 2008, Stefan Grosse wrote: On Wed, 16 Apr 2008 13:20:35 +0200 Agustin Lobo wrote: ALI have problems with emf files because OpenOffice does ALa poor job at importing figures in this format. On the other ALhand

[R] Subset: data frames and factor levels

? Thanks! Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: [EMAIL PROTECTED] http://www.ija.csic.es/gt/obster __ R-help@r-project.org

Re: [R] Subset: data frames and factor levels

Thanks, problem solved thanks to your hints, but, within subset(), drop=T still keeps the unused levels (I'm using 2.6.1 on win). Agus Stefan Grosse escribió: On Tuesday 01 April 2008 09:43:00 am Agustin Lobo wrote: AL I'm doing: AL etni - subset(etni, NAMECOM!=Maniquisito) AL AL where etni

Re: [R] barplot as Trellis graphic

to Trellis graphics? Agus Charilaos Skiadas escribió: On Mar 27, 2008, at 1:47 PM, Agustin Lobo wrote: Thanks, it was a matter of reshaping the data matrix as I usually have it, ie: datos - data.frame(x=abs(round(rnorm(100,10,5))),y=abs(round(rnorm(100,2,1))),f=factor(round(runif(100,1,3

Re: [R] barplot as Trellis graphic

[,3],datos[,3]))) and then require(lattice) barchart(V1~VAR|f,data=datos2) I get horizontal lines in the bars that I do not understand, though. Agus Deepayan Sarkar escribió: On 3/26/08, Agustin Lobo [EMAIL PROTECTED] wrote: Dear list, Is there any way of making barplots as a Trellis

[R] barplot as Trellis graphic (corrected)

In the code of my previous message the barplot should be: barplot(apply(sel[,-1],2,mean)) instead of barplot(sel) Sorry for the confusion. Agus Mensaje original Asunto: barplot as Trellis graphic Fecha: Wed, 26 Mar 2008 22:24:04 +0100 De: Agustin Lobo [EMAIL PROTECTED

Re: [R] Fwd: cut.dendrogram and cutree

be accessed in a binary tree fashion, use '1' for the tree on the left and '2' for the tree on the right. hope the above helps. yaomin -- Yaomin Xu Statstical Computing and Bioinformatics The Cleveland Clinic Foundation -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume

Re: [R] cut.dendrogram and cutree

, Agustin Lobo [EMAIL PROTECTED] wrote: Thanks, But this is not solving my problem, let me further explain by following the example in the help page for dendrogram: require(graphics); require(utils) hc - hclust(dist(USArrests), ave) dend1 - as.dendrogram(hc) dend2 - cut(dend1, h=70) plot(dend2

[R] cut.dendrogram and cutree

= 7.65, x.member = 20 ) correspond to the classes 1, 2, 3... in dendclass, where dendclas - cutree(hc,h=70) ? If they do not, how can I get a vector with the original individuals for each branch? Thanks Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra Jaume Almera (CSIC) LLuis Sole