0

伙计们，我正在 Mac 中使用 R 处理一些数据，首先让我向您展示我的代码。

require(spatstat)
require(maptools)
setwd("/Whereever") # in that folder I have the files sis1993.txt and colombia.shp...

mydata <- read.table("sis1993.txt", header = TRUE)
attach(mydata)
summary(mydata) 
datos=read.table("sis1993.txt", header=T, dec=",", sep="\t")
summary(datos)
attach(datos)

S=readShapePoly("colombia.shp")
SP=as(S,"SpatialPolygons")
W=as(S,"owin")

sis1993=ppp(datos$x, datos$y, window=W)
unitname(sis1993)="meter"

sis1993=ppp(datos$x, datos$y, window=W, marks=m )

unitname(sis1993)="meter"
summary(sis1993)

Kenv <- envelope(sis1993,fun="Kest",nsim=199,nrank=5)
plot(Kenv,xlab="r",ylab="Khat(r)",cex.lab=1.6,cex.axis=1.5,main="K-Hat", cex.main=1.5)

然后，作为函数的结果plot，我看到了这个：

在此处输入图像描述

现在，让我解释一下我的问题......这只是一个测试，它只有2分，但我面临的任务有1000-1200分左右，我也必须做十次左右，多亏了我无法使用我的计算机进行此模拟，因为它会花费大量时间，而不是我计划租用以下服务之一：Amazon EC2（Amazon Elastic Compute Cloud）。

在那里，我必须通过 Linux/UNIX 外壳（所以：Ubuntu，我猜。），现在我的问题与我的代码中的最后两行有关，我如何导出Kenv到文件或其他东西以便导入它，然后plot它在我的电脑。

我希望你能理解我，如果你有任何问题，请告诉我。

4

1 回答 1

1

请参阅help(save)R。这会将内容保存到文件中。

然后将文件从亚马逊实例复制到您的PC（超出我的范围）。

然后看看help(load)如何重新加载它。

于 2013-07-22T16:08:39.953 回答