similar to: incorrect number of subscriptions error in user-written function

Dear R-help, I wonder if you could give me some suggestions in how to do a union join of two data frames as follow: -> union join the common column, and insert a 0 if one is missing. I made a function to perform the following, and I know it may not that quite welly written, but it works. Any suggestions are welcome, many thanks. Anthony > q1 =

Hi everybody, I have problem with samba and really don't know how to solve it. Seems as if samba ignored write list and read list options in smb.conf. Because when I have writeable = yes read list = dummy user dummy can still write/delete. When I try to do it opposite, i.e. read only = yes write list = dummy then dummy cannot write. The point is, I have bunch of users specified by valid

Has the algorithm used by batchSOM for hexagonal topology changed in some way between R 1.9.1 and 2.0.1? In comparing between the two versions, I get identical results for rectangular topology, but very different results for hexagonal topology. (I tried to find release notes for the VR bundle that might explain the changes, but was unable to locate these anywhere.) Here's the script I

I can do this by writing (and reading) the file according to some format of my own devising, but I'm wondering if there is a built-in way to write and then restore a matrix with not only the dimnames (which write.table/read.table can preserve) but also the names(dimnames)? Example: > M <- matrix(1:4, 2, 2) > dimnames(M) <- list(xdim=c("a", "b"),

Hello, I have a problem writing a variable to an existing file. Below is a part of my script and how it fails. I can't find "create.var.ncdf" in help Thanks for any help. Mark nc <- open.ncdf(ncname, readunlim=FALSE, write=TRUE ) missing <- 1.e+30 xdim <- nc$dim[["west_east"]] ydim <- nc$dim[["south_north"]] tdim <- nc$dim[["Time"]]

Good day everyone. I have a question concerning *ggobi* and *rggobi*package. I am using R version 2.13.2 and Ubuntu 11.04. I am new to ggobi and I'm reading the book *Interactive and Dynamic Graphics for Data Analysis *. The book was published 2007. My question is: how can I export a graph for *rggobi* to pdf to include it in latex. For example: I have lots of variables in ggobi and I need to

Good day everyone. I have a question concerning *ggobi* and *rggobi*package. I am using R version 2.13.2 and Ubuntu 11.04. I am new to ggobi and I'm reading the book *Interactive and Dynamic Graphics for Data Analysis *. The book was published 2007. My question is: how can I export a graph for *rggobi* to pdf to include it in latex. For example: I have lots of variables in ggobi and I need to

Hi all, I originally posted this to the bioconductor group, but maybe it's better suited to the r-help... I'm using som() to partition samples of gene expression data into clusters. The point is to classify control vs. experimental cases (sample clustering). The original matrix was 22283 x 8. The 8 samples have 4 controls and 4 experimentals. I transposed the matrix so that its dim

Dear All, Would you please tell me how to display the sample No. on the map ? ---Below commands don't display the sample No.(from 1 to 150).--- library(som) library(klaR) iris.som3 <- som(iris[,1:4], xdim = 14,ydim = 6) library(klaR); opar<- par(xpd = NA) shardsplot(iris.som3, data.or = iris,label = TRUE) legend(3.5,14.3, col = rainbow(3), xjust =0.5, yjust = 0,legend =

Hello R users, I'm using SOM() to cluster a gene expression data set the syntax i used was dataGrid <- c(somgrid(xdim = 3, ydim = 3, topo = c("rectangular","hexagonal"))) dataClusters <- SOM(dataMatrix, grid = dataGrid) plot(dataClusters) it seems that this works just fine but the thing i can't figure out is how to determine where each data point has been

On Mon, Aug 14, 2017 at 5:29 AM, <raphael.felber at agroscope.admin.ch> wrote: Dear all > > I'm a newbie regarding netcdf data. Today I realized that I maybe do not > understand some basics of the netcdf. I want to create a *.nc file > containing three variables for Switzerland. All data outside of the country > are NAs. The third variable is calculated from the first two

Dear all I'm a newbie regarding netcdf data. Today I realized that I maybe do not understand some basics of the netcdf. I want to create a *.nc file containing three variables for Switzerland. All data outside of the country are NAs. The third variable is calculated from the first two variables. Basically there is no problem to do that. I copy the file with the data of the first variable,

I am trying to understand how the SOM algorithm works using library(class) SOM function. I have a 1000*10 matrix and I want to be able to summarize the different types of 10-element vectors. In my real world case it is likely that most of the 1000 values are of one kind the rest of other (this is an oversimplification). Say for example: InputA<-matrix(cos(1:10),nrow=900,ncol=10,byrow=TRUE)

biplot could benefit from allowing colour palettes, eg. allowing col=list(xcol= xxpal, ycol= yypal) where xxpal and yypal are palettes. For an example, see the robust pricipal components in the principal componetent section of <http://www.statlab.uni-heidelberg.de/data/olive/olive10.pdf> The required change is to replace col[1L] resp. col[2L] by col[1L]] reep. col[[2L]] for example in

Hi List, I am trying to plot a grid with an overlayed height. I have a dataframe with four variables: x,y,gridvalue,height. The dataframe has 2.5mio observations (ie grid points), I assign colors through the gridvalue using map_color_gradient thus producing: x,y,gridvalue,height,gridcol as variables of the dataframe. The grid dimensions are 1253 x 2001 (=2507253 data points). My attempts with

Hi, Excuse me for asking this silly question. But I really couldn't understand why cov() and ccov() don't work for my calculation of covariance matrix. a <- matrix(1:8, 2, 4) a [,1] [,2] [,3] [,4] [1,] 1 3 5 7 [2,] 2 4 6 8 > ccov(a) Error in solve.default(cov, ...) : Lapack routine dgesv: system is exactly singular I also tried colume bind, but it

Hi, Given > sec98 <- factor(rep(1:2,3), labels=c("A", "B")) > sec99 <- factor(rep(2:1,3), labels=c("A", "B")) > sec99[c(2,5)] <- NA > sec00 <- factor( c( rep(1,3), rep(2,3) ), labels=c("A", "B")) > sec00[c(2,4)] <- NA > sec1 <- ifelse(!is.na(sec99), sec99, ifelse(!is.na(sec00), sec00, NA

Can anyone help me resolve this? A part of the R function looks like this: print(is.loaded('merge_xtabs_patterns_file')) print(is.loaded('merge_xtabs_patterns_file_')) .Fortran('merge_xtabs_patterns_file_',ydim[1],ydim[2],x=as.integer(as.matrix(y)),na=as.integer(c), maxD=as.integer(maxD),lrowmem=length(rowmem),rowmem=as.integer(rowmem),

I have some elements that could show up on the same page multiple times. Can anyone recommend a good way to do a replace_html in this case? I may have, for example, 3 posts that have been repeated on a page, all with the same id (e.g, <div id="post-11">) and I want to perform a replace_html against all of them. Doing it the normal way just replaces the first. Can someone help?

Dear all, I'm used to draw bubble plots (gstat) to get a first view of my spatial marked data, but I can't find a way to label legend with marks after having replaced them by a numeric scale. I have a data frame with numeric coords and factors: > names(data) [1] "x" "y" "bloc" "sub" "inoc"