

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

