皆さん、私は 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)
次に、 function の結果として、次のようplotに表示されます。
さて、私の問題を説明させてください...これは単なるテストであり、2ポイントしかありませんが、私が直面しているタスクには約1000〜1200ポイントがあり、約10回も実行する必要があります。私のコンピューターを使用してこのシミュレーションを行うことはできません。時間がかかるため、代わりに次のサービスのいずれかをレンタルする予定です: Amazon EC2 (Amazon Elastic Compute Cloud) .
そこでは、Linux/UNIX シェル (SO: Ubuntu だと思います) を介してそれを行う必要があります。私の質問は、コードの最後の 2 行に関連しています。Kenvインポートするためにファイルまたは何かにエクスポートするにはどうすればよいですか?それplotを私のコンピューターに入れます。
ご理解いただけたかと思いますが、ご不明な点がございましたらお知らせください。