116 résultats trouvés

Aller à la recherche avancée

par Matthieu Lesnoff
25 Mar 2011, 15:27
Forum : Questions en cours
Sujet : Moyenne ajustée
Réponses : 10
Vues : 15563

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 : 6351

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 : 3979

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 : 2408

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 : 3130

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 : 10
Vues : 9588

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 : 10
Vues : 9588

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 : 3664

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 : 3664

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 : 2580

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 : 3967

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 : 4490

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 : 2901

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...
par Matthieu Lesnoff
09 Fév 2010, 14:14
Forum : Questions en cours
Sujet : format date
Réponses : 4
Vues : 7845

Tu n'as pas reproduit exactement le code que j'avais mis. Il faut d'abord passer ton z en caractères (as.character(z)), puis utiliser as.Date (et non as.date) puis format (attention de bien respecter l'argument format)
M.
par Matthieu Lesnoff
04 Fév 2010, 11:46
Forum : Questions en cours
Sujet : format date
Réponses : 4
Vues : 7845

Re: format date

comment transforme t on une date de ce format 20090729 en un format du style 7Juil2009 ? je pense qu'il s'agit plutôt de 29Juil2009. Quelque chose qui y ressemble (sans utiliser le package date) : > z <- as.Date("20090729", format = "%Y%m%d") > z [1] "2009-07-29&quo...

Aller à la recherche avancée