Wo liegt der Fehler?

Wie rufe ich R-Funktionen auf, wie selektiere ich Daten, ich weiß nicht genau ....

Moderatoren: EDi, jogo

RCarlos

Wo liegt der Fehler?

Beitrag von RCarlos »

Hallo,
ich habe mit der Programmierung siehe unten Probleme. Die Warnmeldung verstehe ich nicht.
Ich versuche aus meiner Tabelle d Spalte 12 und 13 den mean zu bilden. Mit den einzelnen Werten funktioniert es, aber wenn ich beide Werte als Vektor nehme, dann gibt er die Fehlermeldung raus. Was mache ich falsch?

danke
Gruß

Code: Alles auswählen

> de<-d[,12]
> tapply(de,d$Geschlecht,mean)
       F        M 
5.473684 4.671429 
> de<-d[,13]
> tapply(de,d$Geschlecht,mean)
       F        M 
6.605263 5.828571 
> de<-d[,c(12,13)]
> by(de,d$Geschlecht,mean)
d$Geschlecht: F
[1] NA
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 
d$Geschlecht: M
[1] NA
Warnmeldungen:
1: In mean.default(data[x, , drop = FALSE], ...) :
  Argument ist weder numerisch noch boolesch: gebe NA zurück
2: In mean.default(data[x, , drop = FALSE], ...) :
  Argument ist weder numerisch noch boolesch: gebe NA zurück
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Wo liegt der Fehler?

Beitrag von jogo »

Hallo RCarlos,

Du hättest aggregate() nehmen sollen:

Code: Alles auswählen

aggregate( as.matrix(d[12:13]) ~ Geschlecht, data=d, FUN=mean)
Der Fehler in Deinem Code liegt darin, dass nach

Code: Alles auswählen

de <- d[,c(12,13)]
das Objekt de ein Dataframe ist. Die Funktion by() splittet den Dataframe entsprechend der Gruppierung. Dabei entstehen kleinere Dataframes. Dataframes sind Listen und auf Listen kann mean() nicht arbeiten:

Code: Alles auswählen

mean(list(a=1:4, b=12:14))
Gruß, Jörg
p.s.:
Bitte lesen und beherzigen: viewtopic.php?f=20&t=29
RCarlos

Re: Wo liegt der Fehler?

Beitrag von RCarlos »

Vielen Dank.

zum p.s.: ich verstehe nicht wie man das macht...also wie man die Code so in die Anfrage stellt, wie ihr das macht...
RCarlos

Re: Wo liegt der Fehler?

Beitrag von RCarlos »

ich habe das jetzt mit aggregate gemacht und es klappt. danke.
Warum klappt es nicht mit standardabweichung? sd?
Danke
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Wo liegt der Fehler?

Beitrag von jogo »

RCarlos hat geschrieben: Mi Jun 21, 2017 10:28 am ich habe das jetzt mit aggregate gemacht und es klappt. danke.
Warum klappt es nicht mit standardabweichung? sd?
Danke
klappt doch:

Code: Alles auswählen

aggregate(as.matrix(iris[3:4]) ~ Species, data=iris, FUN=sd)
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Wo liegt der Fehler?

Beitrag von jogo »

RCarlos hat geschrieben: Mi Jun 21, 2017 10:25 am Vielen Dank.

zum p.s.: ich verstehe nicht wie man das macht...also wie man die Code so in die Anfrage stellt, wie ihr das macht...
in der Zeile oberhalb des Editierbereiches gibt es die Symbole für die Formatierung:
den Text markieren und einfach mal rumprobieren. Der Knopf für die Code-Tags ist </>
Gruß, Jörg
bigben
Beiträge: 2781
Registriert: Mi Okt 12, 2016 9:09 am

Re: Wo liegt der Fehler?

Beitrag von bigben »

RCarlos hat geschrieben: Mi Jun 21, 2017 10:25 amich verstehe nicht wie man das macht...also wie man die Code so in die Anfrage stellt, wie ihr das macht...
https://bbcodes.de/
---
Programmiere stets so, dass die Maxime Deines Programmierstils Grundlage allgemeiner Gesetzgebung sein könnte
RCarlos

Re: Wo liegt der Fehler?

Beitrag von RCarlos »

Code: Alles auswählen

> de<-aggregate( as.matrix(d[c(12,19,26,33,40)]) ~ Geschlecht, data=d, FUN=mean)
> de
  Geschlecht  preVASR   poVASR  po3VASR  po6VASR po12VASR
1          F 5.666667 3.000000 2.833333 3.083333 3.083333
2          M 3.666667 1.111111 2.000000 1.555556 1.666667
> de<-aggregate( as.matrix(d[c(12,19,26,33,40)]) ~ Geschlecht, data=d, FUN=sd)
Fehler in get(as.character(FUN), mode = "function", envir = envir) : 
  Objekt 'FUN' mit Modus 'function' nicht gefunden
bei mir klappt das mit sd nicht...
Gruß
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Wo liegt der Fehler?

Beitrag von jogo »

ja, das ist interessant. Damit wir den Fehler reproduzieren können, kopiere bitte den Output von

Code: Alles auswählen

dput(d)
in Deine nächste Nachricht.

Gruß, Jörg
RCarlos

Re: Wo liegt der Fehler?

Beitrag von RCarlos »

ok, na dann:

Code: Alles auswählen

