補間された地球化学マップを作成するプロセスを自動化しようとしています。基本的に列 #13 から始まり、#67 までループするループを作成しました。
ただし、IDW コードの一部では、それぞれの地球化学パラメーターをマップするために、現在の列のヘッダーが必要です。
例えば; 列 #13 のヘッダーは「Ag_ppm」であるため、元のコードは次のようになります。
LogSr.idw = idw(log10("Ag_ppm") ~ 1 , locations=NGSA.SPDF, newdata=NGSA.grid,
maxdist=15000, nmin=4)
次を使用して、各ループ内で対象の列のヘッダーを取得できました (ここで、i はループごとに 1 ずつ増加します。13、14、15 など)。
coln <- colnames(NGSA.df[i])
ただし、単に "Ag_ppm" を coln に置き換えると、行は特定のエラーで失敗します。貼り付けを含むさまざまなアプローチを試しましたが、すべて同じエラーが発生します。
LogSr.idw = idw(log10(coln) ~ 1, locations=NGSA.SPDF, newdata=NGSA.grid,
maxdist=15000, nmin=4)
log10(coln) のエラー: 数学関数への数値以外の引数
現在の方法を使用した合理的に単純なアプローチはありますか?