118 résultats trouvés

Aller à la recherche avancée

par Matthieu Lesnoff
17 Juin 2011, 17:52
Forum : Questions en cours
Sujet : Calcul intervalle de confiance de paramètres estimés
Réponses : 5
Vues : 5567

Re: Calcul intervalle de confiance de paramètres estimés

A partir de régressions non linéaires, j'essaye de calculer les intervalles de confiance de mes paramètres estimés. Pour calculer les intervalles de confiance je fais: confint(model_pocmat2) La fonction confint calcule des IC de type "profile likelihood" ; si besoin et en alternat...
par Matthieu Lesnoff
04 Avr 2011, 18:33
Forum : Questions en cours
Sujet : objet formula et modification
Réponses : 8
Vues : 4052

les méthodes terms, all.vars, etc fournissent la totalités des termes, y en a t'il une spécifique pour récupérer ces noms de variables? je ne pense pas (?) car "1|id" est aussi une variable, mais vous pouvez y arriver facilement par : > f <- formula(y~x+z+e+(1|id)) > z <- ...
par Matthieu Lesnoff
25 Mar 2011, 15:27
Forum : Questions en cours
Sujet : Moyenne ajustée
Réponses : 10
Vues : 17770

Erreur dans marg_means(formula = ~flour, fm = fm) : impossible de trouver la fonction "ctabs" Désolé je n'avais pas vérifié ce vieux code qui ne marche plus. Essayez avec : - Fonction emm (à sourcer) : emm <- function(fm, formula, digits = 4, ...) { CALL <- match.call(...
par Matthieu Lesnoff
21 Mar 2011, 21:20
Forum : Questions en cours
Sujet : Porgrammation des moyennes ajusées Ex : ANOVA à 2 facteurs
Réponses : 4
Vues : 7220

Re: Porgrammation des moyennes ajusées Ex : ANOVA à 2 facteu

Sur SAS, il s'agit de la fonction LSMEAN Avec l'utilitaire "Rechercher" de ce forum et en tapant "lsmeans", vous tombez notamment sur http://forums.cirad.fr/logiciel-R/viewtopic.php?t=625&highlight=lsmeans . La fin du post propose une fonction de calcul des lsmeans ("ma...
par Matthieu Lesnoff
17 Mar 2011, 07:16
Forum : Questions en cours
Sujet : Graphique avec tracé d'isolignes
Réponses : 4
Vues : 4655

Re: Graphique avec tracé d'isolignes

