Bonjour,

Je souhaite créer une nouvelle variable nommée eruptionsnew à partir d'une variable existante nommée eruptions. eruptionsnew ne prendra en compte que les eruptions de plus de 3 minutes. Cette étude prend part dans le cadre d'une analyse sur les geysers, leurs durées et leurs temps d'attente.

Voila le code de l'analyse :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
require(stats); require(graphics)
f.tit <-  "faithful data: Eruptions of Old Faithful"
 
ne60 <- round(e60 <- 60 * faithful$eruptions)
all.equal(e60, ne60)             # relative diff. ~ 1/10000
table(zapsmall(abs(e60 - ne60))) # 0, 0.02 or 0.04
faithful$better.eruptions <- ne60 / 60
te <- table(ne60)
te[te >= 4]                      # (too) many multiples of 5 !
plot(names(te), te, type="h", main = f.tit, xlab = "Eruption time (sec)")
 
 
X11 (); plot(faithful[, -3], main = f.tit,
     xlab = "Eruption time (min)",
     ylab = "Waiting time to next eruption (min)")
lines(lowess(faithful$eruptions, faithful$waiting, f = 2/3, iter = 3),
      col = "red")
J'ai testé différentes choses mais rien n'a réellement fonctionné.
Merci par avance de votre aide.