> dput(d)
structure(list(ID = c(104591566, 104594930, 104601588, 104604650, 
104610250, 104617123, 104639470, 104639542, 104649718, 104667813, 
104669921, 104671837, 104747442, 104805680, 104873227, 104887325, 
104891586, 104927378, 104947182, 104950485, 104977367, 105004656, 
105007809, 105017529, 105043554, 105047258, 105053037, 105065469, 
105071060, 105079532, 105099215, 105168977, 105179952, 105204086, 
105208120, 105217110, 105229755, 105229852, 105232160, 105232381, 
105232403, 105244843, 105248105, 105281374, 105282044, 105299443, 
105302193, 105303076, 105310161, 105318286, 105318316, 105325002, 
105325053, 105325207, 105342691, 105348428, 105354266, 105356625, 
105371500, 105418442, 105419520, 105423640, 105431910, 105445474, 
105469730, 105491557, 105491646, 105510586, 105533179, 105561040, 
105577850, 105583095, 105592990, 105593031, 105627246, 105628340, 
105636172, 105638973, 105639449, 105644523, 105677189, 105688962, 
105692358, 105705751, 105720637, 105733780, 105754604, 105769890, 
105780460, 105792760, 105808461, 105829493, 105834454, 105848773, 
105859473, 105868847, 105884796, 105983760, 106007586, 106008000, 
106011079, 106011761, 106044791, 106047030, 106061432, 106066175, 
106066221, 106066310, 106083010, 106113351, 106117217, 106123004, 
106136076, 106149976, 106153620, 106154708, 106161631, 106161712, 
106175004, 106200025, 106209537, 106219087, 106220891, 106237395, 
106268983, 106269211, 106275408, 106282145, 106282200, 106282331, 
106283856, 106294607, 106307199, 106312656, 106326037, 106334137, 
106348243, 106360553, 106364710, 106397236, 106402647, 106411441, 
106428379, 106433453, 106446458, 106456089), Alter = c(49, 33, 
81, 36, 74, 46, 39, 51, 70, 53, 75, 49, 42, 22, 63, 45, 79, 64, 
75, 55, 23, 27, 39, 75, 53, 30, 39, 41, 71, 66, 61, 53, 56, 45, 
61, 72, 65, 58, 41, 71, 34, 76, 55, 40, 69, 50, 46, 62, 36, 86, 
72, 59, 74, 75, 73, 71, 75, 50, 30, 70, 51, 47, 44, 71, 48, 74, 
56, 65, 61, 19, 47, 22, 74, 51, 72, 73, 42, 43, 78, 84, 74, 69, 
46, 56, 44, 63, 64, 61, 47, 50, 27, 53, 32, 72, 35, 77, 73, 69, 
43, 72, 47, 63, 35, 51, 73, 51, 80, 73, 54, 55, 64, 66, 35, 55, 
62, 66, 51, 34, 31, 50, 57, 46, 47, 71, 68, 48, 46, 53, 24, 60, 
47, 50, 87, 46, 66, 57, 42, 44, 41, 43, 38, 35, 51, 35, 45, 46
), Geschlecht = structure(c(2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 
1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 
2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 
2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 
2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 
1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 
2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 
1L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 
2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 
1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("F", "M"), class = "factor"), 
    Etage = structure(c(29L, 30L, 14L, 17L, 24L, 30L, 17L, 30L, 
    11L, 21L, 21L, 19L, 17L, 30L, 21L, 30L, 14L, 18L, 5L, 30L, 
    17L, 30L, 30L, 17L, 17L, 30L, 30L, 30L, 30L, 18L, 15L, 19L, 
    23L, 7L, 17L, 17L, 30L, 21L, 8L, 7L, 30L, 34L, 30L, 17L, 
    32L, 30L, 17L, 7L, 30L, 14L, 17L, 30L, 34L, 7L, 4L, 19L, 
    14L, 17L, 17L, 30L, 33L, 17L, 17L, 22L, 12L, 14L, 17L, 14L, 
    17L, 17L, 17L, 30L, 9L, 30L, 8L, 3L, 8L, 17L, 18L, 21L, 14L, 
    21L, 19L, 31L, 17L, 1L, 17L, 7L, 26L, 7L, 30L, 33L, 30L, 
    21L, 30L, 13L, 21L, 6L, 28L, 17L, 30L, 17L, 17L, 30L, 14L, 
    17L, 20L, 14L, 7L, 30L, 21L, 14L, 25L, 18L, 30L, 7L, 21L, 
    30L, 17L, 16L, 17L, 17L, 18L, 21L, 18L, 17L, 30L, 10L, 27L, 
    17L, 21L, 17L, 2L, 7L, 9L, 14L, 30L, 30L, 30L, 30L, 17L, 
    30L, 17L, 30L, 30L, 1L), .Label = c("LW2/3 (BSV)", "LW2/3 (SKS)", 
    "LW2/3 (SYN)", "LW2/3, 3/4 (SKS)", "LW2/3, 3/4, 4/5 (SKS)", 
    "LW2/3, 4/5 (SKS)", "LW3/4 (BSV)", "LW3/4 (ELP)", "LW3/4 (SKS)", 
    "LW3/4 (SYN)", "LW3/4, 4/5 (ELP)", "LW3/4, 4/5 (PSL, PLIF)", 
    "LW3/4, 4/5 (SKS COSMIC)", "LW3/4, 4/5 (SKS)", "LW3/4, 4/5 (SKS+BSV)", 
    "LW3/4, 4/5(BSV)", "LW4/5 (BSV)", "LW4/5 (ELP)", "LW4/5 (PLIF PSL)", 
    "LW4/5 (PSL, PLIF)", "LW4/5 (SKS) ", "LW4/5 (SKS+BSV)", "LW4/5 (SYN)", 
    "LW4/5 (TLIF PSL)", "LW4/5, 5/1 (BSV)", "LW4/5, 5/1 (SKS+BSV)", 
    "LW4/5, LW5/SW1 (BSV)", "LW5/1 (BSV)", "LW5/1 (PLIF SLV)", 
    "LW5/SW1 (BSV)", "LW5/SW1 (ELP)", "LW5/SW1 (ELP, COSMIC)", 
    "LW5/SW1 (PLIF SLV)", "LW5/SW1 (SYN)"), class = "factor"), 
    Diagnose = c(5, 1, 2, 1, 5, 1, 1, 1, 1, 2, 2, 5, 1, 1, 2, 
    1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 5, 3, 1, 
    1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 5, 1, 1, 1, 1, 2, 1, 1, 3, 
    1, 2, 5, 2, 1, 1, 1, 5, 1, 1, 4, 5, 2, 1, 2, 1, 1, 1, 1, 
    2, 1, 1, 3, 1, 1, 1, 2, 2, 2, 5, 1, 1, 1, 1, 1, 4, 1, 1, 
    5, 1, 2, 1, 5, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 5, 2, 1, 1, 
    2, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 
    1, 2, 1, 2, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Eingriff = c(1, 
    2, 2, 2, 1, 2, 2, 2, 3, 2, 2, 1, 2, 2, 2, 2, 2, 3, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 1, 2, 2, 2, 2, 2, 2, 3, 
    2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 
    2, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 3, 
    2, 3, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 3, 2, 
    2, 2, 2, 2, 2, 2, 2, 3, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2), Etagen = c(1, 1, 2, 
    1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 3, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 
    1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 
    1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 
    1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Sequestrektomie = c(3, 2, 
    3, 2, 3, 2, 2, 2, 1, 3, 3, 3, 2, 2, 3, 2, 3, 1, 3, 2, 1, 
    2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 3, 3, 1, 2, 2, 2, 3, 1, 1, 
    1, 3, 2, 2, 2, 2, 2, 1, 2, 3, 1, 1, 3, 2, 3, 3, 3, 1, 2, 
    2, 3, 2, 2, 1, 3, 3, 1, 3, 1, 2, 2, 2, 3, 2, 2, 3, 1, 2, 
    1, 3, 3, 3, 3, 1, 2, 1, 2, 2, 2, 1, 2, 3, 1, 3, 2, 3, 3, 
    3, 2, 2, 2, 2, 2, 2, 3, 2, 3, 3, 1, 2, 3, 3, 2, 1, 2, 2, 
    3, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 3, 2, 1, 3, 1, 3, 2, 3, 
    3, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2), WHOStufe = c(1, 1, 1, 2, 
    1, 1, 1, 2, 2, 1, 0, 0, 1, 1, 2, 1, 1, 1, 0, 1, 1, 1, 2, 
    0, 3, 1, 3, 2, 0, 0, 1, 1, 1, 2, 1, 0, 1, 3, 3, 3, 3, 1, 
    1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 0, 2, 1, 2, 3, 1, 2, 1, 1, 
    1, 1, 2, 1, 3, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 2, 2, 1, 
    0, 0, 0, 3, 1, 1, 2, 3, 1, 3, 1, 3, 1, 0, 0, 1, 1, 0, 0, 
    1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 0, 2, 1, 1, 
    1, 1, 1, 3, 1, 1, 3, 2, 1, 3, 1, 2, 0, 3, 3, 1, 1, 1, 3, 
    1, 1, 1, 2, 2, 1, 1, 1, 1), AD = c(2, 2, 1, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 
    2, 2, 2, 2, 2), preBMI = c(23.8, 38.7, 26.5, 34.7, 25.2, 
    39, 24.4, 26.7, 24.4, 34.7, 26.5, 23.4, 31.2, 31.8, 35.2, 
    28.4, 36, 32.8, 29.6, 20.9, 28.4, 22.1, 26.9, 29.4, 25.7, 
    20.1, 19.3, 29.7, 22.9, 24.1, 32.8, 20.7, 37.1, 20.4, 23.3, 
    30.9, 27.5, 25.6, 29.1, 35.3, 37, 25.2, 22.3, 35.8, 26.9, 
    25.8, 24.9, 28.7, 36.6, 26.9, 25.7, 26.4, 25.1, 23, 32.2, 
    29.4, 37.6, 23.5, 37, 23.4, 28.7, 25.7, 20.7, 25.1, 28.7, 
    34.2, 25.1, 30.1, 30.4, 29.7, 22, 34.4, 32.1, 27.9, 27.4, 
    24.8, 30.8, 24.2, 17.2, 26.2, 26.5, 25.8, 22.9, 30.4, 25.7, 
    27.7, 23.2, 33.2, 30.1, 25.3, 24.9, 25.5, 25.1, 26.8, 26, 
    27.5, 24.8, 30.5, 25.1, 29.6, 24.7, 36.7, 44.8, 22, 24.5, 
    27.6, 29.7, 24.2, 25.8, 28.5, 30.9, 39.1, 33.6, 27.3, 26.4, 
    28.7, 28.3, 26.2, 27.2, 32.7, 33.4, 39.6, 22.5, 29.1, 17.3, 
    28.8, 25.7, 36.7, 27.8, 24.5, 28.1, 24.2, 25.4, 35.8, 29.8, 
    30.1, 22.9, 18.8, 22.9, 31.9, 35.7, 45.9, 28.4, 22.9, 24.5, 
    23.5), preVASR = c(6, 7, 9, 2, 7, 8, 7, 3, 1, 5, 3, 6, 7, 
    4, 4, 5, 7, 7, 6, 3, 1, 6, 5, 6, 8, 2, 9, 5, 3, 6, 3, 8, 
    1, 4, 5, 6, 6, 2, 4, 6, 7, 2, 7, 3, 8, 2, 5, 5, 7, 3, 5, 
    4, 6, 5, 2, 2, 7, 9, 4, 6, 5, 3, 5, 7, 3, 10, 4, 5, 4, 6, 
    4, 4, 9, 2, 7, 6, 5, 3, 6, 6, 7, 5, 2, 6, 7, 9, 7, 9, 5, 
    3, 9, 7, 2, 7, 7, 9, 9, 2, 7, 3, 3, 6, 9, 3, 4, 4, 5, 6, 
    4, 5, 2, 9, 7, 4, 4, 7, 8, 0, 2, 5, 2, 1, 2, 4, 9, 8, 7, 
    6, 4, 1, 6, 8, 6, 3, 1, 3, 10, 8, 1, 2, 5, 5, 3, 2, 3, 8), 
    preVASB = c(6, 8, 9, 6, 7, 5, 6, 4, 8, 2, 8, 7, 7, 3, 8, 
    5, 4, 3, 6, 5, 6, 6, 7, 5, 3, 8, 9, 8, 5, 6, 7, 5, 7, 7, 
    7, 8, 7, 4, 7, 10, 7, 2, 7, 5, 8, 7, 6, 6, 5, 8, 10, 2, 8, 
    9, 2, 8, 7, 9, 7, 9, 6, 3, 7, 7, 4, 5, 6, 7, 2, 7, 6, 8, 
    2, 7, 7, 6, 5, 7, 3, 4, 7, 6, 4, 8, 9, 9, 8, 9, 5, 6, 10, 
    0, 4, 8, 7, 1, 9, 5, 7, 10, 6, 8, 7, 8, 4, 3, 7, 8, 5, 3, 
    6, 3, 7, 8, 7, 7, 8, 6, 2, 8, 8, 5, 4, 5, 9, 5, 4, 6, 6, 
    9, 8, 9, 5, 2, 6, 10, 9, 8, 8, 4, 5, 7, 2, 8, 7, 8), preWeg = structure(c(31L, 
    25L, 9L, 20L, 6L, 34L, 25L, 20L, 23L, 34L, 12L, 5L, 25L, 
    34L, 5L, 13L, 24L, 19L, 11L, 34L, 25L, 25L, 5L, 25L, 14L, 
    9L, 3L, 1L, 34L, 24L, 5L, 6L, 34L, 9L, 11L, 4L, 34L, 24L, 
    4L, 2L, 3L, 30L, 34L, 29L, 12L, 34L, 11L, 2L, 1L, 11L, 4L, 
    34L, 21L, 5L, 20L, 4L, 5L, 24L, 4L, 11L, 6L, 34L, 11L, 25L, 
    20L, 2L, 12L, 12L, 20L, 6L, 25L, 25L, 27L, 34L, 12L, 24L, 
    13L, 9L, 2L, 26L, 17L, 11L, 15L, 24L, 12L, 12L, 4L, 16L, 
    5L, 22L, 34L, 4L, 34L, 30L, 4L, 5L, 5L, 12L, 17L, 34L, 24L, 
    11L, 24L, 12L, 33L, 6L, 24L, 8L, 5L, 28L, 28L, 3L, 6L, 4L, 
    34L, 6L, 4L, 18L, 17L, 5L, 29L, 4L, 34L, 20L, 2L, 5L, 30L, 
    5L, 34L, 34L, 7L, 10L, 2L, 34L, 13L, 25L, 2L, 23L, 4L, 34L, 
    34L, 11L, 34L, 20L, 34L, 32L), .Label = c("", "0", "1", "10", 
    "100", "1000", "12", "15", "150", "2", "20", "200", "2000", 
    "215", "225", "25", "250", "3", "30", "300", "40", "400", 
    "5", "50", "500", "60", "600", "700", "75", "750", "800", 
    "90", "900", "max"), class = "factor"), preODI = c(46, 50, 
    71, 50, 53, 50, 54, 32, 42, 10, 44, 40, 56, 18, 64, 46, 38, 
    38, 51, 30, 50, 40, 44, 34, 94, 54, 70, 44, 32, 49, 48, 48, 
    56, 58, 16, 56, 13, 34, 60, 69, 88, 20, 44, 28, 30, 58, 66, 
    24, 36, 51, 64, 16, 39, 51, 14, 67, 44, 36, 76, 56, 14, 20, 
    48, 43, 16, 73, 44, 36, 32, 48, 72, 54, 38, 42, 46, 54, 40, 
    70, 56, 70, 56, 24, 50, 72, 48, 67, 46, 49, 22, 36, 46, 80, 
    22, 40, 68, 49, 28, 33, 44, 12, 56, 53, 58, 82, 35, 40, 60, 
    37, 40, 44, 32, 65, 36, 54, 22, 36, 78, 53, 22, 54, 64, 62, 
    29, 27, 98, 54, 26, 50, 48, 40, 60, 56, 71, 18, 8, 38, 92, 
    73, 84, 16, 32, 54, 33, 42, 47, 72), preHADSA = c(9, 5, 13, 
    5, 9, 9, 5, 9, 8, 3, 3, 3, 13, 4, 11, 3, 2, 7, 9, 8, 1, 6, 
    12, 10, 10, 7, 15, 5, 2, 15, 9, 16, 7, 10, 1, 3, 4, 2, 17, 
    2, 4, 4, 12, 2, 3, 1, 4, 2, 8, 1, 12, 11, 6, 7, 3, 7, 6, 
    5, 9, 7, 1, 6, 8, 1, 2, 4, 8, 6, 7, 10, 14, 16, 1, 8, 5, 
    4, 1, 6, 3, 4, 2, 0, 11, 9, 5, 3, 8, 5, 5, 5, 8, 15, 5, 16, 
    12, 8, 11, 4, 7, 11, 9, 3, 10, 11, 2, 6, 8, 7, 2, 5, 6, 10, 
    13, 4, 2, 5, 12, 12, 6, 2, 5, 0, 8, 0, 13, 1, 10, 7, 10, 
    6, 10, 9, 13, 3, 7, 4, 2, 13, 11, 6, 5, 4, 11, 11, 8, 17), 
    preHADSD = c(10, 3, 8, 4, 6, 5, 2, 2, 5, 2, 4, 3, 10, 4, 
    10, 1, 6, 6, 10, 6, 1, 6, 3, 6, 5, 4, 4, 8, 0, 9, 7, 10, 
    9, 10, 3, 11, 3, 7, 8, 1, 2, 2, 11, 1, 2, 2, 2, 2, 8, 6, 
    4, 7, 3, 6, 3, 5, 3, 5, 9, 5, 3, 2, 14, 1, 1, 9, 4, 6, 6, 
    11, 11, 5, 1, 5, 7, 7, 0, 7, 3, 9, 3, 2, 8, 8, 4, 3, 6, 2, 
    0, 5, 5, 9, 1, 14, 3, 7, 10, 7, 3, 14, 13, 5, 15, 9, 4, 8, 
    11, 4, 1, 2, 5, 13, 6, 4, 1, 1, 15, 13, 5, 7, 8, 9, 2, 4, 
    16, 3, 4, 3, 7, 8, 11, 5, 7, 3, 4, 2, 1, 11, 11, 2, 3, 1, 
    10, 3, 5, 13), poBMI = c(23.8, 38.7, 26.5, 34.7, 25.2, 39, 
    24.4, 26.7, 24.4, 34.7, 26.5, NA, 31.2, 31.8, 35.2, 28.4, 
    36, NA, 29.6, NA, NA, NA, NA, NA, 25.7, 20.1, NA, 29.7, 22.9, 
    NA, 32.8, 20.7, NA, 20.4, 23.3, 30.9, 27.5, 25.6, 29.1, 35.3, 
    NA, NA, 22.3, 35.8, 26.9, NA, NA, NA, 36.6, 26.9, NA, 26.4, 
    25.1, 23, 32.2, 29.4, 37.6, 23.5, 37, NA, 28.7, 25.7, NA, 
    NA, 28.7, 34.2, 25.1, 30.1, NA, 29.7, 22, 34.4, NA, 27.9, 
    27.4, 24.8, 30.8, NA, 17.2, NA, 26.5, 25.8, 22.9, NA, NA, 
    27.7, NA, NA, NA, NA, 24.9, 25.5, NA, 26.8, 26, 27.5, 24.8, 
    30.5, NA, NA, NA, NA, NA, NA, 24.5, NA, NA, 24.2, NA, NA, 
    NA, NA, 33.6, NA, NA, NA, NA, 26.2, NA, NA, 33.4, NA, NA, 
    NA, NA, 28.8, 25.7, NA, NA, NA, NA, 24.2, 25.4, 35.8, 29.8, 
    30.1, 22.9, NA, 22.9, NA, NA, NA, 28.4, 22.9, NA, NA), poVASR = c(7, 
    5, 8, 2, 3, 4, 2, 0, 0, 3, 2, NA, 3, 4, 2, 3, 7, NA, 4, NA, 
    NA, NA, NA, NA, 3, 0, NA, 2, 0, NA, 1, 6, NA, 2, 2, 2, 2, 
    2, 2, 2, NA, NA, 5, 1, 1, NA, NA, NA, 3, 2, NA, 2, 3, 5, 
    4, 2, 7, 1, 2, NA, 0, 1, NA, NA, 2, 2, 3, 2, NA, 2, 3, 1, 
    NA, 0, 3, 2, 0, NA, 2, NA, 4, 3, 3, NA, NA, 2, NA, NA, NA, 
    NA, 2, 8, NA, 4, 1, 4, 2, 1, NA, NA, NA, NA, NA, NA, 1, NA, 
    NA, 2, NA, NA, NA, NA, 5, NA, NA, NA, NA, 3, NA, NA, 0, NA, 
    NA, NA, NA, 2, 3, NA, NA, NA, NA, 0, 3, 2, 2, 3, 1, NA, 1, 
    NA, NA, NA, 2, 2, NA, NA), poVASB = c(1, 0, 8, 0, 3, 0, 2, 
    1, 1, 2, 4, NA, 3, 2, 3, 3, 6, NA, 4, NA, NA, NA, NA, NA, 
    3, 0, NA, 1, 0, NA, 1, 2, NA, 2, 2, 1, 2, 2, 0, 2, NA, NA, 
    0, 0, 1, NA, NA, NA, 0, 2, NA, 0, 3, 1, 2, 2, 2, 0, 1, NA, 
    0, 0, NA, NA, 0, 2, 2, 0, NA, 0, 5, 1, NA, 0, 3, 2, 0, NA, 
    3, NA, 2, 1, 0, NA, NA, 2, NA, NA, NA, NA, 1, 1, NA, 2, 1, 
    2, 2, 1, NA, NA, NA, NA, NA, NA, 1, NA, NA, 3, NA, NA, NA, 
    NA, 5, NA, NA, NA, NA, 3, NA, NA, 1, NA, NA, NA, NA, 0, 2, 
    NA, NA, NA, NA, 0, 2, 0, 3, 2, 1, NA, 2, NA, NA, NA, 2, 1, 
    NA, NA), poWeg = structure(c(6L, 3L, 3L, 16L, 4L, 13L, 4L, 
    14L, 16L, 16L, 4L, 1L, 8L, 16L, 9L, 7L, 13L, 1L, 3L, 1L, 
    1L, 1L, 1L, 1L, 6L, 16L, 1L, 16L, 16L, 1L, 1L, 14L, 1L, 14L, 
    16L, 1L, 16L, 16L, 16L, 16L, 1L, 1L, 14L, 16L, 16L, 1L, 1L, 
    1L, 10L, 16L, 1L, 8L, 16L, 3L, 16L, 1L, 1L, 14L, 3L, 1L, 
    16L, 7L, 1L, 1L, 16L, 12L, 14L, 16L, 1L, 16L, 11L, 1L, 1L, 
    16L, 16L, 16L, 5L, 1L, 13L, 1L, 1L, 1L, 4L, 1L, 1L, 4L, 1L, 
    1L, 1L, 1L, 16L, 2L, 1L, 1L, 14L, 16L, 4L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 16L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 16L, 1L, 1L, 
    1L, 1L, 15L, 1L, 1L, 14L, 1L, 1L, 1L, 1L, 13L, 2L, 1L, 1L, 
    1L, 1L, 16L, 12L, 1L, 14L, 14L, 16L, 1L, 16L, 1L, 1L, 1L, 
    9L, 16L, 1L, 1L), .Label = c("", "10", "100", "1000", "10000", 
    "200", "2000", "250", "300", "40", "400", "5", "50", "500", 
    "600", "max"), class = "factor"), poODI = c(67, 67, 44, 27, 
    45, 60, 46, 20, 7, 8, 16, NA, 53, 18, 42, 40, 36, NA, 38, 
    NA, NA, NA, NA, NA, 44, 9, NA, 2, 9, NA, 12, 67, NA, 42, 
    6, NA, 32, 16, 22.5, 47, NA, NA, 42, 6, 12, NA, NA, NA, 44, 
    17, NA, 42, 4, 51, 18, 25, 36, 12, 53, NA, 2, 24, NA, NA, 
    20, 76, 40, 10, NA, 11, 54, 30, NA, 24, 24, 8, 18, NA, 42, 
    NA, 30, 29, 7, NA, NA, 36, NA, NA, NA, NA, 15, 82, NA, 15, 
    46, 65, 23, 25, NA, NA, NA, NA, NA, NA, 7, NA, NA, 44, NA, 
    NA, NA, NA, 33, NA, NA, NA, NA, 37, NA, NA, 14, NA, NA, NA, 
    NA, 27, 31, NA, NA, NA, NA, 32, 70, 20, 28, 16, 22, NA, 49, 
    NA, NA, NA, 47, 29, NA, NA), poHADSA = c(7, 1, 12, 0, 6, 
    7, 3, 7, 5, 6, 4, NA, 12, 5, 3, 3, 6, NA, 4, NA, NA, NA, 
    NA, NA, 9, 1, NA, 0, 0, NA, 1, 13, NA, 6, 1, 1, 5, 1, 10, 
    0, NA, NA, NA, 3, 1, NA, NA, NA, 9, 1, NA, 10, 4, 6, 5, 3, 
    7, 3, 8, NA, 0, 1, NA, NA, 3, 2, 7, 2, NA, 6, 13, 4, NA, 
    1, 6, 0, 0, NA, 3, NA, 2, 2, 0, NA, NA, 0, NA, NA, NA, NA, 
    4, 12, NA, 13, 4, 10, 7, 4, NA, NA, NA, NA, NA, NA, 1, NA, 
    NA, 10, NA, NA, NA, NA, 10, NA, NA, NA, NA, 13, NA, NA, 2, 
    NA, NA, NA, NA, 0, 4, NA, NA, NA, NA, 1, 9, 3, 7, 1, 1, NA, 
    5, NA, NA, NA, 7, 3, NA, NA), poHADSD = c(5, 2, 8, 2, 4, 
    3, 2, 0, 4, 2, 1, NA, 7, 4, 2, 1, 8, NA, 3, NA, NA, NA, NA, 
    NA, 6, 0, NA, 1, 2, NA, 1, 8, NA, 7, 1, 7, 1, 3, 2, 1, NA, 
    NA, NA, 0, 3, NA, NA, NA, 3, 3, NA, 9, 1, 2, 3, 4, 5, 3, 
    5, NA, 0, 1, NA, NA, 1, 10, 6, 1, NA, 3, 11, 0, NA, 1, 4, 
    0, 0, NA, 3, NA, 1, 1, 0, NA, NA, 1, NA, NA, NA, NA, 1, 7, 
    NA, 11, 0, 4, 5, 3, NA, NA, NA, NA, NA, NA, 2, NA, NA, 4, 
    NA, NA, NA, NA, 6, NA, NA, NA, NA, 12, NA, NA, 3, NA, NA, 
    NA, NA, 3, 2, NA, NA, NA, NA, 1, 11, 2, 3, 1, 1, NA, 6, NA, 
    NA, NA, 8, 1, NA, NA), po3BMI = c(24.4, 38.7, 25.7, 35.6, 
    24.5, 35.9, 24.6, 27.7, 25.3, NA, NA, NA, 31.4, 29.8, 36.3, 
    27.7, NA, NA, NA, 21.3, NA, NA, 26.3, 27, NA, NA, NA, NA, 
    23.1, 29, NA, 21.8, 37.4, 21.1, 23.3, NA, 29.6, 25.6, 28.4, 
    34.3, NA, 25.2, 23.1, NA, 26.6, 25.8, 24.2, NA, NA, 27.3, 
    25.7, NA, 25.3, NA, 31.8, 29.7, NA, 23.5, NA, NA, 28.1, 25.7, 
    NA, 24.7, 28.4, 34.9, NA, NA, 27.7, 29.7, 22, 34.6, 31.1, 
    27.9, NA, 24.8, NA, NA, 18.6, 25.4, 26.5, 25.8, 23.7, 30.1, 
    22.4, 27.7, 22.7, NA, 29.9, 24.7, NA, 25.2, NA, 29.1, 26.2, 
    NA, 25.2, 30.5, NA, NA, 24.3, NA, NA, 22.4, 26.8, NA, NA, 
    NA, NA, NA, 29, NA, NA, 27.3, NA, 28.1, 27.6, NA, NA, NA, 
    33.3, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), po3VASR = c(5, 
    7, 3, 1, 6, 3, 2, 3, 1, NA, NA, NA, 3, 2, 3, 3, NA, NA, NA, 
    3, NA, NA, 0, 4, NA, NA, NA, NA, 3, 5, 2, 3, 1, 1, 4, NA, 
    2, 4, 0, 1, NA, 5, 1, NA, 2, 0, 2, NA, NA, 2, 3, NA, 2, NA, 
    2, 1, NA, 1, NA, NA, 0, 1, NA, 5, 2, 6, NA, NA, 4, 7, 3, 
    5, 3, 0, NA, 0, NA, NA, 3, 2, 1, 0, 3, 6, 4, 3, 5, NA, 2, 
    4, NA, 4, NA, 4, 4, NA, 9, 1, NA, NA, 3, NA, NA, 1, 2, NA, 
    NA, NA, 5, NA, 4, NA, NA, 0, NA, 4, 7, NA, NA, NA, 1, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po3VASB = c(0, 4, 3, 
    0, 6, 0, 2, 6, 3, NA, NA, NA, 6, 2, 3, 2, NA, NA, NA, 1, 
    NA, NA, 0, 3, NA, NA, NA, NA, 2, 4, 2, 4, 0, 1, 4, NA, 2, 
    4, 1, 2, NA, 4, 2, NA, 2, 1, 4, NA, NA, 2, 1, NA, 2, NA, 
    3, 1, NA, 2, NA, NA, 0, 1, NA, 1, 0, 6, NA, NA, 3, 6, 2, 
    5, 3, 6, NA, 0, NA, NA, 3, 2, 2, 0, 2, 1, 3, 2, 3, NA, 1, 
    5, NA, 1, NA, 9, 2, NA, 0, 2, NA, NA, 3, NA, NA, 0, 1, NA, 
    NA, NA, 4, NA, 4, NA, NA, 2, NA, 4, 8, NA, NA, NA, 3, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po3Weg = structure(c(7L, 
    7L, 15L, 22L, 17L, 22L, 22L, 5L, 22L, 1L, 1L, 1L, 7L, 7L, 
    10L, 22L, 1L, 1L, 1L, 22L, 1L, 1L, 22L, 16L, 1L, 1L, 1L, 
    1L, 22L, 19L, 9L, 21L, 22L, 22L, 22L, 1L, 22L, 4L, 22L, 9L, 
    1L, 6L, 22L, 1L, 22L, 9L, 22L, 1L, 1L, 22L, 3L, 1L, 22L, 
    1L, 12L, 13L, 1L, 16L, 1L, 1L, 22L, 22L, 1L, 18L, 22L, 11L, 
    1L, 1L, 16L, 2L, 22L, 4L, 20L, 22L, 1L, 22L, 1L, 1L, 10L, 
    16L, 9L, 22L, 22L, 8L, 17L, 9L, 15L, 1L, 22L, 21L, 1L, 4L, 
    1L, 16L, 11L, 1L, 7L, 16L, 1L, 1L, 22L, 1L, 1L, 22L, 22L, 
    1L, 1L, 1L, 14L, 1L, 22L, 1L, 1L, 22L, 1L, 7L, 16L, 1L, 1L, 
    1L, 16L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("", 
    "0", "100", "1000", "1500", "1750", "2000", "250", "2500", 
    "300", "3000", "350", "3500", "45", "450", "500", "5000", 
    "650", "75", "800", "900", "max"), class = "factor"), po3ODI = c(33, 
    30, 32, 4, 31, 18, 12, 26, 2, NA, NA, NA, 57, 16, 36, 0, 
    NA, NA, NA, 10, NA, NA, 0, 42, NA, NA, NA, NA, 0, 49, 0, 
    40, 30, 6, 6, NA, 2, 34, 8, 18, NA, 18, 16, NA, 0, 0, 10, 
    NA, NA, 13, 55, NA, 4, NA, 4, 16, NA, 20, NA, NA, 0, 2, NA, 
    36, 14, 60, NA, NA, 32, 56, 18, 18, 27, 20, NA, 2, NA, NA, 
    31, 49, 16, 0, 0, 36, 8, 30, 54, NA, 12, 24, NA, 46, NA, 
    53, 30, NA, 26, 22, NA, NA, 26, NA, NA, 4, 7, NA, NA, NA, 
    22, NA, 0, NA, NA, 6, NA, 22, 50, NA, NA, NA, 23, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA), po3HADSA = c(7, 5, 11, 3, 
    6, 2, 3, 7, 5, NA, NA, NA, 13, 5, 5, 0, NA, NA, NA, 2, NA, 
    NA, 1, 13, NA, NA, NA, NA, 0, 14, 2, 7, 6, 5, 0, NA, 5, 2, 
    11, 1, NA, 2, 5, NA, 2, 1, 2, NA, NA, 1, 11, NA, 4, NA, 3, 
    1, NA, 5, NA, NA, 0, 1, NA, 1, 5, 9, NA, NA, 10, 11, 7, 14, 
    1, 3, NA, 0, NA, NA, 1, 6, 0, 0, 0, 1, 5, 9, 12, NA, 1, 3, 
    NA, 10, NA, 17, 9, NA, 11, 3, NA, NA, 7, NA, NA, 6, 2, NA, 
    NA, NA, 3, NA, 2, NA, NA, 2, NA, 7, 7, NA, NA, NA, 1, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po3HADSD = c(4, 4, 7, 
    2, 6, 1, 1, 2, 2, NA, NA, NA, 9, 4, 2, 0, NA, NA, NA, 1, 
    NA, NA, 0, 3, NA, NA, NA, NA, 1, 14, 1, 6, 7, 2, 3, NA, 1, 
    3, 5, 1, NA, 4, 5, NA, 0, 1, 2, NA, NA, 4, 14, NA, 3, NA, 
    1, 1, NA, 2, NA, NA, 0, 0, NA, 1, 3, 6, NA, NA, 7, 13, 5, 
    8, 2, 2, NA, 1, NA, NA, 3, 11, 2, 1, 0, 7, 3, 8, 8, NA, 1, 
    4, NA, 8, NA, 17, 1, NA, 7, 2, NA, NA, 4, NA, NA, 1, 0, NA, 
    NA, NA, 2, NA, 2, NA, NA, 2, NA, 2, 6, NA, NA, NA, 5, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po6BMI = c(NA, 39.4, 
    26.5, 34.8, 25.4, 35.2, 25.4, 26.4, 25.9, NA, NA, 23.6, NA, 
    NA, 35.2, 26.4, 34.8, NA, 33.8, 22, NA, NA, 26.7, 27, NA, 
    NA, 19.6, NA, 22.8, 27.4, NA, 23.7, 34.6, NA, 23.7, 28, 27.7, 
    NA, 28.8, 34.3, NA, 25.2, NA, NA, 26.9, 26.1, 24.2, NA, NA, 
    26.6, NA, 25, 24.7, 22.2, 31.8, 28.1, NA, NA, NA, NA, 27.5, 
    26, 20.7, 25.1, 28.4, 33, NA, 30.1, 29.1, NA, 21.6, 34.1, 
    30.4, 28.2, NA, 24.8, NA, NA, NA, 26.2, NA, 25.8, NA, 30.4, 
    25.7, 27.4, NA, NA, NA, 25, NA, 25.8, NA, 29.1, NA, 22.7, 
    25.2, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA), po6VASR = c(NA, 3, NA, 4, 5, 3, 
    5, 1, 0, NA, NA, 2, NA, NA, 3, 2, 6, NA, 6, 2, NA, NA, 0, 
    6, NA, NA, 7, NA, 2, 7, NA, 2, 1, NA, 1, 4, 2, NA, 2, 3, 
    NA, 3, NA, NA, 1, 0, 2, NA, NA, 3, NA, 4, 3, 5, 4, 2, NA, 
    NA, NA, NA, 0, 2, 3, 5, 2, 4, NA, 1, 2, NA, 3, 4, 5, 0, NA, 
    2, NA, NA, NA, 2, NA, 0, NA, 5, 2, 2, NA, NA, NA, 1, NA, 
    4, NA, 6, NA, 9, 9, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po6VASB = c(NA, 3, NA, 
    1, 7, 0, 4, 4, 2, NA, NA, 5, NA, NA, 3, 2, 2, NA, 4, 2, NA, 
    NA, 0, 4, NA, NA, 1, NA, 2, 4, NA, 0, 1, NA, 3, 8, 2, NA, 
    1, 3, NA, 3, NA, NA, 1, 4, 5, NA, NA, 2, NA, 4, 3, 6, 2, 
    2, NA, NA, NA, NA, 0, 1, 2, 3, 0, 4, NA, 2, 2, NA, 3, 4, 
    5, 4, NA, 2, NA, NA, NA, 0, NA, 0, NA, 0, 2, 2, NA, NA, NA, 
    5, NA, 1, NA, 9, NA, 3, 9, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), po6Weg = structure(c(1L, 
    16L, 16L, 16L, 4L, 16L, 4L, 9L, 16L, 1L, 1L, 9L, 1L, 1L, 
    12L, 16L, 3L, 1L, 8L, 16L, 1L, 1L, 16L, 3L, 1L, 1L, 16L, 
    1L, 16L, 14L, 1L, 16L, 16L, 1L, 16L, 2L, 16L, 1L, 16L, 13L, 
    1L, 7L, 1L, 1L, 16L, 10L, 8L, 1L, 1L, 9L, 1L, 11L, 16L, 6L, 
    12L, 10L, 1L, 1L, 1L, 1L, 16L, 16L, 16L, 6L, 16L, 2L, 1L, 
    4L, 10L, 1L, 16L, 4L, 4L, 16L, 1L, 4L, 1L, 1L, 1L, 15L, 1L, 
    16L, 1L, 12L, 5L, 10L, 1L, 1L, 1L, 6L, 1L, 10L, 1L, 12L, 
    1L, 6L, 4L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("", "10", "100", 
    "1000", "10000", "1500", "1750", "200", "2000", "3000", "400", 
    "500", "5000", "750", "800", "max"), class = "factor"), p6ODI = c(NA, 
    10, 24, 6, 31, 18, 48, 16, 8, NA, NA, 14, NA, NA, 36, 0, 
    40, NA, 51, 6, NA, NA, 0, 51, NA, NA, 34, NA, 2, 53, NA, 
    24, 0, NA, 8, 28, 4, NA, 14, 13, NA, 14, NA, NA, 0, 0, 14, 
    NA, NA, 20, NA, 52, 10, 47, 2, 20, NA, NA, NA, NA, 0, 8, 
    30, 42, 8, 69, NA, 11, 13, NA, 6, 26, 31, 16, NA, 2, NA, 
    NA, NA, 53, NA, 0, NA, 36, 0, 31, NA, NA, NA, 22, NA, 30, 
    NA, 50, NA, 20, 30, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA), po6HADSA = c(NA, 1, 
    8, 3, 8, 1, 2, 8, 5, NA, NA, 1, NA, NA, 1, 0, 4, NA, 5, 1, 
    NA, NA, 0, 15, NA, NA, 11, NA, 0, 17, NA, 6, 0, NA, 1, 10, 
    5, NA, 10, 0, NA, 1, NA, NA, 2, 0, 1, NA, NA, 1, NA, 13, 
    3, 3, 1, 1, NA, NA, NA, NA, 0, 0, 9, 4, 0, 4, NA, 0, 11, 
    NA, 5, 14, 3, 4, NA, 0, NA, NA, NA, 11, NA, 0, NA, 1, 3, 
    6, NA, NA, NA, 4, NA, 12, NA, 18, NA, 15, 15, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), po6HADSD = c(NA, 1, 6, 1, 5, 2, 2, 3, 4, NA, NA, 1, 
    NA, NA, 1, 0, 8, NA, 9, 1, NA, NA, 0, 7, NA, NA, 5, NA, 1, 
    12, NA, 4, 0, NA, 1, 10, 1, NA, 8, 1, NA, 3, NA, NA, 0, 0, 
    1, NA, NA, 4, NA, 12, 2, 4, 1, 1, NA, NA, NA, NA, 0, 1, 8, 
    1, 1, 11, NA, 0, 6, NA, 2, 10, 2, 2, NA, 0, NA, NA, NA, 13, 
    NA, 0, NA, 7, 2, 5, NA, NA, NA, 4, NA, 8, NA, 16, NA, 13, 
    6, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA), po12BMI = c(NA, 38.7, NA, 34.5, 25.7, 
    35.9, 24.4, 27.3, 26.2, 36.3, NA, 23.9, NA, NA, 37.1, 27.4, 
    32, NA, NA, 21.6, NA, NA, 26.9, 27, NA, NA, NA, NA, 22.5, 
    26, NA, 23.7, 30, 21.1, 23, NA, 28.3, NA, NA, 34.3, NA, 25.2, 
    22.9, NA, 26.6, 26.1, 24.2, NA, NA, 27, NA, 24.3, 23.6, NA, 
    30.4, 30.5, NA, NA, NA, NA, 27.8, 26, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA), po12VASR = c(NA, 1, NA, 2, 6, 4, 5, 3, 1, 
    5, NA, 2, NA, NA, 5, 3, 7, NA, NA, 2, NA, NA, 0, 7, NA, NA, 
    NA, NA, 3, 5, NA, 2, 1, 0, 2, NA, 2, NA, NA, 2, NA, 4, 3, 
    NA, 1, 1, 2, NA, NA, 2, NA, 6, 3, NA, 2, 2, NA, NA, NA, NA, 
    0, 1, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), po12VASB = c(NA, 
    1, NA, 1, 8, 4, 4, 1, 1, 2, NA, 3, NA, NA, 3, 3, 2, NA, NA, 
    2, NA, NA, 0, 7, NA, NA, NA, NA, 6, 3, NA, 2, 1, 0, 2, NA, 
    1, NA, NA, 2, NA, 4, 3, NA, 1, 0, 2, NA, NA, 1, NA, 6, 4, 
    NA, 2, 2, NA, NA, NA, NA, 0, 1, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA), po12Weg = structure(c(1L, 9L, 1L, 9L, 2L, 9L, 4L, 
    5L, 9L, 9L, 1L, 4L, 1L, 1L, 7L, 9L, 2L, 1L, 1L, 9L, 1L, 1L, 
    9L, 7L, 1L, 1L, 1L, 1L, 9L, 6L, 1L, 9L, 9L, 9L, 9L, 1L, 9L, 
    1L, 1L, 4L, 1L, 8L, 9L, 1L, 9L, 9L, 9L, 1L, 1L, 4L, 1L, 3L, 
    9L, 1L, 7L, 9L, 1L, 1L, 1L, 1L, 9L, 9L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L), .Label = c("", "100", "200", "2000", "2500", 
    "300", "500", "750", "max"), class = "factor"), po12ODI = c(NA, 
    8, NA, 4, 36, 30, 42, 12, 4, 4, NA, 16, NA, NA, 40, 2, 53, 
    NA, NA, 4, NA, NA, 0, 64, NA, NA, NA, NA, 16, 56, NA, 16, 
    0, 0, 2, NA, 2, NA, NA, 12, NA, 26, 12, NA, 0, 0, 10, NA, 
    NA, 16, NA, 60, 4, NA, 2, 22, NA, NA, NA, NA, 0, 2, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA), po12HADSA = c(NA, 1, NA, 1, 
    8, 4, 4, 7, 6, 6, NA, 1, NA, NA, 6, 0, 9, NA, NA, 2, NA, 
    NA, 0, 12, NA, NA, NA, NA, 1, 20, NA, 7, 0, 3, 0, NA, 6, 
    NA, NA, 3, NA, 3, 5, NA, 2, 0, 1, NA, NA, 1, NA, 16, 4, NA, 
    1, 1, NA, NA, NA, NA, 0, 2, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA), po12HADSD = c(NA, 1, NA, 0, 5, 4, 1, 3, 2, 1, NA, 1, 
    NA, NA, 4, 0, 7, NA, NA, 2, NA, NA, 0, 9, NA, NA, NA, NA, 
    1, 11, NA, 3, 0, 1, 1, NA, 2, NA, NA, 1, NA, 5, 5, NA, 0, 
    0, 4, NA, NA, 7, NA, 18, 2, NA, 1, 1, NA, NA, NA, NA, 0, 
    0, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA)), .Names = c("ID", 
"Alter", "Geschlecht", "Etage", "Diagnose", "Eingriff", "Etagen", 
"Sequestrektomie", "WHOStufe", "AD", "preBMI", "preVASR", "preVASB", 
"preWeg", "preODI", "preHADSA", "preHADSD", "poBMI", "poVASR", 
"poVASB", "poWeg", "poODI", "poHADSA", "poHADSD", "po3BMI", "po3VASR", 
"po3VASB", "po3Weg", "po3ODI", "po3HADSA", "po3HADSD", "po6BMI", 
"po6VASR", "po6VASB", "po6Weg", "p6ODI", "po6HADSA", "po6HADSD", 
"po12BMI", "po12VASR", "po12VASB", "po12Weg", "po12ODI", "po12HADSA", 
"po12HADSD"), row.names = c("1", "2", "3", "4", "5", "6", "7", 
"8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", 
"19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", 
"30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", 
"41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", 
"52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", 
"63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", 
"74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", 
"85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", 
"96", "97", "98", "99", "100", "101", "102", "103", "104", "105", 
"106", "107", "108", "109", "110", "111", "112", "113", "114", 
"115", "116", "117", "118", "119", "120", "121", "122", "123", 
"124", "125", "126", "127", "128", "129", "130", "131", "132", 
"133", "134", "135", "136", "137", "138", "139", "140", "141", 
"142", "143", "144", "145", "146"), class = "data.frame")
und mit

Code: Alles auswählen

> de<-aggregate( as.matrix(d[c(15,22,29,36,43)]) ~ Geschlecht, data=d, FUN=sd, na.rm=TRUE)
Fehler in get(as.character(FUN), mode = "function", envir = envir) : 
  Objekt 'FUN' mit Modus 'function' nicht gefunden
geht es auch nicht.

danke für deine Mühen
Antworten