Je souhaiterais représenter l'évolution de la salinité d'un estuaire en fonction de la profondeur pour plusieurs points. Voir aussi fonction contourplot du package lattice. Un exemple fictif ci-dessous : require(lattice) tmp <- expand.grid(x = 1:80, y = 1:80) tmp$z <- rnorm(nrow...
par Matthieu Lesnoff
04 Mar 2011, 19:03
Forum : Questions en cours
Sujet : tester la différence de deux moyennes issues d'un bootstrap
Réponses : 4
Vues : 2944

Maxime Hervé a écrit :je ne suis pas vraiment calé en théorie


Pour une très bonne intro (avec des exemples, dont le cas évoqué dans le message), voir :

Manly, B.F., 2007. Randomization, bootstrap and Monte Carlo methods in biology 3rd ed., Boca Raton: Chapman & Hall/CRC.
par Matthieu Lesnoff
04 Mar 2011, 18:58
Forum : Questions en cours
Sujet : RStudio et R
Réponses : 4
Vues : 4085

Bonjour, Vu le peu d'options modifiables du logiciel toujours pas trouvé le moyen de changer ça, mais je le conseille il est très bien fait ! oui moi aussi j'ai beaucoup aimé. Installation et lien avec R simplissimes (en tout cas sous XP) et interface performante. Le problème d'accents a été soulev...
par Matthieu Lesnoff
01 Fév 2011, 23:58
Forum : Questions en cours
Sujet : Courbe de survie
Réponses : 12
Vues : 12010

personnellemet je ne vois pas trop le point de discussion. La "fonction de survie" définie pour une variable durée de vie "T" est bien-sûr décroissante. Mais la fonction de mortalité cumulée (= fonction de répartion de T) est le complément à 1 de cette fonction de survie. Je pens...
par Matthieu Lesnoff
01 Fév 2011, 14:33
Forum : Questions en cours
Sujet : Courbe de survie
Réponses : 12
Vues : 12010

Oui mais je ne sais pas à quel endroit écrire cette instruction, j'ai un code du type: plot(survfit(Surv(donnees$temps.avant.succes, donnees$succes)~donnees$sexe)) Une autre piste possible : test1 <- list(time=c(4,3,1,1,2,2,3), status=c(1,1,1,0,1,1,0), x=...
par Matthieu Lesnoff
22 Oct 2010, 12:01
Forum : Questions en cours
Sujet : fonction predict avec option newdata sur betabin
Réponses : 8
Vues : 4039

Re: matrice de variance covariance des estimations des coefs

Quelqu'un saurait il comment on obtient cette matrice de variance-covariance ? C'est le slot @varparam > fm <- betabin(formula = cbind(y, n - y) ~ seed * root, random = ~ 1, data = orob2) > round(fm@varparam, digits = 4) (Intercept) seedO75 rootCUCUMBER seedO75:rootC...
par Matthieu Lesnoff
21 Oct 2010, 23:14
Forum : Questions en cours
Sujet : fonction predict avec option newdata sur betabin
Réponses : 8
Vues : 4039

Valdimir, difficile d'identifier le pb sans connaitre ton data.frame, le modèle que tu as ajusté et ton newdata Ton plan est il complet ou incomplet ? Dans le dernier cas, certaines combinaisons de modalités peuvent ne pas être estimables. Mais je pense que le pb est ailleurs. Avec un newdata, le co...
par Matthieu Lesnoff
21 Oct 2010, 12:13
Forum : Questions en cours
Sujet : gérer les années de 53 semaines
Réponses : 2
Vues : 3272

Sûrement pas optimal, mais voici la fonction que j'avais construite pour découper un cycle annuel (démarrant n'importe quel jour de l'année) en un nombre quelconque de "phases" (par exemple 52 "semaines") : fcal <- function(d1, nbphase) { if(is.character(d1&#...
par Matthieu Lesnoff
19 Oct 2010, 16:47
Forum : Questions en cours
Sujet : série temporelle, superposition de graph - décomposition
Réponses : 5
Vues : 4896

Re: série temporelle, superposition de graph - décomposition

Comment pourrais-je réaliser ce graph? pas sûr d'avoir bien compris la question mais peut être utiliser xyplot(lattice) et panel.superpose(lattice) après avoir créé une variable "groupe" (dans le data.frame qui contient les données) qui définit les périodes. Exemples simulés ci-dessous : ...
par Matthieu Lesnoff
02 Avr 2010, 11:57
Forum : Questions en cours
Sujet : Plot de moyennes avec erreurs standards
Réponses : 3
Vues : 5238

Re: Plot de moyennes avec erreurs standards

Je cherche a dessiner un graphique du type de celui présenté ici: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=72 mais où sont portées les moyennes et leurs erreurs standards au cours du temps. Lise en vrac, quelques exemples de graphiques avec plot, xyplot, dotplot et représentati...
par Matthieu Lesnoff
27 Mar 2010, 14:18
Forum : Questions en cours
Sujet : Remplissage intervalle de confiance de prévision (ARIMA)
Réponses : 1
Vues : 3296

Re: Remplissage intervalle de confiance de prévision (ARIMA)

Bonjour, Comment dois-je faire pour colorier l'intérieur de l'intervalle de confiance de prévision à 95% de mon modèle sans avoir utilisé les commandes armaSim ou armaFit ? après avoir utilisé la fonction plot, vous pouvez par exemple utiliser la fonction polygon (voir les exemples dans l'aide de l...

Aller à la recherche avancée