Я застрял с функцией dcast; Я пытаюсь создать таблицу сумм для особей многих видов за год.

У меня есть фрейм данных с 3 столбцами: (1) год (фактор), (2) названия видов (фактор) и (3) счет (числовой).

Year    Species Counts
2002    SP1 2
2002    SP1 3
2004    SP1 2
2002    SP2 8
2002    SP2 2
2002    SP3 1
2002    SP3 1
2003    SP3 2
2004    SP3 1

Я пытаюсь получить такую таблицу с суммами:

    2002    2003    2004
SP1 5   0   2
SP2 10  0   0
SP3 2   2   1

Агрегат не делает то, что я хочу. Я используюdcast функционировать так:

dcast( DF, Species ~ Year , sum)

Что бы я ни пытался, суммы не работают с факторами. Я получаю это сообщение об ошибке:

(Error in Summary.factor(integer(0), na.rm = FALSE): sum not meaningful for factors 

Когда я пробую функцию по умолчанию (length), я получаю сумму подсчетов строк вместо сумм отдельных лиц. Когда я пытаюсь установить сумму для работы со своим столбцом «Счет» вместо других факторов, он также не работает, или я получаю то же сообщение об ошибке.

Как я могу получить такой стол с суммами графов?

РЕДАКТИРОВАТЬ:

Функция рассчитывает ложные суммы и генерирует лиц. Здесь я покажу вам сравнение сумм, сделанных с помощью перекрестной таблицы в Excel, а затем в R:

EXCEL                   2003    2004    2005    2006    2007    2008    2009
Anthocharis_cardamines  1   0   2   2   0   1   0
Apatura_ilia            0   0   0   0   1   0   0
Aporia_crataegi         2   3   29  26  6   4   3
Brintesia_circe         126 217 199 303 64  99  55


DCAST                   2003    2004    2005    2006    2007    2008    2009
Anthocharis_cardamines  2   0   4   4   0   2   0
Apatura_ilia            0   0   0   0   2   0   0
Aporia_crataegi         4   6   258 205 25  8   6
Brintesia_circe         883 1334    1050    1770    490 848 354

Вычисленные числа даже не соответствуют суммам линий, которые у меня были вчера. Как работают эти суммы?

РЕДАКТИРОВАТЬ 2:

>dput(head(counts, 10)
structure(list(year = structure(c(16L, 16L, 16L, 16L, 16L, 16L, 
16L, 16L, 16L, 15L), .Label = c("1994", "1995", "1996", "1997", 
"1998", "1999", "2000", "2001", "2002", "2003", "2004", "2005", 
"2006", "2007", "2008", "2009"), class = "factor"), species = structure(c(146L, 
146L, 146L, 146L, 146L, 146L, 146L, 146L, 146L, 146L), .Label = c("Aglais_urticae", 
"Anthocharis_cardamines", "Anthocharis_euphenoides", "Apatura_ilia", 
"Apatura_iris", "Aphantopus_hyperantus", "Aporia_crataegi", "Araschnia_levana", 
"Arethusana_arethusa", "Argynnis_adippe", "Argynnis_aglaja", 
"Argynnis_paphia", "Aricia_agestis", "Boloria_dia", "Boloria_euphrosyne", 
"Boloria_selene", "Brenthis_daphne", "Brenthis_ino", "Brintesia_circe", 
"Callophrys_rubi", "Carcharodus_alceae", "Carcharodus_floccifera", 
"Carcharodus_lavatherae", "Carterocephalus_palaemon", "Celastrina_argiolus", 
"Charaxes_jasius", "Chazara_briseis", "Clossiana_dia", "Coenonympha_arcania", 
"Coenonympha_dorus", "Coenonympha_glycerion", "Coenonympha_oedippus", 
"Coenonympha_pamphilus", "Coenonympha_tullia", "Colias_alfacariensis", 
"Colias_croceus", "Colias_hyale", "Colias_palaeno", "Cupido_alcetas", 
"Cupido_argiades", "Cupido_minimus", "Cupido_osiris", "Diacrisia_sannio", 
"Erebia_aethiops", "Erebia_euryale", "Erebia_ligea", "Erebia_medusa", 
"Erebia_meolans", "Erynnis_tages", "Euchloe_crameri", "Euclidia_glyphica", 
"Euphydryas_aurinia", "Euplagia_quadripunctaria", "Everes_argiades", 
"Fabriciana_adippe", "Glaucopsyche_alcon", "Glaucopsyche_alexis", 
"Glaucopsyche_arion", "Glaucopsyche_melanops", "Glaucopsyche_nausithous", 
"Glaucopsyche_teleius", "Gonepteryx_cleopatra", "Gonepteryx_rhamni", 
"Hamearis_lucina", "Hesperia_comma", "Heteropterus_morpheus", 
"Hipparchia_fidia", "Hipparchia_semele", "Hyles_euphorbiae", 
"Hyponephele_lupinus", "Inachis_io", "Iphiclides_podalirius", 
"Issoria_lathonia", "Lampides_boeticus", "Lasiommata_maera", 
"Lasiommata_megera", "Leptidea_sinapis", "Leptotes_pirithous", 
"Libelloides_coccajus", "Libelloides_longicornis", "Limenitis_camilla", 
"Limenitis_populi", "Limenitis_reducta", "Lopinga_achine", "Lycaena_alciphron", 
"Lycaena_dispar", "Lycaena_helle", "Lycaena_phlaeas", "Lycaena_tityrus", 
"Macroglossum_stellatarum", "Maculinea_arion", "Maniola_jurtina", 
"Melanargia_arge", "Melanargia_galathea", "Melanargia_lachesis", 
"Melanargia_occitanica", "Melitaea_cinxia", "Melitaea_diamina", 
"Melitaea_didyma", "Melitaea_phoebe", "Mesoacidalia_aglaja", 
"Minois_dryas", "Neohipparchia_statilinus", "Neozephyrus_quercus", 
"Nymphalis_antiopa", "Nymphalis_polychloros", "Ochlodes_sylvanus", 
"Ochlodes_venatus", "Palaeochrysophanus_hippothoe", "Papilio_machaon", 
"Pararge_aegeria", "Pieris_napi", "Plebeius_agestis", "Plebeius_argyrognomon", 
"Polygonia_c-album", "Polyommatus_bellargus", "Polyommatus_coridon", 
"Polyommatus_escheri", "Polyommatus_icarus", "Polyommatus_semiargus", 
"Polyommatus_thersites", "Pontia_daplidice", "Pseudopanthera_macularia", 
"Pseudophilotes_baton", "Pseudotergumia_fidia", "Pyrgus_malvae", 
"Pyronia_bathseba", "Pyronia_cecilia", "Pyronia_tithonus", "Quercusia_quercus", 
"Satyrium_acaciae", "Satyrium_esculi", "Satyrium_ilicis", "Satyrium_pruni", 
"Satyrium_spini", "Satyrium_w-album", "Smerinthus_ocellatus", 
"Speyeria_aglaja", "Spialia_sertorius", "Thecla_betulae", "Thymelicus_acteon", 
"Thymelicus_lineola", "Thymelicus_sylvestris", "Vanessa_atalanta", 
"Vanessa_cardui", "Zerynthia_polyxena", "Zygaena_carniolica", 
"Zygaena_ephialtes", "Zygaena_erythrus", "Zygaena_fausta", "Zygaena_filipendulae", 
"Zygaena_hilaris", "Zygaena_loti", "Zygaena_occitanica", "Zygaena_purpuralis", 
"Zygaena_sarpedon", "Zygaena_transalpina"), class = "factor"), 
    Counts = c(2, 2, 2, 2, 2, 17, 52, 2, 2, 17)), .Names = c("year", 
"species", "Counts"), row.names = 5479:5488, class = "data.frame")

> str(counts)
'data.frame':   3161 obs. of  3 variables:
 $ year   : Factor w/ 16 levels "1994","1995",..: 16 16 16 16 16 16 16 16 16 15 ...
 $ species: Factor w/ 157 levels "Aglais_urticae",..: 146 146 146 146 146 146 146 146 146 146 ...
 $ Counts : num  2 2 2 2 2 17 52 2 2 17 ...

Я надеюсь, что это помогает ...