Passage entre script et fonction

Version imprimable

Bonjour,

J'ai un soucis car je n'arrive pas à faire passer mes deux objets data et modele dans une fonction que j'ai écrite.

Alors qu'en copiant les lignes de la fonction dans la commande, tout marche.

Voici le main :

Code:

1
2
3
4
5
6
7
source("mon chemin de la fonction backtest")	
 
data_modelisation=read.table("mon chemin",header=TRUE, sep="\t", na.strings="NA", dec=",", strip.white=TRUE)	
modele=glm(RA ~ varA varB ,family=binomial(logit), data=data_modelisation)
 
data=read.table("chemin",header=TRUE, sep="\t", na.strings="NA", dec=",", strip.white=TRUE)
backtest(data,modele)

Et la fonction :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
Backtest = function (modele,data)
{
 
		#Prévision des probabilités
		prevision=predict(modele,newdata=data,type="response")
		data$prevision <- with(data, prevision)
		summary(data)
 
		#Prévision des modalités
		pred.moda <- factor(ifelse(prevision > 0.5,1,0))
		data$pred.moda <- with(data, pred.moda)
 
		RA_nb_model=sum(data$prevision)/(12*length(data$prevision)) 
		RA_nb_reel=sum(data$RA)/length(data$prevision)
 
		resultat =c(RA_nb_model,RA_nb_reel)
		return(resultat)
 
 
 
}

Merci à ceux qui pourront m'aider

Bonjour,

Citation:

Envoyé par Stepsbysteps

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
Backtest = function (modele,data)
{
 
		#Prévision des probabilités
		prevision=predict(modele,newdata=data,type="response")
		<s>data$prevision <- with(data, prevision)</s>		summary(data)
 
		#Prévision des modalités
		pred.moda <- factor(ifelse(prevision > 0.5,1,0))
		<s>data$pred.moda <- with(data, pred.moda)</s>
		RA_nb_model=sum(<s>data$</s>prevision)/(12*length(<s>data$</s>prevision)) 
		RA_nb_reel=sum(data$RA)/length(<s>data$</s>prevision)
 
		resultat =c(RA_nb_model,RA_nb_reel)
		return(resultat)
 
 
}

Merci à ceux qui pourrront m'aider

De rien

12/03/2013, 14h38
Stepsbysteps
Merci mais j'ai encore cette erreur,

Code:

pas de méthode pour 'predict' applicable pour un objet de classe "data.frame"

Et pourquoi ne peut-on pas mettre ces lignes dans le code ?

J'avais aussi cette manipulation que je n'avais pas affichée :
Code:

1 2 #Passe les variables en facteur for (i in 4:16) data[,i] = factor(data[,i])
12/03/2013, 17h17
Sengar
Citation:

Envoyé par Stepsbysteps

Merci mais j'ai encore cette erreur,

Code:

pas de méthode pour 'predict' applicable pour un objet de classe "data.frame"

Cette erreur signifie comme elle le décrit bien qu'au lieu d'appliquer la fonction predict à un modèle, tu l'a appliqué à... un data.frame!
La seule explication est que tu as du te tromper dans l'ordre de tes paramètres quand tu as appellé la fonction! Je regarde et....

Effectivement je voir l'erreur dans ta main que je n'avais pas regardé :mouarf:

La declaration de la fonction est:
Code:

1 2 Backtest = function (modele,data)
Or dans ta main tu appelles backtest par :
Code:

1 2 backtest(data,modele)
Donc tu t'es trompé dans le sens des paramètres, et R comprend :
Code:

1 2 Backtest(modele=data,data=modele)
D'où l'erreur.

Citation:

Et pourquoi ne peut-on pas mettre ces lignes dans le code ?

Parce qu'elles sont inutiles et que l'emploi de with est totalement injustifié.
Cela dit après coup peut être que ce n'était pas ca qui faisait bugger mais juste le changement de sens des paramètres donc autant pour moi!

Mais il faut bien comprendre que with sert à travailler sur les variables d'un environnement (ici ton data frame data).
Donc quand tu fais with(data, unautretrucpasdansdata) au mieux ca equivaut à unautretrucpasdansdata, au pire ca peut te faire bugger ton programme.

edit: Au passage attention, car nommer tes paramètres appellés comme les paramètres de ta fonction (data=data) est susceptible de produire des bugs si tu retravailles ce paramètre après dans ta fonction ( ca peut parfois poser des problèmes d'environement quand tu imbriques des fonctions entre elles, et c'est d'ailleurs pour ca que je me suis dis que le with(data,) arrangeait pas trop les choses ^^)
12/03/2013, 17h45
Stepsbysteps

Ah! Super merci, j'ai été un peu découragé mais pourtant j'étais pas loin de mon coup.

Le with marche, merci de ces bonnes explications