headers
Emiliano Grossman | 4 Feb 19:26
Picon
Gravatar

exclusion de multiples observations

Bonjour, 


j'ai un jeu de données avec les scores de certains partis aux élections dans une variété de pays (chaque ligne représente un parti dans un pays dans une élection). 
Or la moitié de ces pays ne m'intéressent pas. Je voudrais donc les exclure, ce que je fais d'habitude avec

data2<-data[which(country!="countryname"),]

Ma question est la suivante: y aurait-il moyen d'exclure d'un coup dix ou vingt pays?
Ce que je fais actuellement, c'est que rajoute une ligne pour chaque pays que je veux exclure. Et ça me dérange pas de continuer à le faire. Je me demande juste s'il y a un moyen plus simple / élégant de le faire. 

bien à vous,


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 1 comment |
headers
Pierre de Larminat | 3 Feb 17:20
Picon
Gravatar

Superposer un histogramme, une courbe et leurs étiquettes

Bonjour à toutes et à tous,
Quelqu'un saurait-il faire le graphique ci-joint?
On y trouve le "nombre" en histogramme, avec les valeurs numériques en étiquettes,
et le "volume" en courbe (comme un "plot" de type "b"), avec aussi les valeurs numériques en étiquettes.

Dans un monde idéal, "volume" et "nombre" seraient représentés avec des couleurs différentes.

#Voici les données:
annee<-factor(2003:2010)
nombre<-c(44, 106, 143, 115, 130, 106, 85, 89)
volume<-c(6.8, 22.7, 14.6, 10.8, 26.9, 18.2, 19.4, 30.5)
data<-data.frame(annee, nombre, volume) 
str(data)
data

#Voici ma tentative sur la première partie du graphique:
nombre<-data$nombre
names(nombre) <- data$annee
bar<-barplot(nombre, ylim = c(0, max(nombre)*1.2))
text(bar, nombre+3, labels = nombre)

# Pour la suite, mes tentatives restent infructueuses.

Merci.
Pierre.

--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 5 comments |
headers
Jean-Michel Le Bot | 1 Feb 14:22
Picon

inverse de mls.eclate.multi ?

Bonjour
Existe-t-il une fonction qui ferait l'inverse de mls.eclate.multi ou 
sinon comment l'écrire ?
Ce que je voudrais faire est un tri à plat "à la Modalisa" (% de 
réponses et/ou de répondants) à partir d'une question à réponses 
multiples déjà dichotomisée par Limesurvey en autant de questions 
binaires que de modalités (ce que fait si je comprends bien 
mls.table.multi, mais à partir du format d'export de ces questions par 
Modalisa).
Merci
-- 
J e a n  -  M i c h e l   L e  B o t
Maître de conférences de sociologie
Université Rennes 2 - Département AES
Place du recteur Henri Le Moal
CS 24307
35043 Rennes Cedex
Chercheur au LAS/LARES (Université Rennes 2)
http://www.univ-rennes2.fr/las/
Tél. (33) 02 23 46 14 70


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 1 comment |
headers
Jean Michel Le Bot | 12 Jan 09:34
Picon

Calcul impôt selon revenu et quotient familial

Bonjour
Éternel débutant avec R, les discussions politiques actuelles sur le 
quotient familial m'ont donné l'idée, en guise d'exercice, de tenter de 
reproduire avec R des graphiques comme ceux-ci :
http://fr.wikipedia.or/wiki/Fonction_affine_par_morceaux#Un_exemple_:_la_courbe_d.E2.80.99imposition
http://fr.wikipedia.org/wiki/Imp%C3%B4t_progressif#Progressivit.C3.A9_par_tranches
J'y suis assez bien parvenu avec le code suivant :
x <- seq(5964,150000,1000) # x est le revenu ; on peut bien sûr aller de 
100 en 100...
# formule pour une part (barêmes 2011):
fx <- ifelse(x >= 5964 & x < 11897, (x*0.055)-327.97,
       ifelse(x >= 11897 & x < 26421,  (x*0.14)-1339.13,
       ifelse(x >= 26421 & x < 70831, (x*0.3)-5566.33,
       ifelse(x >= 70831, (x*0.4)-13357.63, NA))))
taux <- fx/x*100
et comme je voulais superposer sur un seul graphique les courbes 
correspondant à différents quotients familiaux j'ai renouvelé les 
calculs pour différents nombres de parts. Exemple:
# formule pour deux parts:
fx2 <- ifelse(x/2 >= 5964 & x/2 < 11897, (x*0.055)-(327.97*2),
       ifelse(x/2 >= 11897 & x/2 < 26421,  (x*0.14)-(1339.13*2),
       ifelse(x/2 >= 26421 & x/2 < 70831, (x*0.3)-(5566.33*2),
       ifelse(x/2 >= 70831, (x*0.4)-(13357.63*2), NA))))
#etc
ensuite évidemment utilisation de plot() avec par(new=TRUE) ou matplot()
Je peux aussi rassembler mes calculs dans un seul data frame avec 
data.frame(x,fx,fx2...) et comparer en ligne l'impôt pour un même revenu 
en utilisant l'indexation. Ça fonctionne, mais c'est peut-être un peu 
basique.
Question : pouvait-on mieux faire ?

-- 
J e a n  -  M i c h e l   L e  B o t
Maître de conférences de sociologie
Département AES
Chercheur au LAS/LARES (Université Rennes 2)
http://www.uhb.fr/sc_humaines/las/
Tél. (33) 02 23 46 14 70


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 2 comments |
headers
Joel GOMBIN | 22 Dec 01:38
Picon

Lissage spatial de données

Bonjour,


