bonjour, suite un petit test et à la relecture du code initialement proposé, j'ai ajouté deux modificatiuons: 1- j'avais oublié de multiplier la matrice kinship par 2. 2- la fonction chol renvoie le triangle supérieur ... c'est le triangle inférieur qui nous intéresse donc il faut transposer le résu...
bonjour, Random effects: ~1 | id Variance list: list(matkin) id resid Standard Dev: 26.12692 1.023894e-03 % Variance: 1.00000 1.535795e-09 ici ... % Variance: 1.00000 1.535795e-09 Les résultats différents quand même pas mal des autres méthodes... J'ai pas encore compris tous le code (ni les ...
rebonjour, après une petite recherche dans le code, on trouve ceci: pour lmekin (coxme), on a ceci: yhat <- as.vector(zstar1 %*% rcoef + X %*% fcoef) ... residuals <- Y - yhat et pour lmekin (kinship), on a ceci: fitted <- c(X %*% lfit$coef) residuals <- Y - fitted donc, les résidus ...
bonjour, la function coxme me semblerait plus 'fiable' (package de l'auteur de lmekin). Mais bon, ... les variances du modèle semblent quand même un peu bizarre! Random effects Group Variable Std Dev Variance id Vmat.1 4.401188e+27 1.937045e+55 Residual error= 3.979124e-07 en plus, tous les résidus ...
rebonjour, pas tout à fait exact: input arguments and form the random effects, comparison of its output with lme helped validate those operations. It is possible to specify some models in this framwork that can not be fit with lme, Effectivement, il y a certains modèles qui ne sont pas réalisable av...
bonjour, la premier version (dans kinship) n'est plus distribuée (installation possible à partir d'archive en local), la seconde (dans coxme) est plus flexible ... mais, je n'ai pas réussi à la faire tourner correctemment :/ sinon, on peut également utiliser la function lme avec un effet random de t...
bonjour, le première argument est un objet de type formula. On peut générer un nouvelle formule avec les fonctions paste et as.formula . R> dat <- data.frame(A=1:5, B=factor(LETTERS[5:1])) R> mavariable <- "B" R> gx.sort.df(as.formula(paste("~",mavaria...
bonjour, Ce qui m'interpelle, c'est qu'un package tel que vegan avec des applications spécifiques pour les analyses multivariées n'aborde pas le thème des proportions alors que c'est quand même une forme commune de données en écologie. ce type de données est peut-être plus observé en génétique des p...
Bonjour, je ne dis pas que c'est pas bon (le cylindre, si il a le bon diamètre, peut rentrer dans le carré), mais il faut juste bien savoir pourquoi on le fait (?) En plus, l'utilisation de cette distance n'est pas forcément anodine ... c'est une semimétrique, non-euclidienne (on doit pouvoir la ren...
bonjour, pourquoi vouloir utiliser la distance de Bray-Curtis à tout prix? la méthode d'analyse ne doit pas forcer les données. en principe, on fonctionne dans l'autre sens (?) ça me fait penser aux jeux d'encastrement (jeu d'éveil) ... lorsque le cylindre rentre pas dans la forme d'étoile ... Une m...
bonjour, Pour l'utilisation/la construction d'une dll simple , il me semble que l'échange R -> C implique l'utilisation d'objet de type vector (c'est à vérifier quand même, je suis un peu rouillé et mes compétences en C restent très approximatives :D). donc, pour manipuler un tableau dans ton code C...