Tester un modèle

Bonjour,

Je souhaiterais avoir un peu d'aide pour mon programme star.

A partir de ma base de donnée, j'ai effectué un sous échantillonnage aléatoire créant ainsi deux jeux de données : Echantillon1 et Echantillon2

A partir de échantillon 1, j'ai pu réaliser un modèle linéaire multiple de ce type :
model.regmult=lm(FishCd ~ FishTaille+FishMG+WaterCd,data=Echantillon1)

je voudrais confronter mon modèle aux données du sous-échantillon 2

et faire un graphique où sont représentées les données réelles de FishCd de l'échantillon 2 et celles prédites par mon modèle.

et je voudrais savoir si mon modèle colle bien.

Quelqu'un pourrait me dire comment écrire ce type de programme ?

merci à vous.

Bonjour,

Il faut utiliser la fonction predict() pour générer les prédiction du modèle. Ensuite tu peux afficher ça sur un plot comme tu le veux.
Il n'y a pas vraiment de programmation à faire, tout existe déjà.

j'ai essayé cette formulation avec predict

newdata<-subset(Echantillon2B, select= c(FishTaille, FishMG, WaterCd))
print(newdata)
predict(model.regmult, newdata, se.fit=FALSE, scale = NULL, df = Inf, interval = "confidence", level = 0.95)

comment faire le plot des données réelles "FishCd" d'Echantillon2B en x versus les valeurs "FishCd" prédites selon le modèle à partir de l'échantillon 2B ?

Merci pour votre aide