connaissez-vous une fonction sous R permettant de lisser facilement des données spatiales, sous forme de polygones ? J'ai trouvé pas mal de choses sur le krieging de points, mais rien sur des polygones...

Pour que les choses soient plus concrètes, j'ai des données au niveau des 36.000 communes, que je souhaiterais lisser pour dégager les grandes tendances, de sorte que la valeur que prend la variable d'intérêt d'une commune soit une moyenne locale (en configurant la distance à prendre en compte et l'intensité du lissage...). 

Une telle fonction existe-t-elle ? Ou dois-je la coder à la main, à partir de la fonction loess ?

Merci,

--
Joël Gombin

Doctorant en science politique / PhD candidate in political science
CURAPP - Université de Picardie Jules Verne

Tel : +33 (0)6 61 55 22 41



--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 1 comment |
headers
Mehrshad Ghaffary | 17 Dec 15:19
Picon

Question méthodologique

Bonjour,

J'ai une série de variables dépendantes discrètes qui peuvent prendre 9 modalités [0;8]. Je dois effectué des modèles de régressions avec une liste de critères explicatives majoritairement qualitatives.

Quelle méthode de régression je dois choisir sachant que ma variable expliquée a une limite supérieure ?

Merci,

MG



--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 3 comments |
headers
Denis Colombi | 31 Oct 17:16
Picon

Variable pondération dans une régression

Bonjour à tous,

Je travaille sur des données de l'Insee qui sont fournies avec une variable de pondération "poids" : chaque individu se voit attribué une valeur qui permet de reconstituer l'ensemble de la population française (pour indication, les premiers valeurs sont : 4766.865 12381.590 14685.431  9516.050  3925.908  9423.811). Entre autres choses, j'essaye de faire des régressions logistiques sur des variables dichotomiques avec glm(). Mais lorsque j'essaye d'intégrer la pondération à l'aide de l'argument weights, j'obtiens les messages d'erreurs suivants :

Messages d'avis :
1: In eval(expr, envir, enclos) :
  nombre de succès non entier dans un glm binomial !
2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart,  :
  l'algorithme n'a pas convergé
3: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart,  :
  des probabilités ont été ajustées numériquement à 0 ou 1

Je n'arrive pas à bien comprendre d'où vient le problème. Est-ce que quelqu'un aurait une idée et pourrait m'aider ?

Merci par avance !


--
Denis Colombi
http://uneheuredepeine.blogspot.com/


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 0 comments |
headers
Olivier Godechot | 22 Oct 12:04
Picon
Picon
Favicon

Regénérer les étiquettes d'une base de données

Bonjour,

Quand on tronque une base de donnée (sélection d'une partie des 
individus selon un critère quelconque), cela engendre souvent des 
problèmes dans les traitements liés à des variables dont certaines 
modalités sont vides.

Je connais bien la technique pour nettoyer chaque variable en la 
"refactorisant" si on me permet l'expression.

he$cdip_b<-factor(he$cdip_b)

Je me demandais s'il existait une fonction pour l'ensemble d'une table, 
ce qui évite de reprendre une à une les variables.

Amicalement
Olivier

-- 
Olivier Godechot

Centre Maurice Halbwachs - CNRS
Ecole Normale Supérieure
Bureau 16 - Bâtiment B 2ème étage
48 boulevard Jourdan
75014 Paris
01 43 13 62 25

Olivier.Godechot@...
http://olivier.godechot.free.fr



--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 2 comments |
headers
sylvie.celerier | 23 Sep 19:52
Picon

solution trouvée

je viens de poster un message sur l'incompatibilité du package rgrs dans la
version 2.13.1. de R. Or, il suffit de se rendre sur le site CRAN
(http://cran.r-project.org/) et de cliquer sur le lien "contrib" pour trouver
tous ses bonheurs en les téléchargeant et en installant via le menu Packages
de R et par la fonction installer la fonction depuis les packages zip.

je donne cette solution pour éviter une réponse inutile et, peut-être, aider
qu-un qui rencontrerait le même problème.


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 1 comment |
headers
sylvie.celerier | 23 Sep 19:40
Picon

version R 2.13.1

Bonjour,
ma question est sans doute aussi naïve qu'est faible mon expérience de R ...
Voilà J'ai dû réinstaller R (version 2.13.1) et voulant charger rgrs j'obtiens
la réponse suivante :  package 'rgrs' is not available (for R version 2.13.1).
me voilà bloquée .. auriez-vous la gentillesse de m'aider ? merci d'avance


--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 0 comments |
headers
Lallemand Carine | 14 Sep 15:07
Picon

Participez a une expérience en sciences humaines !

Pourriez-vous diffuser ce message sur votre liste ?
D'avance un grand merci,

Carine Lallemand


Bonjour à tous,

Dans le cadre de nos recherches sur les Interactions Homme-Machine, Carine Lallemand et moi-même réalisons actuellement une expérience (dont les objectifs ne peuvent pas encore être dévoilés) sous la forme d'un petit jeu de mémorisation.

Nous vous serons ainsi très reconnaissants de bien vouloir y participer !
L'expérience ne vous prendra pas plus de 5 minutes. Elle est accessible depuis ce lien :
http://www.laxlinux.net/ihm

Si vous êtes intéressés par les résultats de cette étude, n'hésitez pas à me contacter et je vous les transmettrai dès qu'ils seront analysés.

Merci par avance pour votre participation !

Guillaume Gronier
guillaume.gronier-bqDqoxiz9wg@public.gmane.org
Dr. en Psychologie ergonomique
Centre de Recherche Public Henri Tudor

--

-- 
Liste des utilisateurs de R en sciences sociales
Informations / options / désabonnement :
https://listes.cru.fr/sympa/info/r-soc
Permalink | Reply | 0 comments |

Gmane