Bonjour, n'étant pas une pro du codage sous R, je voulais savoir s'il y avait des améliorations à apporter à mon code pour gagner en temps de calcul lorsque je fais tourner mon programme ! L'objectif de ce code est de faire une validation croisée de mon modèle. J'ai 3 années de données, j'en sélecti...
J'ai résolu mon problème de message d'erreur en changeant l'écriture de nlminb : regroupement des paramètres en vecteur et appel de mon jeux de données observées avec un nom différent de celui des fonctions et ça tourne!
Bonjour, Je cherche à optimiser 3 paramètres d'un modèle : r0, Pc0, n_cum. Seulement la fonction nlminb me renvoie l'erreur suivante : Error in cbind(myData, Pcum = c(rep(0, ceiling(n_cum)), rep(NA, length(myData[, : argument "n_cum" is missing, with no defa...
Génial, c'est exactement ce qu'il me fallait, un grand merci! (il y avait juste un "n" qui trainait en trop dans le code de la fonction) J'ai un petit souci maintenant dans le code pour optimiser mon paramètre r0 : i_reel=c(0,0.1,0.15) s2_B <- function(myData,r0,i_reel) ...
Bonjour, j'aurais besoin d'aide dans la construction de ma fonction. Je cherche à simuler l'apparition de tache (n) en fonction de 2 variables (Temp et Pluie) et d'un paramètre (r0, à optimiser, valant 0.5 environ) mais elle dépend également de la valeur de n à la date précédente (appelée n_t_1). J'...
Bonjour, j'ai cherché à savoir si certains de mes capteurs enregistraient une humidité plus importante que d'autres, j'ai donc fait un test de Friedman qui m'a donné un résultat significatif. Je voulais ensuite effectuer un friedmanmc pour pouvoir classer mes capteurs. > donneeHR date PV INRA ID2LGM...
Bonjour, merci pour l'info c'est exactement ce qui me convient! Par contre, j'ai un peu de mal à coder la partie aléatoire de mon modèle. Tout d'abord, j'ai un jeu de données comme ceci : > donneeLAI Bloc Parcelle Dens Con date LAI Variete 1 3 15 D2 I 8/6 2.01 ES 2 3 14 D2 I 8/6 1.77 NK 3 3 13 D2 I ...
Bonjour, je me permets de faire remonter un peu mon sujet, j'ai fait des recherches et découvert l'existence des anovas à mesures répétées, seulement je n'ai trouvé aucune doc pertinente sur leurs programmations sous R, qqun aurait-il le nom de la fonction ou de la doc qui lui est associée? Merci d'...
Bonjour, J'ai utilisé un dispositif en split-splot à 4 facteurs pour relever une variable quantitative LAI : Bloc : 3 blocs Con : une parcelle irriguée, l'autre sèche Dens : une densité D1 et l'autre D2 Variete : 5 variétés Je me suis aidée d'un topic du forum sur lequel j'ai trouvé cette doc : http...
Je suis allée fouillée dans les famille "apply" et j'ai trouvé mon bonheur avec un sapply, plus d'erreur et ca m'évite de réécrire ma fonction (qui est bien plus compliquée que l'exemple donné!)
Dsl si mon exemple n'était pas suffisamment clair, je le refais: xvirt_calc_B=function(d_B, N_B){ if (d_B > (10 + N_B)) { return(0) } else if (d_B >= 8){ return ((sqrt(100) - d_B + N_B)^2) } else if (d_B >6&...