similar to: se.contrast

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

Hi all, Regular and avid readers of this column will know that Don Driscoll and I have recently posted two messages requesting assistance concerning an apparent failure of "se.contrast" to produce an se for a contrast. So far, an ominous silence rings in our ears, but read on Gentle Reader, and see if even the machinations of "debug" doesn't stimulate you to respond with a

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

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

Hello! I'm having some problems with  pit function in ensembleBMA. When I run on mu data I have the same warnings as below, when I type the example(pit). Any suggestion on how to resolve it? > library(ensembleBMA) Loading required package: chron > example(pit) pit>   data(ensBMAtest) pit>   ensMemNames <-

Hello Chris,  I 've found two other issues  with MAE and CRPS, giving warning  when running examples. I've the same issue on my data.  Hope that you could find some time to take a look here. Thank you Anna > library(ensembleBMA) Loading required package: chron > example(MAE) MAE>   data(ensBMAtest) MAE>   ensMemNames <-

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:

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

hello, I have a graph and i use qgraph package to calculate centrality parameters. Now I want to know the maximum value of shortest path for each vertex with discarding the Inf value in short pathes. For this I use the ShortestPathLengths of centrality function in qgraph. but when I want to get the maximum the result is wrong. here is my code: cen<-centrality(Q)