Hi there,
Just wondering why my post was rejected?
cheersRachel

Subject: repeating autocovariate functions
From: [email protected]
To: [email protected]
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: [email protected]
To: [email protected]
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]]

______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to