Hi there,
Just wondering why my post was rejected?
cheersRachel
Subject: repeating autocovariate functions
From: r-help-owner@r-project.org
To: moyble@hotmail.com
Date: Sat, 2 Feb 2013 02:56:27 +0100
Message rejected by filter rule match
--Forwarded Message Attachment--
Date: Fri, 1 Feb 2013 17:56:14 -0800
From: moyble@hotmail.com
To: r-help@r-project.org
Subject: repeating autocovariate functions
Hi there,
I would like to repeat an autocovariate term calculation using 30 different
neighborhood sizes. Then I would like to run a nominal logistic regression
on the generated autocovariate values and their respective neighborhood
sizes to see which would be most appropriate to use in the final calculation
of my autocovariate term.
I have a matrix of x,y values:
x y
174.7173 -35.967
174.7166 -35.9649
174.7174 -35.968
174.7418 -35.9678
174.741 -35.9672
174.7395 -35.9671
(and 150 more)
To calculate the autocovariate terms my code is as follows:
library(spdep)
Taranga <- read.csv(file.choose()) #contains xy coordinates
xy <- cbind(Taranga$x, Taranga$y)
acinvb <- autocov_dist(Taranga$burrows.pa, xy, nbs=1,zero.policy = TRUE,
type="inverse")
datfrm <- data.frame(autocov=acinvb, nbs=1)
write.table(dfrm, file = 'results.csv',sep=",",row.names = FALSE)
acinva <- autocov_dist(Taranga$burrows.pa, xy, nbs=100,zero.policy = TRUE,
type="inverse")
dfrm <-data.frame(autocov=acinva, nbs=100)
names(datfrm) <- NULL
write.table(datfrm, file = 'results.csv',sep=",",append=TRUE, row.names =
FALSE, col.names=FALSE)
I want to repeat this function, each time adding 100 to "nbs"
Then I will run a model from the resulting "results.csv"
thanks in advance for any insight!
Rachel
--
View this message in context: http://r.789695.n4.nabble.com/repeating-autocovariate-functions-tp4657353.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]