[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

9902, Adjusted R-squared: 0.9902 F-statistic: 1.009e+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 re

[R] Friedman test

76,0.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? Than

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

d not attached): [1] colorspace_1.2-4 dichromat_2.0-0digest_0.6.4 grid_3.0.2 gtable_0.1.2 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.co

[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

x86_64-apple-darwin10.8.0 (64-bit) locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] lubridate_1.3.0 loaded via a namespace (and not attached):

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

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

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

ets methods base other attached packages: [1] sp_0.9-99 loaded via a namespace (and not attached): [1] grid_2.15.1lattice_0.20-6 tools_2.15.1 -- Dr. Agustin Lobo Institut de Ciencies de la Terra "Jaume Almera" (CSIC) Lluis Sole Sabaris s/n 08

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 wrote: > On 09/10/2012 09:54, Agustin Lobo wrote: >> >> If I do: >>> >>> a >> >> [1]

[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 o

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

rmutation of the characters "d", "m", or "y" delimited by a separator (possibly null), e.g., "m/d/y", "d-m-y", "ymd", are all valid; " Thanks! Agus On Fri, Jun 8, 2012 at 12:54 PM, Gabor Grothendieck wrote: > On Fri, Jun 8, 20

[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-

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

Month" "Day" "Hour" "doy2" ... Index: POSIXct[1:12001], format: "2007-03-04 00:00:00" "2007-03-04 00:30:00" "2007-03-04 01:00:00" "2007-03-04 01:30:00" "2007-03-04 02:00:00" "2007-03-04 02:30:00"

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

) : Incompatible methods ("Ops.POSIXt", "Ops.Date") for "<=" > str(delme) ‘zoo’ series (without observations) > delme <- window(alyL32007z, start = c(2007, 3), end = c(2007, 12)) > str(delme) ‘zoo’ series (without observations) Is this a problem with t

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 wrote: > On Mon, May 28, 2012 at 5:35 AM, Agustin Lobo > wrote: >> I'm do

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

dow(...) : invalid 'ylim' value Any help 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. 34 934110012 e-mail agustin.l...@ictja.csic.es https://sites.google.com

[R] Problem with zoo::window()

Error in plot.window(xlim = xlim, ylim = ylim, log = log, yaxs = pars$yaxs) : need finite 'ylim' values Calls: plot ... boxplot -> boxplot.default -> do.call -> bxp -> plot.window why do I get the warning in window()? My command follows what is mentioned in the manual:

[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

[R] as.numeric() and POSIXct format

which the inverse 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 F

[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, Colv=NA,c

[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

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

rptime’ does not validate such times (it does not assume a >      specific timezone), but conversion by ‘as.POSIXct’) will do so. > > > > On Wed, May 18, 2011 at 3:53 PM, Agustin Lobo > wrote: >> >> I have a problem with duplicated date_time stamps that I do not

[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

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

B 5 5 142 95 112 176 161.5 133 0.42 B 6 6 89 113 112 125 128.0 195 0.42 B ... so that I could use xyplot(RefVal ~ XR | Band, data=datosX) but do not get to that. Don't even know if the problem is at the melt() or cast() level. I've also tried reshape(). Any hel

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

lazuanna wrote: Try this: mapply(function(x, y)mitest[x, y], c("a", "b", "c"), c("b", "c", "b")) or diag(`[`(mitest, i = c("a", "b", "c"), j = c("b", "c", "b"))) O

[R] Problem at subsetting matrix by using dimnames

8 12 16 I can do: > mitest[cbind(c(1,2,3),c(2,3,2))] [1] 5 10 7 but using the names does not work: > mitest[cbind(c("a","b","c"),c("b","c","b"))] [1] NA NA NA NA NA NA despite the fact that > mitest["a","b"]

[R] Problem with area()

Could this be a problem with an inappropriate use of approxfun() ? Data can be found here: http://sites.google.com/site/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

[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.

[R] robustbase: cannot install

linux-gnu/4.2 -lgfortran -lm -L/usr/lib/R/lib -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 Ciencie

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

I wish to write "" using paste(), but > paste("\"","Hola","\"",sep="") [1] "\"Hola\"" > while the same approach works with '' > paste("\'","Hola","\'",

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

Thu, 19 Mar 2009 20:04:48 +0100 From: Agustin Lobo 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)

[R] Matching failure in merge()

ly appreciated, I'm really astonished. I think it might involve an encoding problem with the non-ascii characters, but don't get to see it. I'm using R 2.8.1 on ubuntu 8.04 (in english; And R is in English too) Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terr

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

(centered and possibly scaled) data matrix..." The parenthesis implies some ambiguity, but I do interpret the sentence as indicating 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 implie

[R] reading files (readOGR) and non-ascii characters

2:34 +0100 From: Agustin Lobo <[EMAIL PROTECTED]> Reply-To: [EMAIL PROTECTED] To: sig-geo <[EMAIL PROTECTED]> Hi! I have problems at reading a shp file with readOGR(), names with non-ascii characters (accents ',`) in the dbf table are imported as i.e. Barber\x85 del Vall\x8as

[R] Problem with data.frame()

79 215 84 58 127 134 15 122 101 99 ... $ 5-9 : num 76 180 32 56 81 106 10 112 128 96 ... $ 10-14: num 68 145 39 46 78 81 8 92 142 107 ... $ 15-19: num 73 170 49 52 103 77 10 116 129 129 ... > str(data.matrix(x2)) num [1:5284, 1:4] 79 215 84 58 127 134 15 122 101 99 ... - attr(*, "

[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" (C

[R] Component merge in hclust objects

elements 6,10 and 3 Is this right? Am I understanding correctly? 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 ___

Re: [R] R on a eeepc 901 ?

rwiki/doku.php?id=getting-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 T

Re: [R] R on a eeepc 901 ?

I'll be erasing 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 t

[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

" package: http://cran.r-project.org/web/packages/animation/index.html http://animation.yihui.name/ 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&#

[R] Greyed text in the background of a plot

videoplay?docid=4237353244338529080&sourceid=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] htt

[R] Converting from char to POSIX:

h the same length as delme2 so that I can put it in a data.frame? 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

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 d

[R] Recoding

aux <- origen 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

[R] subscripting in data frames with NA

: Factor w/ 4 levels "0","1","2","3": 4 1 4 4 3 NA 3 3 3 4 ... > table(b2$p_desc_no3) 0 1 2 42 44 66 and > levels(b2$p_desc_no3)[levels(b2$position_geo)=="3"] <- NA which does not result into error but leaves b2$p_desc_no

[R] reformatting R scripts for htmlize()

() Any ideas on how 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. Than

[R] Subset in cast: compact form?

="ph" variable=="temperature"| 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.

Re: [R] Proportional axes

liver 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 1

[R] Proportional axes

m 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 93411001

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

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

ore same result. 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: Yo

Re: [R] Log or diary file

Windows) user 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

[R] Log or diary file

atlab, 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

[R] Startin R from .RData in linux

hanks 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.o

[R] Problem with help system on linux

ate your help. 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-

[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 Cienci

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

problems at importing emf files by OpenOffice have been 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 080

Re: [R] space between plots

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 >&

[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 Instit

[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

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 i

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 summariz

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

;> I've seen so many problems with SVG renderers (especially related to >> fonts) that I would prefer not to rely on it. >> >> If this is Windows, the preferred route seems to be to use EPS with a >> preview (which GSView can add for you). >> >> On We

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

e emf as vector option, and emf files are poorly imported by OO. Agus Stefan Grosse escribió: > On Wed, 16 Apr 2008 13:20:35 +0200 Agustin Lobo wrote: > AL>I have problems with emf files because OpenOffice does > AL>a poor job at importing figures in this format. On the other > AL&g

[R] SVG format from any R graphic GUI?

hanks! 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.o

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, NAMECO

[R] Subset: data frames and factor levels

e care of recalculating the factor levels automatically? 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

Re: [R] barplot as Trellis graphic

guide 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

Re: [R] barplot as Trellis graphic

00),rep("y",100)),f=factor(c(datos[,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: >> De

[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 PROTEC

[R] barplot as Trellis graphic

umn of mydata would be the class. But I'm looking for a nicer code and a nicer graphic using the approach of the lattice package. Thanks! -- Dr. Agustin Lobo Institut de Ciencies de la Terra "Jaume Almera" (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fa

Re: [R] cut.dendrogram and cutree

;Arkansas" "Tennessee" "Georgia" > "Colorado" "Texas" >> labels(dend2$lower[[4]]) > [1] "Idaho" "Nebraska" "Kentucky" "Montana" > "Ohio" "

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

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 Alm

[R] cut.dendrogram and cutree

ch 3" (h= 44.8 midpoint = 6.8, x.member = 14 ) `--leaf "Branch 4" (h= 54.7 midpoint = 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 eac