Le portail francophone de la géomatique

Bonjour,
je vous explique mon travail en 2 mots,

Je souhaiterais réaliser une carte de probabilité avec partir de 2 shapefiles.
Donc j'ai réalisé dans GRASS une classification du sol (6 classes en tout, donc une classe eau et une classe rizière). Cette classification, je l'ai importé dans R sous forme de raster. Ensuite j'ai importé deux shapefiles qui me calcule pour les fermes de ma région, la distance à l'eau et la distance à la rizière la plus proche.

library(raster)
library(rgdal)
library(maptools)
myInShapeRice = readShapeSpatial ("C.....shp")
myInShapeDistWater=readShapeSpatial("C:...shp")
myDuRaster = raster("C:/temp/essai.gtif")

myProbR = log10(myDuRaster) + 2


J'ai ensuite réalisée une regression pour réaliser ma carte


myTempR = 0.523 + 5.329*(myInShapeDistWater) + 3.154*(myInShapeRice)
myProbR = exp(myTempR)/(1+exp(myTempR))

Le premier probleme que j'ai rencontré c'est qu'il ne veut pas lire mes shape directement, j'ai donc du passer par la création de data.frame
vecteur1=date.frame(myInShapeDistWater)

Le vrai probleme que je rencontre, c'est maintenant j'aimerais bien passer à un raster
writeRaster(myProbR, "C:...gtif")
Mais on me dit que je ne peux pas utiliser de data.frame
Pourriez-vous m'aider?