7

アフリカの背景地図を使用して、研究サイトのポイントをプロットしようとしています。2 つを個別に作成することはできますが、それらを重ね合わせるのに苦労しています。

私が使用しているアフリカの地図は、maplibrary.org の Esri シェープファイルです。https://www.dropbox.com/s/etqdw3nky52czv4/Africa%20map.zipのドロップボックスから入手できます。ポイントはテキスト ファイルに保存されており、ドロップ ボックスからも入手できます。 https://www.dropbox.com/s/scvymytjsr5pvaf/SPM-437-22Nov12.txt . 彼らは、マラリア原虫の分子薬剤耐性に関する研究に言及しています。色が薬剤耐性遺伝子マーカーを持つ寄生虫の割合で、サイズがテストされた寄生虫の数になるようにプロットしたいと思います。

ポイントを個別にプロットする:

qplot(Longitude, Latitude, data = d.spm.437, colour = Frc437, size = Tot437)

アフリカの地図をプロットする:

library(maptools)
africa = readShapePoly("Africa.shp")
africa.map = fortify(africa, region="COUNTRY")
qplot(long, lat, data = africa.map, geom="path", group=group)

ポイントの表示を維持しながら、これら 2 つを一緒に配置するための助けをいただければ幸いです。

4

1 に答える 1