0

私は R を初めて使用する Python ユーザーです。現在、R パッケージ GWmodel を扱っています。

基本的な GWR の関数を見ると、次のようになります。

gwr.res <- gwr.basic(GenEl2004 ~ DiffAdd + LARent + SC1 + Unempl + LowEduc + Age18_24 + Age25_44 + Age45_64, data = Dub.voter, bw = 100,  kernel = "bisquare", adaptive = TRUE, F123.test = TRUE)

必要なのは、各変数の推定パラメーターの平均を収集し、それを bw (帯域幅) の任意の値のリストに追加することです。

Python 用語では、これは次のようになります。

LARentMean = []
SC1Mean = []
UnenmplMean = []
LowEducMean = []
Age18_24Mean = []
Age25_44Mean = []
Age45_64Mean = []


for i in range (20,400):
    gwrres = gwr.basic(GenEl2004 ~ DiffAdd + LARent + SC1 + Unempl + LowEduc + Age18_24 + Age25_44 + Age45_64, data = Dub.voter, bw = i,  kernel = "bisquare", adaptive = TRUE, F123.test = TRUE)
    a = gwrres(LARent).mean()    #a <- mean(gwrres$SDF$LARent)
    b = gwrres(SC1).mean()       #b <- mean(gwrres$SDF$SC1)
    c = gwrres(Unenmpl).mean()   #c <- mean(gwrres$SDF$Unempl)
    d = gwrres(lowEduc).mean()   #d <- mean(gwrres$SDF$LowEduc)
    e = gwrres(Age18_24).mean()  #e <- mean(gwrres$SDF$Age18_24)
    f = gwrres(Age25_44).mean()  #f <- mean(gwrres$SDF$Age25_44)
    g = gwrres(Age45_64).mean()  #g <- mean(gwrres$SDF$Age45_64)
    LARentMean.append(a)
    SC1Mean.append(b)
    UnenmplMean.append(c)
    LowEducMean.append(d)
    Age18_24Mean.append(e)
    Age25_44Mean.append(f)
    Age45_64Mean.append(g)
4

1 に答える 1

0

リストにループする lapply を使用できます。

l = lapply(20:400, function(i){
    gwr.basic(GenEl2004 ~ DiffAdd + LARent + SC1 + Unempl + LowEduc + 
        Age18_24 + Age25_44 + Age45_64, data=Dub.voter, bw=i,
        kernel="bisquare", adaptive=T, F123.test=T)
})

gwr.basic が出力として何を生成するかはわかりません。そのため、追加の行で平均値のみを取得する必要がある場合があります (より多くの情報を出力する場合)。

于 2015-08-22T11:17:22.960 に答える