Modérateur : Groupe des modérateurs
Code : Tout sélectionner
# ne connaissant pas vraiment tes données je te propose ça
# pour le mois de juillet, data etant le tablmeau de données
plot(data$dateJ[format(data$dateJ,"%m")=="07"],data$donne[format(data$dateJ,"%m")=="07"])
# pour le premier Juillet
plot(data$dateH[format(data$dateJ,"%d/%m")=="01/07"], data$donne[format(data$dateJ,"%d/%m")=="01/07"])
# pour le complet
plot(data$dateH,data$donne)
Code : Tout sélectionner
ozone <- read.table("ozone.txt", header = TRUE, sep = " ")
ozone$Jour <- as.Date(as.character(ozone$JOUR) , format = "%d/%m/%y")
ozone$Heure <- format(ozone$HEURE, format = "%H:%M")
Time <-as.POSIXct(paste(ozone$Jour, ozone$Heure))
> ozone
JOUR HEURE Ozone Jour Heure
1 01/07/03 00:15 79 2003-07-01 00:15
2 01/07/03 00:30 77 2003-07-01 00:30
3 01/07/03 00:45 77 2003-07-01 00:45
4 01/07/03 01:00 77 2003-07-01 01:00
5 01/07/03 01:15 77 2003-07-01 01:15
6 01/07/03 01:30 76 2003-07-01 01:30
7 01/07/03 01:45 75 2003-07-01 01:45
8 01/07/03 02:00 77 2003-07-01 02:00
9 01/07/03 02:15 75 2003-07-01 02:15
10 01/07/03 02:30 71 2003-07-01 02:30
11 01/07/03 02:45 73 2003-07-01 02:45
12 01/07/03 03:00 69 2003-07-01 03:00
13 01/07/03 03:15 3 2003-07-01 03:15
14 01/07/03 03:30 40 2003-07-01 03:30
15 01/07/03 03:45 13 2003-07-01 03:45
16 01/07/03 04:00 17 2003-07-01 04:00
17 01/07/03 04:15 73 2003-07-01 04:15
18 01/07/03 04:30 71 2003-07-01 04:30
19 01/07/03 04:45 69 2003-07-01 04:45
20 01/07/03 05:00 63 2003-07-01 05:00
21 01/07/03 05:15 66 2003-07-01 05:15
22 01/07/03 05:30 67 2003-07-01 05:30
library(lattice)
library(gam)
xyplot(Ozone ~ Time, data = ozone, aspect = "xy",
panel = function(x, y){
# trace une grille
panel.grid()
# trace une ligne joignant les points
panel.xyplot(x, y, type = "l", lty = "dotdash")
# trace les points
panel.xyplot(x, y, cex = 1.2, pch = 21, fill = "light blue", col = "blue")
# lissage spline
m <- mean(x)
X <- x - m
fm1 <- gam(y ~ s(X))
llines(x, fitted(fm1), col = "red")
# lissage loess
fm2 <- gam(y ~ lo(X))
llines(x, fitted(fm2), col = "green")
},
key = list(space = "right",
text = list(c("data", "spline", "loess")),
lines = list(col = c("blue", "red", "green"),
lty = c("dotdash", "solid", "solid")))
)
ensuite il faudrait que j applique la regression grace a la fonction loess et peut etre spline. pa dans R mais j avoue que j ai un peu de mal car je n ai pas l 'habitude , je ne travail pas sur ce langage .
Geissmann Gael a écrit :Merci pour ces réponse Juste quelque petites questions :
-pouvait vous me donner une url du package pour R et windows
-Y a t il uine autre solution pour réaliser cette regression non para. sans package ?
Merci beacoup
Retourner vers « Questions en cours »
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité