J'ai un petit problème de la représentation de mon axe des abscisses. Je représente des classes sédimentaires donc j'aimerais qu'elles soient classées de manière croissante.
Donc voici mon code pour mettre les bon "levels dans mon datafit pour les sédiments :
Code : Tout sélectionner
datafit$CODE_SEDIM <- factor(datafit$CODE_SEDIM, levels = c("M", "FS", "CS", "G", "P"))
J'ai essayé aussi avec ça :
Code : Tout sélectionner
levels(datafit$CODE_SEDIM) <-c("M", "FS", "CS", "G", "P")
Ensuite je fais un gam:
Code : Tout sélectionner
full1_gam<- gam(TDI ~ s(FA, k=3) + s(Salinity, k=3) + Ti + s(wave_stress, k=3)
+ s(friction_velo, k=3) + s(bed_stress, k=3)
+ CODE_SEDIM+ Year, data=datafit, family = gaussian)
Puis je veux des graphes avec mes résidus partiels pour chaque variable donc je fais tourner le code suivant:
Code : Tout sélectionner
library(mgcViz)
df_fig <- getViz(full1_gam)
fig_sediment <- plot(df_fig,
rug = F,
residuals = F,
select = 7) +
l_ciBar(lwd = 1, lty = 1, width = 0.2)+
l_fitPoints()+
labs(x = "Sediment type",y = "Partial residuals")
Et c'est là où j'ai un bug, mes "sédiments" s'affichent par ordre alphabétique donc CS, FS, G, M, P et non M, FS, CS, P, G
Pourtant quand je fait juste un plot "normal" les sédiments sont bien "rangés" dans le sens que je veux...
Code : Tout sélectionner
plot(TDI~CODE_SEDIM, data=datafit)
Je comprends vraiment pas d'où viens le souci.
Quelqu'un aurais une idée pour que j'obtiennes ce que je souhaites?
Cordialement
Cyrielle