search for: grum

...bset(x, var1!="NA") and repeat for each variable. It would be nice to be able to do it in one fell swoop. Also, surprisingly, it doesn't always work. Sometimes I'm left with an empty set even though not all rows have missing values for the variable. Cheers, mikkel Mikkel Grum, PhD Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** c/o ICRAF PO Box 30677 Nairobi, Kenya Tel: 254 2 524505 / 524500 Fax: 254 2 524501 / 524001 m.grum at cgiar.org ipgri-kenya at cgiar.org www.ipgri.cgiar.org -.-.-.-.-.-.-.-.-.-.-.-...

Is there a way of extracting mean distance or dissimilarity for a given strata from a 'dist' or 'dissimilarity' object, e.g. extract mean distances for each species in Anderson's iris data? data(iris) iris.dist<-dist(iris[,1:4]) then what? Mikkel Grum, PhD Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** PO Box 30677 Nairobi, Kenya Tel: +254 2 524505(direct)/524500(IPGRI) Fax: +254 2 524501(IPGRI)/524001(ICRAF) m.grum at cgiar.org www.ipgri.org

...:4[codes(xx$TYPBOTA)],pch=c(4,3,1,2)[codes(xx$TYPBOTA)]) grid() I get the following error message Error in 1:4[codes(xx$TYPBOTA)] : NA/NaN argument In addition: Warning message: Numerical expression has 162 elements: only the first used in: 1:4[codes(xx78$TYPBOTA)] Best wishes, Mikkel Mikkel Grum, PhD Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** c/o ICRAF PO Box 30677 Nairobi, Kenya Tel: 254 2 524505 / 524500 Fax: 254 2 524501 / 524001 IVDN Tel: 645 4505 / 4500 IVDN Fax: 645 4501 / 4001 m.grum at cgiar.org ipgri-kenya at cg...

...th 03, ie. 2May03) and spaced with the observations (every seventh day) rather than the 10 days that seem to be default?? > Week [1] 2May3 9May3 16May3 23May3 30May3 7Jun3 14Jun3 ALSO, why won't the three last lines of code in the top example show on the graph? cheers, Mikkel Mikkel Grum International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** PO Box 30677 00100 Nairobi, Kenya Tel: +254 20 524505(direct)/524500(IPGRI) Fax: +254 20 524501(IPGRI)/524001(ICRAF) m.grum at cgiar.org www.ipgri.org

Is there any function in R like is.not.found(x, y) meaning if you can't find object x, then use object y?? Mikkel Grum

Hi, all: I'm using centos 6.2-x86_64, the built-in kernel version is 2.6.32-220.el6.x86_64. Now I compiled kernel of version of linux-2.6.32-220.23.1.el6.x86_64, which get from kernel-2.6.32-220.23.1.el6.src.rpm. I have put vmlinuz and initramfs into /boot directory, as follow: [kdev at panda boot]$ ll /boot/vmlinuz-2.6.32-* ./initramfs-2.6.32-* -rwxr-xr-x. 1 root

...work fine xx<-subset(x,LOCAL.NAME==c("Chikota"),select=c(13,16,19,23,26,29,30,33:48)) xx<-subset(x,LOCAL.NAME==c("Chifumbata"),select=c(13,16,19,23,26,29,30,33:48 )) but I want a file with both the Chifumbatas and the Chikotas and several others. Cheers, Mikkel Mikkel Grum Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** m.grum at cgiar.org ipgri-kenya at cgiar.org www.ipgri.cgiar.org -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.c...

...4])*(log((xtab[,14]+0.000001)/s um(xtab[,14])))) xshw<-rbind(xshw1,xshw2,xshw3,xshw4,xshw5,xshw6,xshw7,xshw8,xshw9,xshw10,xsh w11,xshw12,xshw13,xshw14) xtap<-cbind(xtap1,xtap2,xtap3,xtap4,xtap5,(xtapx<-(xtap1+xtap2+xtap3+xtap4+x tap5)),xshw) Thanks a lot in advance!! Mikkel Mikkel Grum, PhD Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** c/o ICRAF PO Box 30677 Nairobi, Kenya m.grum at cgiar.org ipgri-kenya at cgiar.org www.ipgri.cgiar.org -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...me(cbind(f, g)) i <- 2 df[df$b %in% df2[i, 2], ] # or df[df$b %in% levels(df2)[as.integer(df2[i, 2])], ] # this is the closest I've got The results I want is given by > df[df$b %in% c(7:9, 12:14), ] a b 7 A7 7 8 A8 8 9 A9 9 12 B2 12 13 B3 13 14 B4 14 Can it be done? Mikkel Grum

