Re: [R-sig-Geo] Selecting cells in a raster

2013-10-15 Thread Hugh Sturrock
Hi Manuel, raster[] will give you a vector of the raster values so #create a new raster of the same red and extent as cob_r forest_raster-cob_r # assuming forest pixels have a value of 1 forest_raster[forest_raster[]!=1]-0 plot(forest_raster) Does that work? Date: Tue, 15 Oct 2013 10:49:50

Re: [R-sig-Geo] output data from R to SaTScan

2013-08-26 Thread Hugh Sturrock
Hi James, Not sure how to do that, but check out SpatialEpi and Dcluster packages, I believe they can do a lot of what SaTScan does.. Cheers, Hugh From: roone...@tcd.ie To: r-sig-geo@r-project.org Date: Mon, 26 Aug 2013 16:42:55 +0100 Subject: [R-sig-Geo] output data from R to SaTScan Hi

Re: [R-sig-Geo] output data from R to SaTScan

2013-08-26 Thread Hugh Sturrock
$KNcluster) print(spplot(ED, KNcluster, main=Kulldorff's method, col.regions=c(gray(1), gray(.5), gray(.8 Thanks, James From: Hugh Sturrock [hughsturr...@hotmail.com] Sent: 26 August 2013 17:10 To: James Rooney; r-sig-geo@r

Re: [R-sig-Geo] nearest neighbour using distance in metres

2013-05-16 Thread Hugh Sturrock
Ross, The geodist function in the package GMT will compute distance in all sorts of units (m, km etc.) from lat longs. Hugh Date: Thu, 16 May 2013 14:16:03 +0100 From: rossah...@googlemail.com To: r-sig-geo@r-project.org Subject: [R-sig-Geo] nearest neighbour using distance in metres Hi

[R-sig-Geo] Spatially autocorrelated random effects in logistic regression.

2013-05-15 Thread Hugh Sturrock
Hi all, I am trying to model binomial data (infection status) of people living in households. I have run a logistic regression with a household level random effect to account for correlation between individuals within households. When I look at the value of the random effect (intercept) at each