similar to: Help with Lattice

Hi: Thanks a lot to Deepayan Sarkar, author of lattice I think, who solved my first query. I am afraid that I have another one. I am plotting several mutipanels boxplots (with one conditioning variable) on one page. The x, y and conditioning variable are all continuous variables. The x and conditioning variables are transformed to shingles before being plotted. The plot looks nice but there

Hi: I have made a lot of progress reading and manipulating large data files, thanks to the help of several of you. I am now stuck with writing the final file with the following error (see also the full transcript below): Error in as.matrix.data.frame(x) : length of dimnames[2] not equal to array extent I do not know what it means and could not find how to get around that in the manual nor

Hi: I am a very recent user of R 1.40 for MacOSX and went through the help files but did not not find the function that I am looking for. I have two comma separated files which I want to partly combine. File 1 has 49670 lines and 26 columns (separated by commas) and looks like this: ...

Hi: I have have taken over from a colleague who prepared an R package and failed to build it on Windows. I am doing this with unix as I am a mac user. Below is the output I get when I use the build command: [gattuso:unix/R/CO2.Rcheck] gattuso% R CMD build CO2 * checking for file 'CO2/DESCRIPTION' ... OK * preparing 'CO2': * checking whether 'INDEX' is up-to-date ...

Hi: I am posting this message for a colleague who has a lot of trouble to build an R package on Windows. He did not find a solution to his problems on the R-help archives and hopes that one of the R gurus will be able to help. He has a directory "CO2" which should contain all the required files and directories: DATA: DESCRIPTION MAN: R: README SRC: The following command is

Hi: I am using the following command: xyplot(dat6$CO3*1e6 ~ dat6$irradiance, data=dat6, group=ref, xlab=list(label=expression(paste("Irradiance (", mu, "mol photons", m^"-2", " ", s^"-1", ")")), cex=1.3), ylab=list(label=expression(paste("Carbonate concentration (x ", 10^"6", " ", kg^"-1",

> Dear all, > I am looking for basic examples of the use of tcltk in R. > Thank you in advance for any information, R-script, ... Most of us started from the handful of examples you can see via demo("tcltk") There are also some good examples in the mailing list archives though some of the R to Tcl/Tk interfaces changed slightly in recent releases. -- According to the

Hello, I do not understand the following behaviour. Could someone explain me what happens? > a <- NULL > a$item <- 1:3 > a$item [1] 1 2 3 > rm(a) > a <- NULL > a[["item"]] <- 1:3 Error: more elements supplied than there are to replace Why do I get an error message using list[["item"]], and not using list$item? Best, Philippe Grosjean

Hi, This is (presumably?) a bug in R 1.7.0 under Windows. I have not tested it on other systems. Attachment of packages is case sensitive but not library(), resulting in multiple loadings of the same package if the library name is spelled differently. The following example loads the `tools' package, once as `tools' and once as `Tools'. This behavior is the same with all packages and

Sorry for this off-topic subject. I am fighting for running SJava under Windows. SJava_0.64 (compiled by Simon Urbanek, thanks), R 1.6.2, Java JDK 1.4.0_02, Windows XP pro: > library(SJava) > .JavaInit() Error in .JavaInit() : Couldn't start Java Virtual Machine: Cannot find the Omegahat interface manager class. Check you classpath! > # And the second time... > .JavaInit() It

Dear R-listers: I am trying to work with a big (262 Mb) file but apparently reach a memory limit using R on a MacOSX as well as on a unix machine. This is the script: > type=list(a=0,b=0,c=0) > tmp <- scan(file="coastal_gebco_sandS_blend.txt", what=type, sep="\t", quote="\"", dec=".", skip=1, na.strings="-99", nmax=13669628)

Hi, I was able to connect to a MySQL database (called "zooscan" and with a table "serie" under Win XP with R 1.6.2 using: > library(DBI) > library(RMySQL) Warning message: DLL attempted to change FPU control word from 8001f to 9001f > drv <- dbDriver("MySQL") > con <- dbConnect(drv, dbname="zooscan") > dbExistsTable(con,

Hi, I run R 2.5.1 on Mac OS 10.4.10 and fail to build a package manual. Below is the transcript. The problem is a "pdflatex: command not found" error but I think that I have a fully working latex install (installed with iInstaller) as I use TeXShop with no problem. Could someone help? Jean-Pierre -----------------------------------------------------------------------------

OK, I have made a little search about this "problem" that apparently occurs only on Windows platform... (but I am sure most of you are already aware of it): the slow down is due to the adoption of a different algorithm for pow in mingw 3.x. This is motivated by some other changes in mingw. Here is a quote of Danny Smith that did this change: >When mingw changed default FPU settings

I do not understand what happens here (under Win XP): a <- abs(matrix(rnorm(800*800)/2, ncol=800, nrow=800)) system.time(b <- a^1000)[3] took about 1 sec on my computer with R 1.7.0 and it takes now 4.59 sec with R 1.7.1 Similarly, phi <- 1.6180339887498949 a <- floor(runif(750000)*1000) system.time(b <- (phi^a - (-phi)^(-a))/sqrt(5))[3] took about 0.9 sec with R 1.7.0, and it

I do not understand what happens here (under Win XP): a <- abs(matrix(rnorm(800*800)/2, ncol=800, nrow=800)) system.time(b <- a^1000)[3] took about 1 sec on my computer with R 1.7.0 and it takes now 4.59 sec with R 1.7.1 Similarly, phi <- 1.6180339887498949 a <- floor(runif(750000)*1000) system.time(b <- (phi^a - (-phi)^(-a))/sqrt(5))[3] took about 0.9 sec with R 1.7.0, and it

Happy New Year, r-helpers! I am using lowess to smooth a scatter plot, xx<-lowess(xinput,f=.04) #defaults for other args followed by turnpoints(xx$y) #defaults for other args I plot the smoothed result as well as turnpoints (using yy$tppos) on top of raw data plot. Result is exactly as expected, graphically. For another purpose, I calcuate the difference between turnpoints (representing

Although the 'spatial' documentation doesn't mention that extractAIC works, it does seem to give an output. I may have misunderstood, but shouldn't the following give at least the same d.f.? > library(spatial) > data(topo, package="MASS") > extractAIC(surf.ls(2, topo)) [1] 46.0000 437.5059 > extractAIC(lm(z ~ x+I(x^2)+y+I(y^2)+x:y, topo)) [1]

Hello, This message is little off-topic in R-help. Sorry for that, but not all interested people are wired yet to r-sig-gui (http://www.stat.math.ethz.ch/mailman/listinfo/r-sig-gui). Thanks for your comprehension. A preview version of SciViews (a Graphical User Interface for R under Windows, http://www.sciviews.org) was released a few weeks ago. Since then, the Web site recorded several

Hello, I suppose this is a general behavior with external function calls, so I do not post (yet) a specific bug report. Could someone explain this? a <- rep(1, 20) + rnorm(20, mean=0.00001, sd=0.0001) b <- embed(a, 3) # I want to know where the item in column 2 is greated than both col 1 and 3 (peak) test1 <- max.col(b) == 2 # ... or I could use a less optimal code test2 <- apply(b,