Groupe des utilisateurs du logiciel R

Dans le cas d'une fonction de variance mu^z (où mu est l'espérance), lorsque z est différent de 0, 1 ou 2 la quasi-vraisemblance s'écrit (mu^-z)*(mu*y/(1-z)-(mu^2)/(2-z)) (McCullagh & Nelder p.326). Y a-t-il un package R qui permette d'estimer les paramètres d'un modèle avec cette fonction? Je ...

Je cherche une fonction R pour maximiser une vraisemblance quand la loi est connue Plusieurs fonctions d'optimisation sont disponibles, par exemple : optim, nlm ou nlminb. Il suffit de specifier la (log-)vraisemblance comme fonction a optimiser. Des exemples de code sont diponibles dans la rubrique...

Je trouve que le graphe suivant est assez clair et permet peut-être de mieux révéler les stuctures (absentes dans cet exemple !) qu'un graphe 3D: oui barchart peut être une bonne alternative mais je cherchais une fonction R équivalente de la fonction bar3 de Matlab, qui ne semble donc pas exister ....

Est-ce que tu peux aussi poster le genre de graphique que tu veux obtenir, quitte a le generer avec matlab ou a le faire a la main ... J'ai placé un exemple sur le lien : ftp://ftp.cirad.fr/pub/group-r/groupe-r/Documents/bar3_ex.bmp Le code Matlab est tout simple : bar3(Z) xlabel('Time') ylabel('Ag...

Renaud Lancelot a écrit :Pour ce genre de graphe je ne ferais pas de 3D mais un graphe trellis. Voir ?barchart dans le package lattice.

Renaud

en fait, j'ai essayé précedemment, ca devient illisible si bcp de structures, alors que la visu est excellente avec bar3.

M.

http://addictedtor.free.fr/graphiques/search.php?q=3d oui j'avais vu la fonction scatterplot3d ou cloud(..., type = h) en parcourant les galeries, mais pour ce que je cherche à représenter (évolution de structures en classes d'âge en fonction du temps) le rendu n'est vraiment pas terrible. Pour une...

Bonjour Je cherche a construire un graphique 3D du meme type que ce que font persp (graphics) et wireframe(lattice), par exemple : # ?wireframe library(lattice) wireframe(volcano, shade = FALSE, aspect = c(61/87, 0.4), light.source = c(10,0,10)) mais avec une represen...

Donc personne ne sait s'il existe un moyen de se connecter à une telle base...une option qui permettrait d'indiquer le fichier workgroup ? En regardant rapidement ds les docs sur CRAN et RODBC, je n'ai rien trouve sur les mdw (mais peut-etre y a t il des infos qq part). Je te conseille de contacter...

Effectivement, comme vous avez dit plus haut dans ce que j'ai mis en gras, ce dont j'ai besoin, c'est de faire des stats sur des données issues de simulation . La simulation est faite par un autre logiciel, lequel produit des fichiers texte ou xml, que je dois exploiter comme entrées. Ces fichiers ...

Bonjour à tous, Auriez-vous connaissance d'une étude comparative entre ces différents outils, ou pourriez-vous me donner vos propres points de vue ? Bonjour, je ne parle que de R et Scilab (je ne connais pas maxima). Scilab est davantage un logiciel de modelisation (systemes dynamiques, optimisatio...

base$obesite[base$obesite<30]<-0 base$obesite[base$obesite>=30]<-1 Bonjour, une approche + simple est d'utiliser ifelse : base$obesite <- ifelse(base$imc < 30, 0, 1) On peut emboiter des ifelse, par exemple : base$obesite <- ifelse( base$imc < 30, 0, ifelse(base$imc >= 30 & base...

Bonjour, Ma question est en fait de savoir à partir de quel effectif doit on utiliser le test de fisher au lieu du test du chi2. Murielle Le pb n'est malheureusement pas simple et comme le dit Pierre Bady les usages varient. En complément de la règle citée par Samir, il y a celle citée dans le livr...

Practical regression and anova using R:

http://cran.r-project.org/doc/contrib/Faraway-PRA.pdf

Super, merci Renaud

Bonjour Je cherche à représenter un lissage séparé de 2 nuages de points avec xyplot (lattice). J'ai le code suivant : xyplot( x = p ~ unitdec, groups = sex, panel = function(x, y, ...){ panel.superpose(x, y, ...) panel.loess(x, y, lty = 2, span = 2 / 4, col = "red&...