280
Le logiciel R - Maˆ ıtriser le langage
5
nbindiv <- nrow (X)
6 browser ()
7
if ( intercept == TRUE ) X <- cbind ( rep (1 , nbindiv ),X)
8
betachap <- solve (t(X) % * % X) % * % t(X) % * % Y
9
curve ( betachap [1]+ betachap [2] * x , add = TRUE )
10
11 return ( betachap )
12 }
Sourcez le fichier contenant le bout de code pr´ ec´ edent (par exemple au
moyen de l’instruction source(file.choose())), puis tapez :
mc(X=cars[,2],Y=cars[,1])
Vous constatez alors que le programme stoppe et que l’on peut examiner
le contenu des variables locales situ´ ees avant browser(). Par exemple, tapez
nbindiv.
Remarque
En entrant la lettre n (pour next) de fa¸ con successive, on peut inspecter
le code et le contenu des variables locales de fa¸ con s´ equentielle. Pour quitter
le mode d’inspection, tapez Q.
Voici un aper¸ cu d’une possible session de d´ ebogage :
mc(X=cars[,2],Y=cars[,1])
Called from: mc(X = cars[, 2], Y = cars[, 1])
Browse[1]>nbindiv
[1] 50
Browse[1]> betachap
Erreur: Objet "betachap" non trouvé
Browse[1]> n
debug: if (intercept == T) X <- cbind(rep(1, nbindiv), X)
Browse[1]> n
debug: betachap <- solve(t(X) % * % X) % * % t(X) % * % Y
Browse[1]> n
debug: curve(betachap[1] + betachap[2] * x, add = T)
Browse[1]> betachap
[,1]
[1,] 8.2839056
[2,] 0.1655676
Browse[1]> Q
>