similar to: se.contrast ....too hard??? .... Too easy????? .....too trivial???? ...... Too boring.....too????????

Hi all, Just to follow up Don Driscoll's earlier post, can anyone please explain why "se.contrast" fails here?? > shp<-factor(rep(c("reserve","strip"),each=96)) > > > site<-factor(rep(c("1g","1p","1t","2g","2p","2t","3g","3p","3t","4g","4p

G'day, I'm fitting a simple one-way nested anova and would like to obtain standard errors or confidence limits. I'm using se.contrast to try to get se for the contrast between the two levels of the main effect, but I get an error message (see below). What is going on? Don > shp<-factor(rep(c("reserve","strip"),each=96)) >

G'day, I'm fitting a simple one-way nested anova and would like to obtain standard errors or confidence limits. I'm using se.contrast to try to get se for the contrast between the two levels of the main effect, but I get an error message (see below). What is going on? Don > shp<-factor(rep(c("reserve","strip"),each=96)) >

With this example set.seed(123) A <- matrix(runif(40), nrow = 8) y <- 1:nrow(A) A.laqr <- qr(A, LAPACK=TRUE) both qr.qy(A.laqr,y) and qr.qty(A.laqr,y) give the respective error messages Error in qr.qy(A.laqr, y) : 'b' must be a numeric matrix Error in qr.qty(A.laqr, y) : 'b' must be a numeric matrix However when Lapack is not used as in A.liqr <- qr(A,

The man page for se.contrast, when describing the contrast.obj parameter, states that "Multiple contrasts should be specified by a matrix as returned by contrasts." When doing an unbalanced single factor ANOVA, using a contrast.obj as returned by contrasts results in the following error from qr.qty when se.contrast is called: Error in qr.qty(object$qr, contrast) : qr and y must have

Hii All, I'm really a linux newbie, I managed to run into a problem with rsync. I want to conduct backups of various files and directories over about miles of Internet, so I chose ssh and rsync for the job. Earlier it used to work fine, but suddenly donno what happened it started to show the following errors sync: mkstemp

Dear All, I have a pair of .cpp and .def file can be compiled using VC++ and works perfectly well in S-PLUS. I wanted to do the same for R; so I followed the guidline given in "Building from a source-code library under Windows" as much as possible and manage to compile them using VC++ and call it from R. But it gives different answer from the one called from S-Plus. I know that I did

The replies so far have helped me see the issues more clearly. Further comments: 1. This issue started with a bug report from a user: library(survival) fform <- as.formula(Surv(time, status) ~ age) myfun <- function(dform, ddata) { predict(coxph(dform, data=ddata), newdata=ddata) } Gabor's suggestion to change the call is a useful idea but not completely relevant: I'm

Given the data structure below and the call to ggplot2, how can I increase the size of the axis scale points, the line weight, and the size of the legend? ddata <-structure(list(year = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("2003", "2007"), class = "factor"), area = structure(c(7L, 6L, 1L, 2L, 3L, 4L, 5L, 7L, 6L, 1L,

Dear R users, I use R 2.4.0 on an iMac running Mac OS X 10.4.8, with a 2.16GHz Intel Core 2 Duo and 2GB 667 MHz DDR2 SDRAM. > sessionInfo() R version 2.4.0 (2006-10-03) i386-apple-darwin8.8.1 locale: es_ES.UTF-8/es_ES.UTF-8/es_ES.UTF-8/C/es_ES.UTF-8/es_ES.UTF-8 attached base packages: [1] "methods" "stats" "graphics" "grDevices"

Hello, I'm having problems reading a shapefile with read.shape (maptools). I'm absolutely sure my file is there, but I get "no such file". The wd is ok, since read.table for example does find the file. > getwd() [1] "D:/somedirectory/R scripts" > read.table("cities.shp") Error in read.table("cities.shp") : empty beginning of file In addition:

The following snippet suggests that there is either a bug in qr(,LAPACK=T), or some bug in my understanding. Note that the detected rank is correct (= 2) using the default LINPACK qr, but incorrect (=3) using LAPACK. This is running on Linux Redhat 9.0, using the lapack library that comes with the Redhat distribution. I'm running R 1.7.1 compiled from the source. If the bug is in my

Hi, I am new to R and am a very basic user. I'm importing SHP files, adding plots of random locations within my polygon (these files have GPS data), and then want to save these plots (intact with added points) as KML files to look at in GoogleEarth (or possibly as SHP files which I can then convert into KML). This is as far as I've gotten (code below) and am lost as to next steps, Any

Hi, I have a data frame listing US counties and a quantity ("number") per county and I have a shapefile of the US with county ID's. I would like to plot the "number" variable on a map (in the shapefile) using a color range per county (e.g. white = min(number) = 2, black = max(number) = 15). Can anyone help me actually plotting the data on the map? This is how far I got.

Dear list, I have found a thread dealing with similar problems with a Mac, but somehow I cannot get it to work and the problem is slightly different: I am using the shell() command to execute a program from within R. The problem is that I need to execute SDKShell.bat file first that sets the environment variables for the program ogr2ogr.exe: shell(cmd =

I want to perform linear regression on groups of consecutive rows--say 5 to 10 such--of two matrices. There are many such potential groups because the matrices have thousands of rows. The matrices are both of the form: > shp[1:5,16:20] SL495B SL004C SL005C SL005A SL017A -2649 1.06 0.56 NA NA NA -2648 0.97 0.57 NA NA NA -2647 0.46 0.30 NA NA

I am attempting to generate a series of plots with cubic B-Splines and the code I am using is generating the following error message that I have been unable to resolve. Error in grid.Call.graphics("L_xspline", x$x, x$y, x$shape, x$open, x$arrow, : add_point - reached MAXNUMPTS (25200) This is the code that generates this error. plot.results <- function(dfn) { dt <-

I am having difficulty loading a points shapefile (1600 records) I generated from ArcGIS as a PointsSpatialDataFrame using the maptools readShapePoints function. library(maptools) surveypts <- readShapePoints(system.file("C:/temp/survey.shp", package="maptools")[1] Error in getinfo.shape(filen) : Error opening SHP file I've tried numerous re-exports and renames of the

Dear all, I've been puzzled why I not able to load a shapefile from a connection. Does anyone here can give a reasonable answer? When I try the following script I got this error: Error in getinfo.shape(fn) : Error opening SHP file #Reproduction temp <- tempfile() download.file("http://metodologiapolitica.com/download/Brasil.zip",temp) shp.br <- readShapePoly(unzip(temp,

Now available on CRAN is a revised version of the shapefiles package for reading and writing shapefiles in R. New additions, courtesy of others, include the ability to convert a simple R data frame of points, polylines or polygons to a shp format list, which can then be written out to a shapefile with write.shp. There is also a function to convert the read.shp shp format list to a simple data