[R] problem of tapply function
Massimiliano Tripoli
mtripoli at istat.it
Sun Mar 30 09:55:53 CEST 2014
Thanks Bert
----- Messaggio originale -----
Da: "Bert Gunter" <gunter.berton at gene.com>
A: "Massimiliano Tripoli" <mtripoli at istat.it>
Cc: r-help at r-project.org
Inviato: Sabato, 29 marzo 2014 15:52:59
Oggetto: Re: [R] problem of tapply function
The argument to tapply is FUN, not fun. R is case sensitive.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
"Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom."
H. Gilbert Welch
On Sat, Mar 29, 2014 at 6:51 AM, Massimiliano Tripoli <mtripoli at istat.it> wrote:
>
> Dear All,
> I have a vector x (as a factor, see below) for which I'd like to know frequencies for each
> level of the factor itself. I generally use "tapply function" for that but it gave to me
> an unespected result as you can see and according my knowledge of R language and if my memory did not fail.
>
> x <-
> c(16056,16056,16056,16056,16056,16056,16056,16057,16057,16057,16058,16060
> ,16060,16060,16061,16061,16062,16062,16062,16062,16062,16063,16064,16065
> ,16073,16076,16076,16076,16076,16077,16078,16079,16080,16083,16083,16084
> ,16086,16086,16087,16087,16088,16088,16088,16088,16088,16088,16088,16088
> ,16089,16089,16089,16089,16090,16090,16090,16090,16092,16092,16092,16092
> ,16092,16092,16093,16093,16093,16093,16093,16094,16094,16094,16094,16094
> ,16094,16094,16094,16095,16095,16095,16095,16095,16095,16095,16095,16095
> ,16095,16095,16095,16095,16095,16095,16096,16096,16096,16096,16096,16097
> ,16097,16097,16097,16098,16098,16098,16098,16098,16099,16099,16099,16100
> ,16100,16101,16101,16103,16106,16109,16110,16112,16114,16116,16117,16117
> ,16119,16119,16120,16121,16121,16121,16122,16122,16122,16122,16122,16123
> ,16123,16123,16123,16123,16123,16124,16124,16124,16125,16125,16125,16126
> ,16126,16126,16127,16127,16128,16128,16129,16129,16129,16129,16130,16131
> ,16131,16131,16131,16131,16131,16132,16132,16133,16133,16133,16133,16134
> ,16134,16134,16134,16134,16135,16137,16137,16137,16137,16137,16137,16137
> ,16138,16138,16138,16138,16138,16138,16138,16138,16138,16139,16139,16139
> ,16139,16139,16140,16140,16140,16140,16140,16142,16142,16143,16143,16143
> ,16143,16144,16144,16144,16144,16144,16144,16145,16145,16145,16146,16146
> ,16146,16146,16146,16147,16147,16147,16150,16151,16151,16153,16154,16155
> ,16155,16155,16155,16157,16158)
>
> x <- factor(x)
>
> tapply(x,x,fun=length)
> [1] 1 1 1 1 1 1 1 2 2 2 3 4 4 4 5 5 6 6 6 6 6 7 8 9 10
> [26] 11 11 11 11 12 13 14 15 16 16 17 18 18 19 19 20 20 20 20 20 20 20 20 21 21
> [51] 21 21 22 22 22 22 23 23 23 23 23 23 24 24 24 24 24 25 25 25 25 25 25 25 25
> [76] 26 26 26 26 26 26 26 26 26 26 26 26 26 26 26 27 27 27 27 27 28 28 28 28 29
> [101] 29 29 29 29 30 30 30 31 31 32 32 33 34 35 36 37 38 39 40 40 41 41 42 43 43
> [126] 43 44 44 44 44 44 45 45 45 45 45 45 46 46 46 47 47 47 48 48 48 49 49 50 50
> [151] 51 51 51 51 52 53 53 53 53 53 53 54 54 55 55 55 55 56 56 56 56 56 57 58 58
> [176] 58 58 58 58 58 59 59 59 59 59 59 59 59 59 60 60 60 60 60 61 61 61 61 61 62
> [201] 62 63 63 63 63 64 64 64 64 64 64 65 65 65 66 66 66 66 66 67 67 67 68 69 69
> [226] 70 71 72 72 72 72 73 74
>
> table(x)
> x
> 16056 16057 16058 16060 16061 16062 16063 16064 16065 16073 16076 16077 16078
> 7 3 1 3 2 5 1 1 1 1 4 1 1
> 16079 16080 16083 16084 16086 16087 16088 16089 16090 16092 16093 16094 16095
> 1 1 2 1 2 2 8 4 4 6 5 8 15
> 16096 16097 16098 16099 16100 16101 16103 16106 16109 16110 16112 16114 16116
> 5 4 5 3 2 2 1 1 1 1 1 1 1
> 16117 16119 16120 16121 16122 16123 16124 16125 16126 16127 16128 16129 16130
> 2 2 1 3 5 6 3 3 3 2 2 4 1
> 16131 16132 16133 16134 16135 16137 16138 16139 16140 16142 16143 16144 16145
> 6 2 4 5 1 7 9 5 5 2 4 6 3
> 16146 16147 16150 16151 16153 16154 16155 16157 16158
> 5 3 1 2 1 1 4 1 1
>>
> # By looking tapply function help why the table's output is different from tapply's output?
> # Anyone could help to me, please?
> Massimiliano
> --
> Massimiliano Tripoli
> Collaboratore T.E.R. scado il 31/12/2014
> ISTAT - DCCN - Direzione Centrale della Contabilità Nazionale
> U.O. Contabilità dei flussi di materia del sistema economico - CSA/C
> Via Depretis, 74/B 00184 Roma
> Tel. 06.4673.3132
> E-mail: mtripoli at istat.it
>
>
> Precedentemente:
> Collaboratore T.E.R.
> Dipartimento DICS Dipartimento per i conti nazionali e le statistiche economiche
> Direzione DCCN Direzione Centrale della contabilità nazionale
> Servizio FIP Servizio Statistiche della Finanza Pubblica
> Ufficio FIP/E U.O. Conti trimestrali delle amministrazioni pubbliche e conti della sanità
> Via Depretis, 74/B 00184 Roma
> Tel. 06.4673.3132
> E-mail: mtripoli at istat.it
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Massimiliano Tripoli
Collaboratore T.E.R. scado il 31/12/2014
ISTAT - DCCN - Direzione Centrale della Contabilità Nazionale
U.O. Contabilità dei flussi di materia del sistema economico - CSA/C
Via Depretis, 74/B 00184 Roma
Tel. 06.4673.3132
E-mail: mtripoli at istat.it
More information about the R-help
mailing list