...t;2017-11-13" [1] "Date" [1] 17484 [1] "numeric" [1] "2017-11-14" [1] "Date" Why is this apparently not a bug? Are there other types that change type when looped over? Kind regards Mikkel _____________________________________________________ *Mikkel Grum* +44 7377337321 (mobile) mikkelgrum (Skype) [[alternative HTML version deleted]]

...0, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: RSQLite DBI "0.4-1" "0.1-10" Mikkel Grum Genetic Diversity International Plant Genetic Resources Institute

...0, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: RSQLite DBI "0.4-1" "0.1-10" Mikkel Grum Genetic Diversity International Plant Genetic Resources Institute

...all the analyses. How do I select just those that I want for the next analysis? I've worked through R for beginners, the import/export manual, an introduction to R, and considerable chunks of the reference manual, as well as a few other texts and I could do with some help. -mikkel Mikkel Grum Genetic Diversity Scientist International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group m.grum at cgiar.org ipgri-kenya at cgiar.org www.ipgri.cgiar.org -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.t...

I run a batch file with the following command in Windows XP: C:\R\R-2.12.1\bin\Rterm.exe --no-save --no-restore <C:\users\me\file.R> C:\users\me\file.out 2>&1 Is there any way to get only the output of R in file.out, without getting all the code from file.R too? Any help greatly appreciated, Mikkel

I have a situation where I'm filling out a dataframe from a database. Sometimes the database query doesn't get anything, so I end up trying to place NULL in the dataframe like below. > temp <- NULL > xmat <- as.data.frame(matrix(NA, 2, 3)) > xmat[1, 2:3] <- temp Error in if (m < n * p && (n * p)%%m) stop(gettextf("replacement has %d items, need %d",

Is the following expected behaviour for a date used in an ifelse function? > date <- Sys.Date() > date [1] "2007-12-30" > ifelse(TRUE, date-1, date) [1] 13876 > ifelse(FALSE, date-1, date) [1] 13877 > ifelse(TRUE, as.character(date-1), date) [1] "2007-12-29" > if (TRUE) {date} [1] "2007-12-30" It would seem more natural to me if a date produced

...agnesmol,cex=0.6) coagnes<-cophenetic(agnesmol) cor(coagnes,daisemol) cohclust<-cophenetic(hclustmol) cor(cohclust,daisemol) I'm using R 1.7.1 (2003-06-16) on Windows XP. Library "cluster" version 1.7.4 and library "mva" version 1.7.1. Best wishes, Mikkel Mikkel Grum International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** PO Box 30677 00100 Nairobi, Kenya Tel: +254 20 524505(direct)/524500(IPGRI) Fax: +254 20 524501(IPGRI)/524001(ICRAF) m.grum@cgiar.org www.ipgri.org ------_=_NextPart_000_01C35A53.75780090 Content-Type: text/plain...

Dear useRs, Why are the rotated blue and yellow boxes in the example below clipped outside of 6 x 6 inch window in the middle of the page?? Where does the 6 x 6 inch window come from? I would like to make use of the entire page. > library(grid) > pdf(file = "FarmMaps.pdf", paper = "a4") > pushViewport(viewport( + width = unit(7.6, "inches"), height =

...) library(maptree) data(flower) dfl2<-daisy(flower,type=list(asymm=c(1,3),ordratio=7)) hdfl2<-hclust(dfl2) prune.dfl2<-prune.clust(hdfl2,k=6) plot(prune.dfl2) Of course the following will display all six clusters, but without a scale: draw.tree(prune.dfl2) Best wishes, Mikkel Mikkel Grum International Plant Genetic Resources Institute (IPGRI) Sub-Saharan Africa Group *** PO Box 30677 00100 Nairobi, Kenya Tel: +254 20 524505(direct)/524500(IPGRI) Fax: +254 20 524501(IPGRI)/524001(ICRAF) m.grum at cgiar.org www.ipgri.org

Dear R-helpers, When I use sqlQuery in the package RODBC, I cannot break the line, but have to write the entire SQL Query on the same line. Is this expected behaviour? It is definitely workable, but makes the queries a slightly difficult to read and edit. I'm using R 1.9.1 and RODBC 1.0-4 on Windows Server 2003 and querying a Sybase database. Best wishes, Mikkel