eric
2015-May-14 23:00 UTC
[R-es] comportamiento de data.table al hacer calculos por grupos
Estimada comunidad tengo un problema del que no encuentro datos que me
ayuden mucho en la web.
Estoy haciendo calculos por grupos con data,table. Tengo un archivo
(zp.res) con tres columnas que clasifican los datos (sol, con, dia) y
una columna de datos numericos (media), de la siguiente forma:
sol con dia media
1: con 0 1 -22.6
2: con 0 1 -36.6
3: con 0 1 -35.6
y quiero calcular el promedio de "media" (la col 4) agrupando con las
variables sol,con,dia. Lo hago de la siguiente forma:
med <- zp.res[, mean(media), by="sol,dia,con"]
cuando reviso "med" esta todo bien, se han calculado las medias y el
objeto tiene solo las filas que corresponden a los promedios con sus
respectivas columnas sol,con,dia que los identifican. Pero como notaran
por el codigo, la columna con el resultado no tiene un nombre asi es que
R la bautiza como V1. Quise ponele un nombre y agregue este al codigo de
la siguiente forma:
med <- zp.res[, "ave":=mean(media), by="sol,dia,con"]
el problema es que ahora el objeto med tiene el mismo numero de filas
que zp.res y repite el promedio para cada dato dentro del grupo obtenido
con by="sol,dia,con". Esto no me sirve pues tengo que graficar los
promedios ...
mmmmm puede que mi explicacion sea algo confusa, espero que me entiendan.
Encontre que luego puedo renombrar la columna, pero no lo quiero hacer,
pues pienso que deberia ser estandar poder ponerle el nombre y que se
construya bien el objeto con los promedios inmediatamente. Ademas el
promedio es solo uno de los calculos que debo hacer y los otros tambien
quedan con el nombre V1 en la mismo data.table.
Alguna idea de como hacerlo ?
Adjunto archivo con datos.
Saludos y gracias, eric.
--
Forest Engineer
Master in Environmental and Natural Resource Economics
Ph.D. student in Sciences of Natural Resources at La Frontera University
Member in AguaDeTemu2030, citizen movement for Temuco with green city
standards for living
Nota: Las tildes se han omitido para asegurar compatibilidad con algunos
lectores de correo.
------------ próxima parte ------------
"","sol","con","dia","media","ds","num","med","ave"
"1","con",0,1,-22.6,-22.6,1,-29.37,-29.37
"2","con",0,1,-36.6,-36.6,1,-29.37,-29.37
"3","con",0,1,-35.6,-35.6,1,-29.37,-29.37
"4","con",0,1,-27.1,-27.1,1,-29.37,-29.37
"5","con",0,1,-25.5,-25.5,1,-29.37,-29.37
"6","con",0,1,-28.9,-28.9,1,-29.37,-29.37
"7","con",0,1,-30.8,-30.8,1,-29.37,-29.37
"8","con",0,1,-31.4,-31.4,1,-29.37,-29.37
"9","con",0,1,-29.3,-29.3,1,-29.37,-29.37
"10","con",0,1,-25.9,-25.9,1,-29.37,-29.37
"11","dec",1,1,-25.1,-25.1,1,-30.35,-30.35
"12","dec",1,1,-23.8,-23.8,1,-30.35,-30.35
"13","dec",1,1,-25.2,-25.2,1,-30.35,-30.35
"14","dec",1,1,-24.2,-24.2,1,-30.35,-30.35
"15","dec",1,1,-24.8,-24.8,1,-30.35,-30.35
"16","dec",1,1,-29.3,-29.3,1,-30.35,-30.35
"17","dec",1,1,-35.5,-35.5,1,-30.35,-30.35
"18","dec",1,1,-38.3,-38.3,1,-30.35,-30.35
"19","dec",1,1,-39.3,-39.3,1,-30.35,-30.35
"20","dec",1,1,-38,-38,1,-30.35,-30.35
"21","dec",3,1,-38.8,-38.8,1,-37.15,-37.15
"22","dec",3,1,-36.9,-36.9,1,-37.15,-37.15
"23","dec",3,1,-39.2,-39.2,1,-37.15,-37.15
"24","dec",3,1,-45.8,-45.8,1,-37.15,-37.15
"25","dec",3,1,-31.7,-31.7,1,-37.15,-37.15
"26","dec",3,1,-35,-35,1,-37.15,-37.15
"27","dec",3,1,-36.4,-36.4,1,-37.15,-37.15
"28","dec",3,1,-34.4,-34.4,1,-37.15,-37.15
"29","dec",3,1,-36.8,-36.8,1,-37.15,-37.15
"30","dec",3,1,-36.5,-36.5,1,-37.15,-37.15
"31","dec",5,1,-30.2,-30.2,1,-31.55,-31.55
"32","dec",5,1,-35.6,-35.6,1,-31.55,-31.55
"33","dec",5,1,-38.9,-38.9,1,-31.55,-31.55
"34","dec",5,1,-20.1,-20.1,1,-31.55,-31.55
"35","dec",5,1,-32,-32,1,-31.55,-31.55
"36","dec",5,1,-27.8,-27.8,1,-31.55,-31.55
"37","dec",5,1,-33.2,-33.2,1,-31.55,-31.55
"38","dec",5,1,-34.1,-34.1,1,-31.55,-31.55
"39","dec",5,1,-35.4,-35.4,1,-31.55,-31.55
"40","dec",5,1,-28.2,-28.2,1,-31.55,-31.55
"41","dol",1,1,-22.5,-22.5,1,-24.02,-24.02
"42","dol",1,1,-20.7,-20.7,1,-24.02,-24.02
"43","dol",1,1,-22.6,-22.6,1,-24.02,-24.02
"44","dol",1,1,-21.2,-21.2,1,-24.02,-24.02
"45","dol",1,1,-21.5,-21.5,1,-24.02,-24.02
"46","dol",1,1,-21.5,-21.5,1,-24.02,-24.02
"47","dol",1,1,-21.9,-21.9,1,-24.02,-24.02
"48","dol",1,1,-26.4,-26.4,1,-24.02,-24.02
"49","dol",1,1,-34.8,-34.8,1,-24.02,-24.02
"50","dol",1,1,-27.1,-27.1,1,-24.02,-24.02
"51","dol",3,1,-23.4,-23.4,1,-32.95,-32.95
"52","dol",3,1,-33.6,-33.6,1,-32.95,-32.95
"53","dol",3,1,-30.8,-30.8,1,-32.95,-32.95
"54","dol",3,1,-26.6,-26.6,1,-32.95,-32.95
"55","dol",3,1,-33.3,-33.3,1,-32.95,-32.95
"56","dol",3,1,-31.7,-31.7,1,-32.95,-32.95
"57","dol",3,1,-38.2,-38.2,1,-32.95,-32.95
"58","dol",3,1,-38.3,-38.3,1,-32.95,-32.95
"59","dol",3,1,-38.3,-38.3,1,-32.95,-32.95
"60","dol",3,1,-35.3,-35.3,1,-32.95,-32.95
"61","dol",5,1,-31.8,-31.8,1,-34.95,-34.95
"62","dol",5,1,-39,-39,1,-34.95,-34.95
"63","dol",5,1,-36.4,-36.4,1,-34.95,-34.95
"64","dol",5,1,-34.7,-34.7,1,-34.95,-34.95
"65","dol",5,1,-32.9,-32.9,1,-34.95,-34.95
"66","dol",5,1,-33.6,-33.6,1,-34.95,-34.95
"67","dol",5,1,-38,-38,1,-34.95,-34.95
"68","dol",5,1,-32.7,-32.7,1,-34.95,-34.95
"69","dol",5,1,-36,-36,1,-34.95,-34.95
"70","dol",5,1,-34.4,-34.4,1,-34.95,-34.95
"71","lim",1,1,-34.5,-34.5,1,-31.55,-31.55
"72","lim",1,1,-30.8,-30.8,1,-31.55,-31.55
"73","lim",1,1,-33.3,-33.3,1,-31.55,-31.55
"74","lim",1,1,-29.6,-29.6,1,-31.55,-31.55
"75","lim",1,1,-21.8,-21.8,1,-31.55,-31.55
"76","lim",1,1,-30.6,-30.6,1,-31.55,-31.55
"77","lim",1,1,-34.9,-34.9,1,-31.55,-31.55
"78","lim",1,1,-32.2,-32.2,1,-31.55,-31.55
"79","lim",1,1,-31.8,-31.8,1,-31.55,-31.55
"80","lim",1,1,-36,-36,1,-31.55,-31.55
"81","lim",3,1,-25.1,-25.1,1,-25.63,-25.63
"82","lim",3,1,-23.2,-23.2,1,-25.63,-25.63
"83","lim",3,1,-24.6,-24.6,1,-25.63,-25.63
"84","lim",3,1,-28.2,-28.2,1,-25.63,-25.63
"85","lim",3,1,-28.5,-28.5,1,-25.63,-25.63
"86","lim",3,1,-24.5,-24.5,1,-25.63,-25.63
"87","lim",3,1,-22,-22,1,-25.63,-25.63
"88","lim",3,1,-24.8,-24.8,1,-25.63,-25.63
"89","lim",3,1,-27.3,-27.3,1,-25.63,-25.63
"90","lim",3,1,-28.1,-28.1,1,-25.63,-25.63
"91","lim",5,1,-31.1,-31.1,1,-32.24,-32.24
"92","lim",5,1,-31.9,-31.9,1,-32.24,-32.24
"93","lim",5,1,-32.4,-32.4,1,-32.24,-32.24
"94","lim",5,1,-35.3,-35.3,1,-32.24,-32.24
"95","lim",5,1,-33.8,-33.8,1,-32.24,-32.24
"96","lim",5,1,-23.2,-23.2,1,-32.24,-32.24
"97","lim",5,1,-33.2,-33.2,1,-32.24,-32.24
"98","lim",5,1,-33.9,-33.9,1,-32.24,-32.24
"99","lim",5,1,-34,-34,1,-32.24,-32.24
"100","lim",5,1,-33.6,-33.6,1,-32.24,-32.24
"101","con",0,2,-26.8,-26.8,1,-31.65,-31.65
"102","con",0,2,-27.4,-27.4,1,-31.65,-31.65
"103","con",0,2,-32.5,-32.5,1,-31.65,-31.65
"104","con",0,2,-30.6,-30.6,1,-31.65,-31.65
"105","con",0,2,-30.2,-30.2,1,-31.65,-31.65
"106","con",0,2,-31.6,-31.6,1,-31.65,-31.65
"107","con",0,2,-32.6,-32.6,1,-31.65,-31.65
"108","con",0,2,-32.1,-32.1,1,-31.65,-31.65
"109","con",0,2,-36.1,-36.1,1,-31.65,-31.65
"110","con",0,2,-36.6,-36.6,1,-31.65,-31.65
"111","dec",1,2,-28.3,-28.3,1,-29.47,-29.47
"112","dec",1,2,-37.3,-37.3,1,-29.47,-29.47
"113","dec",1,2,-32.1,-32.1,1,-29.47,-29.47
"114","dec",1,2,-29.6,-29.6,1,-29.47,-29.47
"115","dec",1,2,-30.8,-30.8,1,-29.47,-29.47
"116","dec",1,2,-25.8,-25.8,1,-29.47,-29.47
"117","dec",1,2,-26.9,-26.9,1,-29.47,-29.47
"118","dec",1,2,-29,-29,1,-29.47,-29.47
"119","dec",1,2,-28.1,-28.1,1,-29.47,-29.47
"120","dec",1,2,-26.8,-26.8,1,-29.47,-29.47
"121","dec",3,2,-28,-28,1,-31.58,-31.58
"122","dec",3,2,-32,-32,1,-31.58,-31.58
"123","dec",3,2,-30,-30,1,-31.58,-31.58
"124","dec",3,2,-32.6,-32.6,1,-31.58,-31.58
"125","dec",3,2,-32.4,-32.4,1,-31.58,-31.58
"126","dec",3,2,-28.7,-28.7,1,-31.58,-31.58
"127","dec",3,2,-34.4,-34.4,1,-31.58,-31.58
"128","dec",3,2,-31.1,-31.1,1,-31.58,-31.58
"129","dec",3,2,-32.1,-32.1,1,-31.58,-31.58
"130","dec",3,2,-34.5,-34.5,1,-31.58,-31.58
"131","dec",5,2,-28.9658811115,-28.9658811115,1,-27.81041153711,-27.81041153711
"132","dec",5,2,-29.2296869504,-29.2296869504,1,-27.81041153711,-27.81041153711
"133","dec",5,2,-29.2296869504,-29.2296869504,1,-27.81041153711,-27.81041153711
"134","dec",5,2,-31.1290889905,-31.1290889905,1,-27.81041153711,-27.81041153711
"135","dec",5,2,-31.4456559972,-31.4456559972,1,-27.81041153711,-27.81041153711
"136","dec",5,2,-29.2824481182,-29.2824481182,1,-27.81041153711,-27.81041153711
"137","dec",5,2,-19.9964825888,-19.9964825888,1,-27.81041153711,-27.81041153711
"138","dec",5,2,-20.8934224411,-20.8934224411,1,-27.81041153711,-27.81041153711
"139","dec",5,2,-25.2725993669,-25.2725993669,1,-27.81041153711,-27.81041153711
"140","dec",5,2,-32.6591628561,-32.6591628561,1,-27.81041153711,-27.81041153711
"141","dol",1,2,-16.9,-16.9,1,-19.36,-19.36
"142","dol",1,2,-16.9,-16.9,1,-19.36,-19.36
"143","dol",1,2,-19.8,-19.8,1,-19.36,-19.36
"144","dol",1,2,-16.1,-16.1,1,-19.36,-19.36
"145","dol",1,2,-17.6,-17.6,1,-19.36,-19.36
"146","dol",1,2,-17.4,-17.4,1,-19.36,-19.36
"147","dol",1,2,-24.3,-24.3,1,-19.36,-19.36
"148","dol",1,2,-24.7,-24.7,1,-19.36,-19.36
"149","dol",1,2,-17.5,-17.5,1,-19.36,-19.36
"150","dol",1,2,-22.4,-22.4,1,-19.36,-19.36
"151","dol",3,2,-26.5,-26.5,1,-30.12,-30.12
"152","dol",3,2,-26,-26,1,-30.12,-30.12
"153","dol",3,2,-24.8,-24.8,1,-30.12,-30.12
"154","dol",3,2,-32.4,-32.4,1,-30.12,-30.12
"155","dol",3,2,-33.3,-33.3,1,-30.12,-30.12
"156","dol",3,2,-29.1,-29.1,1,-30.12,-30.12
"157","dol",3,2,-31.4,-31.4,1,-30.12,-30.12
"158","dol",3,2,-30.2,-30.2,1,-30.12,-30.12
"159","dol",3,2,-33.5,-33.5,1,-30.12,-30.12
"160","dol",3,2,-34,-34,1,-30.12,-30.12
"161","dol",5,2,-39.7,-39.7,1,-41.36,-41.36
"162","dol",5,2,-34.8,-34.8,1,-41.36,-41.36
"163","dol",5,2,-43.4,-43.4,1,-41.36,-41.36
"164","dol",5,2,-46.1,-46.1,1,-41.36,-41.36
"165","dol",5,2,-40.5,-40.5,1,-41.36,-41.36
"166","dol",5,2,-39.3,-39.3,1,-41.36,-41.36
"167","dol",5,2,-39.2,-39.2,1,-41.36,-41.36
"168","dol",5,2,-44.4,-44.4,1,-41.36,-41.36
"169","dol",5,2,-46.8,-46.8,1,-41.36,-41.36
"170","dol",5,2,-39.4,-39.4,1,-41.36,-41.36
"171","lim",1,2,-25.1,-25.1,1,-30.6,-30.6
"172","lim",1,2,-27.8,-27.8,1,-30.6,-30.6
"173","lim",1,2,-29.6,-29.6,1,-30.6,-30.6
"174","lim",1,2,-28.6,-28.6,1,-30.6,-30.6
"175","lim",1,2,-26.5,-26.5,1,-30.6,-30.6
"176","lim",1,2,-33.1,-33.1,1,-30.6,-30.6
"177","lim",1,2,-34,-34,1,-30.6,-30.6
"178","lim",1,2,-34.3,-34.3,1,-30.6,-30.6
"179","lim",1,2,-33.9,-33.9,1,-30.6,-30.6
"180","lim",1,2,-33.1,-33.1,1,-30.6,-30.6
"181","lim",3,2,-24.6,-24.6,1,-26.95,-26.95
"182","lim",3,2,-29.1,-29.1,1,-26.95,-26.95
"183","lim",3,2,-31.4,-31.4,1,-26.95,-26.95
"184","lim",3,2,-29.5,-29.5,1,-26.95,-26.95
"185","lim",3,2,-29.8,-29.8,1,-26.95,-26.95
"186","lim",3,2,-24.4,-24.4,1,-26.95,-26.95
"187","lim",3,2,-28.4,-28.4,1,-26.95,-26.95
"188","lim",3,2,-29.5,-29.5,1,-26.95,-26.95
"189","lim",3,2,-20.6,-20.6,1,-26.95,-26.95
"190","lim",3,2,-22.2,-22.2,1,-26.95,-26.95
"191","lim",5,2,-26.9,-26.9,1,-26.04,-26.04
"192","lim",5,2,-26.5,-26.5,1,-26.04,-26.04
"193","lim",5,2,-26.6,-26.6,1,-26.04,-26.04
"194","lim",5,2,-29.1,-29.1,1,-26.04,-26.04
"195","lim",5,2,-26.7,-26.7,1,-26.04,-26.04
"196","lim",5,2,-21.8,-21.8,1,-26.04,-26.04
"197","lim",5,2,-22.1,-22.1,1,-26.04,-26.04
"198","lim",5,2,-27.4,-27.4,1,-26.04,-26.04
"199","lim",5,2,-25.2,-25.2,1,-26.04,-26.04
"200","lim",5,2,-28.1,-28.1,1,-26.04,-26.04
"201","con",0,3,-26.2,-26.2,1,-28.25,-28.25
"202","con",0,3,-29.5,-29.5,1,-28.25,-28.25
"203","con",0,3,-31,-31,1,-28.25,-28.25
"204","con",0,3,-27.5,-27.5,1,-28.25,-28.25
"205","con",0,3,-27.7,-27.7,1,-28.25,-28.25
"206","con",0,3,-27.2,-27.2,1,-28.25,-28.25
"207","con",0,3,-26.9,-26.9,1,-28.25,-28.25
"208","con",0,3,-29.1,-29.1,1,-28.25,-28.25
"209","con",0,3,-27.8,-27.8,1,-28.25,-28.25
"210","con",0,3,-29.6,-29.6,1,-28.25,-28.25
"211","dec",1,3,-18.8,-18.8,1,-23.13,-23.13
"212","dec",1,3,-23.1,-23.1,1,-23.13,-23.13
"213","dec",1,3,-24.7,-24.7,1,-23.13,-23.13
"214","dec",1,3,-21.8,-21.8,1,-23.13,-23.13
"215","dec",1,3,-22.6,-22.6,1,-23.13,-23.13
"216","dec",1,3,-20.3,-20.3,1,-23.13,-23.13
"217","dec",1,3,-19.5,-19.5,1,-23.13,-23.13
"218","dec",1,3,-22.1,-22.1,1,-23.13,-23.13
"219","dec",1,3,-28.9,-28.9,1,-23.13,-23.13
"220","dec",1,3,-29.5,-29.5,1,-23.13,-23.13
"221","dec",3,3,-28.8,-28.8,1,-26.52,-26.52
"222","dec",3,3,-25.4,-25.4,1,-26.52,-26.52
"223","dec",3,3,-25.7,-25.7,1,-26.52,-26.52
"224","dec",3,3,-31.7,-31.7,1,-26.52,-26.52
"225","dec",3,3,-32.8,-32.8,1,-26.52,-26.52
"226","dec",3,3,-22.4,-22.4,1,-26.52,-26.52
"227","dec",3,3,-23.8,-23.8,1,-26.52,-26.52
"228","dec",3,3,-21.6,-21.6,1,-26.52,-26.52
"229","dec",3,3,-26.1,-26.1,1,-26.52,-26.52
"230","dec",3,3,-26.9,-26.9,1,-26.52,-26.52
"231","dec",5,3,-29.2,-29.2,1,-25.65,-25.65
"232","dec",5,3,-27.7,-27.7,1,-25.65,-25.65
"233","dec",5,3,-27.5,-27.5,1,-25.65,-25.65
"234","dec",5,3,-29.3,-29.3,1,-25.65,-25.65
"235","dec",5,3,-31.8,-31.8,1,-25.65,-25.65
"236","dec",5,3,-20.7,-20.7,1,-25.65,-25.65
"237","dec",5,3,-22.9,-22.9,1,-25.65,-25.65
"238","dec",5,3,-26.8,-26.8,1,-25.65,-25.65
"239","dec",5,3,-22.8,-22.8,1,-25.65,-25.65
"240","dec",5,3,-17.8,-17.8,1,-25.65,-25.65
"241","dol",1,3,-22,-22,1,-23.12,-23.12
"242","dol",1,3,-21.6,-21.6,1,-23.12,-23.12
"243","dol",1,3,-23.1,-23.1,1,-23.12,-23.12
"244","dol",1,3,-23.2,-23.2,1,-23.12,-23.12
"245","dol",1,3,-25.5,-25.5,1,-23.12,-23.12
"246","dol",1,3,-23.4,-23.4,1,-23.12,-23.12
"247","dol",1,3,-23.4,-23.4,1,-23.12,-23.12
"248","dol",1,3,-23.3,-23.3,1,-23.12,-23.12
"249","dol",1,3,-22,-22,1,-23.12,-23.12
"250","dol",1,3,-23.7,-23.7,1,-23.12,-23.12
"251","dol",3,3,-24.3,-24.3,1,-26.37,-26.37
"252","dol",3,3,-25.7,-25.7,1,-26.37,-26.37
"253","dol",3,3,-25.3,-25.3,1,-26.37,-26.37
"254","dol",3,3,-28.6,-28.6,1,-26.37,-26.37
"255","dol",3,3,-29.9,-29.9,1,-26.37,-26.37
"256","dol",3,3,-23.7,-23.7,1,-26.37,-26.37
"257","dol",3,3,-25.7,-25.7,1,-26.37,-26.37
"258","dol",3,3,-25.3,-25.3,1,-26.37,-26.37
"259","dol",3,3,-27.6,-27.6,1,-26.37,-26.37
"260","dol",3,3,-27.6,-27.6,1,-26.37,-26.37
"261","dol",5,3,-40.4,-40.4,1,-45.69,-45.69
"262","dol",5,3,-43.4,-43.4,1,-45.69,-45.69
"263","dol",5,3,-43.8,-43.8,1,-45.69,-45.69
"264","dol",5,3,-46.5,-46.5,1,-45.69,-45.69
"265","dol",5,3,-55.8,-55.8,1,-45.69,-45.69
"266","dol",5,3,-45,-45,1,-45.69,-45.69
"267","dol",5,3,-44,-44,1,-45.69,-45.69
"268","dol",5,3,-46.1,-46.1,1,-45.69,-45.69
"269","dol",5,3,-47.6,-47.6,1,-45.69,-45.69
"270","dol",5,3,-44.3,-44.3,1,-45.69,-45.69
"271","lim",1,3,-27.1,-27.1,1,-28.8,-28.8
"272","lim",1,3,-30,-30,1,-28.8,-28.8
"273","lim",1,3,-28.9,-28.9,1,-28.8,-28.8
"274","lim",1,3,-30.3,-30.3,1,-28.8,-28.8
"275","lim",1,3,-32.8,-32.8,1,-28.8,-28.8
"276","lim",1,3,-26.6,-26.6,1,-28.8,-28.8
"277","lim",1,3,-28.7,-28.7,1,-28.8,-28.8
"278","lim",1,3,-27.7,-27.7,1,-28.8,-28.8
"279","lim",1,3,-26.8,-26.8,1,-28.8,-28.8
"280","lim",1,3,-29.1,-29.1,1,-28.8,-28.8
"281","lim",3,3,-23.6,-23.6,1,-25.93,-25.93
"282","lim",3,3,-24.3,-24.3,1,-25.93,-25.93
"283","lim",3,3,-24.4,-24.4,1,-25.93,-25.93
"284","lim",3,3,-27.9,-27.9,1,-25.93,-25.93
"285","lim",3,3,-30.9,-30.9,1,-25.93,-25.93
"286","lim",3,3,-25,-25,1,-25.93,-25.93
"287","lim",3,3,-22.5,-22.5,1,-25.93,-25.93
"288","lim",3,3,-22,-22,1,-25.93,-25.93
"289","lim",3,3,-28,-28,1,-25.93,-25.93
"290","lim",3,3,-30.7,-30.7,1,-25.93,-25.93
"291","lim",5,3,-24.2,-24.2,1,-25.4,-25.4
"292","lim",5,3,-21.7,-21.7,1,-25.4,-25.4
"293","lim",5,3,-22.6,-22.6,1,-25.4,-25.4
"294","lim",5,3,-26.5,-26.5,1,-25.4,-25.4
"295","lim",5,3,-25.9,-25.9,1,-25.4,-25.4
"296","lim",5,3,-22.2,-22.2,1,-25.4,-25.4
"297","lim",5,3,-25.4,-25.4,1,-25.4,-25.4
"298","lim",5,3,-25.6,-25.6,1,-25.4,-25.4
"299","lim",5,3,-29,-29,1,-25.4,-25.4
"300","lim",5,3,-30.9,-30.9,1,-25.4,-25.4
"301","con",0,4,-24.3,-24.3,1,-26.17,-26.17
"302","con",0,4,-27,-27,1,-26.17,-26.17
"303","con",0,4,-25.5,-25.5,1,-26.17,-26.17
"304","con",0,4,-27.1,-27.1,1,-26.17,-26.17
"305","con",0,4,-28.2,-28.2,1,-26.17,-26.17
"306","con",0,4,-25.3,-25.3,1,-26.17,-26.17
"307","con",0,4,-26.5,-26.5,1,-26.17,-26.17
"308","con",0,4,-26.9,-26.9,1,-26.17,-26.17
"309","con",0,4,-25.3,-25.3,1,-26.17,-26.17
"310","con",0,4,-25.6,-25.6,1,-26.17,-26.17
"311","dec",1,4,-22,-22,1,-25.89,-25.89
"312","dec",1,4,-24.9,-24.9,1,-25.89,-25.89
"313","dec",1,4,-24.4,-24.4,1,-25.89,-25.89
"314","dec",1,4,-26.9,-26.9,1,-25.89,-25.89
"315","dec",1,4,-28.4,-28.4,1,-25.89,-25.89
"316","dec",1,4,-25.5,-25.5,1,-25.89,-25.89
"317","dec",1,4,-26.2,-26.2,1,-25.89,-25.89
"318","dec",1,4,-25.2,-25.2,1,-25.89,-25.89
"319","dec",1,4,-26.6,-26.6,1,-25.89,-25.89
"320","dec",1,4,-28.8,-28.8,1,-25.89,-25.89
"321","dec",3,4,-26.5,-26.5,1,-27.42,-27.42
"322","dec",3,4,-25.8,-25.8,1,-27.42,-27.42
"323","dec",3,4,-26.5,-26.5,1,-27.42,-27.42
"324","dec",3,4,-31.9,-31.9,1,-27.42,-27.42
"325","dec",3,4,-31.6,-31.6,1,-27.42,-27.42
"326","dec",3,4,-24.6,-24.6,1,-27.42,-27.42
"327","dec",3,4,-24.6,-24.6,1,-27.42,-27.42
"328","dec",3,4,-26.6,-26.6,1,-27.42,-27.42
"329","dec",3,4,-27.1,-27.1,1,-27.42,-27.42
"330","dec",3,4,-29,-29,1,-27.42,-27.42
"331","dec",5,4,-30.9,-30.9,1,-30.3,-30.3
"332","dec",5,4,-32.1,-32.1,1,-30.3,-30.3
"333","dec",5,4,-36.6,-36.6,1,-30.3,-30.3
"334","dec",5,4,-30.4,-30.4,1,-30.3,-30.3
"335","dec",5,4,-30.3,-30.3,1,-30.3,-30.3
"336","dec",5,4,-25.7,-25.7,1,-30.3,-30.3
"337","dec",5,4,-32.1,-32.1,1,-30.3,-30.3
"338","dec",5,4,-28.3,-28.3,1,-30.3,-30.3
"339","dec",5,4,-29.5,-29.5,1,-30.3,-30.3
"340","dec",5,4,-27.1,-27.1,1,-30.3,-30.3
"341","dol",1,4,-18,-18,1,-21.37,-21.37
"342","dol",1,4,-16.9,-16.9,1,-21.37,-21.37
"343","dol",1,4,-18.4,-18.4,1,-21.37,-21.37
"344","dol",1,4,-25.2,-25.2,1,-21.37,-21.37
"345","dol",1,4,-28.6,-28.6,1,-21.37,-21.37
"346","dol",1,4,-16.3,-16.3,1,-21.37,-21.37
"347","dol",1,4,-23.2,-23.2,1,-21.37,-21.37
"348","dol",1,4,-20.6,-20.6,1,-21.37,-21.37
"349","dol",1,4,-22.2,-22.2,1,-21.37,-21.37
"350","dol",1,4,-24.3,-24.3,1,-21.37,-21.37
"351","dol",3,4,-26.3,-26.3,1,-26.48,-26.48
"352","dol",3,4,-29.5,-29.5,1,-26.48,-26.48
"353","dol",3,4,-28.9,-28.9,1,-26.48,-26.48
"354","dol",3,4,-32.6,-32.6,1,-26.48,-26.48
"355","dol",3,4,-29.4,-29.4,1,-26.48,-26.48
"356","dol",3,4,-22.2,-22.2,1,-26.48,-26.48
"357","dol",3,4,-22.5,-22.5,1,-26.48,-26.48
"358","dol",3,4,-24.1,-24.1,1,-26.48,-26.48
"359","dol",3,4,-24.2,-24.2,1,-26.48,-26.48
"360","dol",3,4,-25.1,-25.1,1,-26.48,-26.48
"361","dol",5,4,-37.3,-37.3,1,-38.67,-38.67
"362","dol",5,4,-39,-39,1,-38.67,-38.67
"363","dol",5,4,-40.7,-40.7,1,-38.67,-38.67
"364","dol",5,4,-38.2,-38.2,1,-38.67,-38.67
"365","dol",5,4,-35.1,-35.1,1,-38.67,-38.67
"366","dol",5,4,-36.6,-36.6,1,-38.67,-38.67
"367","dol",5,4,-37.9,-37.9,1,-38.67,-38.67
"368","dol",5,4,-41.8,-41.8,1,-38.67,-38.67
"369","dol",5,4,-42.9,-42.9,1,-38.67,-38.67
"370","dol",5,4,-37.2,-37.2,1,-38.67,-38.67
"371","lim",1,4,-22,-22,1,-24.3,-24.3
"372","lim",1,4,-21.1,-21.1,1,-24.3,-24.3
"373","lim",1,4,-20.4,-20.4,1,-24.3,-24.3
"374","lim",1,4,-28.2,-28.2,1,-24.3,-24.3
"375","lim",1,4,-29,-29,1,-24.3,-24.3
"376","lim",1,4,-22.1,-22.1,1,-24.3,-24.3
"377","lim",1,4,-22.9,-22.9,1,-24.3,-24.3
"378","lim",1,4,-26.4,-26.4,1,-24.3,-24.3
"379","lim",1,4,-24.7,-24.7,1,-24.3,-24.3
"380","lim",1,4,-26.2,-26.2,1,-24.3,-24.3
"381","lim",3,4,-20.2,-20.2,1,-20.98,-20.98
"382","lim",3,4,-21.5,-21.5,1,-20.98,-20.98
"383","lim",3,4,-24,-24,1,-20.98,-20.98
"384","lim",3,4,-21.8,-21.8,1,-20.98,-20.98
"385","lim",3,4,-22.5,-22.5,1,-20.98,-20.98
"386","lim",3,4,-17.5,-17.5,1,-20.98,-20.98
"387","lim",3,4,-18.7,-18.7,1,-20.98,-20.98
"388","lim",3,4,-19.5,-19.5,1,-20.98,-20.98
"389","lim",3,4,-22.8,-22.8,1,-20.98,-20.98
"390","lim",3,4,-21.3,-21.3,1,-20.98,-20.98
"391","lim",5,4,-23.9,-23.9,1,-24.86,-24.86
"392","lim",5,4,-24.1,-24.1,1,-24.86,-24.86
"393","lim",5,4,-26.9,-26.9,1,-24.86,-24.86
"394","lim",5,4,-26.7,-26.7,1,-24.86,-24.86
"395","lim",5,4,-30.2,-30.2,1,-24.86,-24.86
"396","lim",5,4,-21.5,-21.5,1,-24.86,-24.86
"397","lim",5,4,-22.8,-22.8,1,-24.86,-24.86
"398","lim",5,4,-24.7,-24.7,1,-24.86,-24.86
"399","lim",5,4,-21.7,-21.7,1,-24.86,-24.86
"400","lim",5,4,-26.1,-26.1,1,-24.86,-24.86
"401","con",0,5,-27.6,-27.6,1,-27.94,-27.94
"402","con",0,5,-27.1,-27.1,1,-27.94,-27.94
"403","con",0,5,-26.1,-26.1,1,-27.94,-27.94
"404","con",0,5,-32,-32,1,-27.94,-27.94
"405","con",0,5,-32.7,-32.7,1,-27.94,-27.94
"406","con",0,5,-24.8,-24.8,1,-27.94,-27.94
"407","con",0,5,-28.1,-28.1,1,-27.94,-27.94
"408","con",0,5,-28.4,-28.4,1,-27.94,-27.94
"409","con",0,5,-26.7,-26.7,1,-27.94,-27.94
"410","con",0,5,-25.9,-25.9,1,-27.94,-27.94
"411","dec",1,5,-20.3,-20.3,1,-25.27,-25.27
"412","dec",1,5,-25.5,-25.5,1,-25.27,-25.27
"413","dec",1,5,-23.4,-23.4,1,-25.27,-25.27
"414","dec",1,5,-23.1,-23.1,1,-25.27,-25.27
"415","dec",1,5,-21.4,-21.4,1,-25.27,-25.27
"416","dec",1,5,-29.3,-29.3,1,-25.27,-25.27
"417","dec",1,5,-32.5,-32.5,1,-25.27,-25.27
"418","dec",1,5,-30.6,-30.6,1,-25.27,-25.27
"419","dec",1,5,-23.3,-23.3,1,-25.27,-25.27
"420","dec",1,5,-23.3,-23.3,1,-25.27,-25.27
"421","dec",3,5,-29.9,-29.9,1,-32.07,-32.07
"422","dec",3,5,-28.3,-28.3,1,-32.07,-32.07
"423","dec",3,5,-30.3,-30.3,1,-32.07,-32.07
"424","dec",3,5,-31.3,-31.3,1,-32.07,-32.07
"425","dec",3,5,-33.7,-33.7,1,-32.07,-32.07
"426","dec",3,5,-30.1,-30.1,1,-32.07,-32.07
"427","dec",3,5,-32,-32,1,-32.07,-32.07
"428","dec",3,5,-33.7,-33.7,1,-32.07,-32.07
"429","dec",3,5,-35.1,-35.1,1,-32.07,-32.07
"430","dec",3,5,-36.3,-36.3,1,-32.07,-32.07
"431","dec",5,5,-27.9,-27.9,1,-31,-31
"432","dec",5,5,-29.2,-29.2,1,-31,-31
"433","dec",5,5,-29.2,-29.2,1,-31,-31
"434","dec",5,5,-31.1,-31.1,1,-31,-31
"435","dec",5,5,-32.1,-32.1,1,-31,-31
"436","dec",5,5,-32,-32,1,-31,-31
"437","dec",5,5,-31.4,-31.4,1,-31,-31
"438","dec",5,5,-29.5,-29.5,1,-31,-31
"439","dec",5,5,-33.6,-33.6,1,-31,-31
"440","dec",5,5,-34,-34,1,-31,-31
"441","dol",1,5,-23.1,-23.1,1,-22.76,-22.76
"442","dol",1,5,-23.8,-23.8,1,-22.76,-22.76
"443","dol",1,5,-23.6,-23.6,1,-22.76,-22.76
"444","dol",1,5,-24.8,-24.8,1,-22.76,-22.76
"445","dol",1,5,-24.7,-24.7,1,-22.76,-22.76
"446","dol",1,5,-18.9,-18.9,1,-22.76,-22.76
"447","dol",1,5,-21.6,-21.6,1,-22.76,-22.76
"448","dol",1,5,-19.1,-19.1,1,-22.76,-22.76
"449","dol",1,5,-24.7,-24.7,1,-22.76,-22.76
"450","dol",1,5,-23.3,-23.3,1,-22.76,-22.76
"451","dol",3,5,-26.2,-26.2,1,-25.24,-25.24
"452","dol",3,5,-27.1,-27.1,1,-25.24,-25.24
"453","dol",3,5,-28,-28,1,-25.24,-25.24
"454","dol",3,5,-29.7,-29.7,1,-25.24,-25.24
"455","dol",3,5,-28.2,-28.2,1,-25.24,-25.24
"456","dol",3,5,-22.2,-22.2,1,-25.24,-25.24
"457","dol",3,5,-20.8,-20.8,1,-25.24,-25.24
"458","dol",3,5,-20.2,-20.2,1,-25.24,-25.24
"459","dol",3,5,-23.4,-23.4,1,-25.24,-25.24
"460","dol",3,5,-26.6,-26.6,1,-25.24,-25.24
"461","dol",5,5,-45.4,-45.4,1,-40.57,-40.57
"462","dol",5,5,-47.7,-47.7,1,-40.57,-40.57
"463","dol",5,5,-48.1,-48.1,1,-40.57,-40.57
"464","dol",5,5,-46.6,-46.6,1,-40.57,-40.57
"465","dol",5,5,-40.5,-40.5,1,-40.57,-40.57
"466","dol",5,5,-34.9,-34.9,1,-40.57,-40.57
"467","dol",5,5,-32.7,-32.7,1,-40.57,-40.57
"468","dol",5,5,-33.1,-33.1,1,-40.57,-40.57
"469","dol",5,5,-38.5,-38.5,1,-40.57,-40.57
"470","dol",5,5,-38.2,-38.2,1,-40.57,-40.57
"471","lim",1,5,-22.9,-22.9,1,-22.8,-22.8
"472","lim",1,5,-19.6,-19.6,1,-22.8,-22.8
"473","lim",1,5,-19.9,-19.9,1,-22.8,-22.8
"474","lim",1,5,-26.2,-26.2,1,-22.8,-22.8
"475","lim",1,5,-26.5,-26.5,1,-22.8,-22.8
"476","lim",1,5,-20.5,-20.5,1,-22.8,-22.8
"477","lim",1,5,-23,-23,1,-22.8,-22.8
"478","lim",1,5,-23.6,-23.6,1,-22.8,-22.8
"479","lim",1,5,-21.6,-21.6,1,-22.8,-22.8
"480","lim",1,5,-24.2,-24.2,1,-22.8,-22.8
"481","lim",3,5,-20.1,-20.1,1,-21.95,-21.95
"482","lim",3,5,-21.6,-21.6,1,-21.95,-21.95
"483","lim",3,5,-23,-23,1,-21.95,-21.95
"484","lim",3,5,-20.7,-20.7,1,-21.95,-21.95
"485","lim",3,5,-23.7,-23.7,1,-21.95,-21.95
"486","lim",3,5,-16.7,-16.7,1,-21.95,-21.95
"487","lim",3,5,-23.2,-23.2,1,-21.95,-21.95
"488","lim",3,5,-23.8,-23.8,1,-21.95,-21.95
"489","lim",3,5,-22.7,-22.7,1,-21.95,-21.95
"490","lim",3,5,-24,-24,1,-21.95,-21.95
"491","lim",5,5,-21.6,-21.6,1,-23.83,-23.83
"492","lim",5,5,-22.2,-22.2,1,-23.83,-23.83
"493","lim",5,5,-25.8,-25.8,1,-23.83,-23.83
"494","lim",5,5,-25.2,-25.2,1,-23.83,-23.83
"495","lim",5,5,-24.2,-24.2,1,-23.83,-23.83
"496","lim",5,5,-19.4,-19.4,1,-23.83,-23.83
"497","lim",5,5,-24.2,-24.2,1,-23.83,-23.83
"498","lim",5,5,-24.6,-24.6,1,-23.83,-23.83
"499","lim",5,5,-24.9,-24.9,1,-23.83,-23.83
"500","lim",5,5,-26.2,-26.2,1,-23.83,-23.83
"501","con",0,6,-25.2,-25.2,1,-28.68,-28.68
"502","con",0,6,-33.6,-33.6,1,-28.68,-28.68
"503","con",0,6,-29,-29,1,-28.68,-28.68
"504","con",0,6,-26,-26,1,-28.68,-28.68
"505","con",0,6,-29.6,-29.6,1,-28.68,-28.68
"506","con",0,6,-25.2,-25.2,1,-28.68,-28.68
"507","con",0,6,-33.6,-33.6,1,-28.68,-28.68
"508","con",0,6,-29,-29,1,-28.68,-28.68
"509","con",0,6,-26,-26,1,-28.68,-28.68
"510","con",0,6,-29.6,-29.6,1,-28.68,-28.68
"511","dec",1,6,-23,-23,1,-29.24,-29.24
"512","dec",1,6,-24.6,-24.6,1,-29.24,-29.24
"513","dec",1,6,-26.4,-26.4,1,-29.24,-29.24
"514","dec",1,6,-27.6,-27.6,1,-29.24,-29.24
"515","dec",1,6,-31,-31,1,-29.24,-29.24
"516","dec",1,6,-33.5,-33.5,1,-29.24,-29.24
"517","dec",1,6,-29.3,-29.3,1,-29.24,-29.24
"518","dec",1,6,-32.2,-32.2,1,-29.24,-29.24
"519","dec",1,6,-32,-32,1,-29.24,-29.24
"520","dec",1,6,-32.8,-32.8,1,-29.24,-29.24
"521","dec",3,6,-27.7,-27.7,1,-28.66,-28.66
"522","dec",3,6,-27.5,-27.5,1,-28.66,-28.66
"523","dec",3,6,-29.9,-29.9,1,-28.66,-28.66
"524","dec",3,6,-29.9,-29.9,1,-28.66,-28.66
"525","dec",3,6,-28.4,-28.4,1,-28.66,-28.66
"526","dec",3,6,-27.4,-27.4,1,-28.66,-28.66
"527","dec",3,6,-29.1,-29.1,1,-28.66,-28.66
"528","dec",3,6,-29.9,-29.9,1,-28.66,-28.66
"529","dec",3,6,-27.6,-27.6,1,-28.66,-28.66
"530","dec",3,6,-29.2,-29.2,1,-28.66,-28.66
"531","dec",5,6,-27.4175545408,-27.4175545408,1,-28,-28
"532","dec",5,6,-28.6960933536,-28.6960933536,1,-28,-28
"533","dec",5,6,-30.8269913749,-30.8269913749,1,-28,-28
"534","dec",5,6,-24.7184170472,-24.7184170472,1,-28,-28
"535","dec",5,6,-17.0471841705,-17.0471841705,1,-28,-28
"536","dec",5,6,-30.4008117707,-30.4008117707,1,-28,-28
"537","dec",5,6,-28.1278538813,-28.1278538813,1,-28,-28
"538","dec",5,6,-29.4063926941,-29.4063926941,1,-28,-28
"539","dec",5,6,-31.1111111111,-31.1111111111,1,-28,-28
"540","dec",5,6,-32.2475900558,-32.2475900558,1,-28,-28
"541","dol",1,6,-21,-21,1,-24.55,-24.55
"542","dol",1,6,-23.6,-23.6,1,-24.55,-24.55
"543","dol",1,6,-22.8,-22.8,1,-24.55,-24.55
"544","dol",1,6,-23.9,-23.9,1,-24.55,-24.55
"545","dol",1,6,-23.8,-23.8,1,-24.55,-24.55
"546","dol",1,6,-23.3,-23.3,1,-24.55,-24.55
"547","dol",1,6,-27.8,-27.8,1,-24.55,-24.55
"548","dol",1,6,-29.1,-29.1,1,-24.55,-24.55
"549","dol",1,6,-28,-28,1,-24.55,-24.55
"550","dol",1,6,-22.2,-22.2,1,-24.55,-24.55
"551","dol",3,6,-23.2,-23.2,1,-26.51,-26.51
"552","dol",3,6,-23.6,-23.6,1,-26.51,-26.51
"553","dol",3,6,-24.6,-24.6,1,-26.51,-26.51
"554","dol",3,6,-27.2,-27.2,1,-26.51,-26.51
"555","dol",3,6,-26.6,-26.6,1,-26.51,-26.51
"556","dol",3,6,-22.6,-22.6,1,-26.51,-26.51
"557","dol",3,6,-26.4,-26.4,1,-26.51,-26.51
"558","dol",3,6,-30.2,-30.2,1,-26.51,-26.51
"559","dol",3,6,-29.9,-29.9,1,-26.51,-26.51
"560","dol",3,6,-30.8,-30.8,1,-26.51,-26.51
"561","dol",5,6,-37,-37,1,-44.08,-44.08
"562","dol",5,6,-37.4,-37.4,1,-44.08,-44.08
"563","dol",5,6,-39.9,-39.9,1,-44.08,-44.08
"564","dol",5,6,-42.4,-42.4,1,-44.08,-44.08
"565","dol",5,6,-42.2,-42.2,1,-44.08,-44.08
"566","dol",5,6,-46.1,-46.1,1,-44.08,-44.08
"567","dol",5,6,-48.4,-48.4,1,-44.08,-44.08
"568","dol",5,6,-48.7,-48.7,1,-44.08,-44.08
"569","dol",5,6,-50.9,-50.9,1,-44.08,-44.08
"570","dol",5,6,-47.8,-47.8,1,-44.08,-44.08
"571","lim",1,6,-20.4,-20.4,1,-21.96,-21.96
"572","lim",1,6,-25.9,-25.9,1,-21.96,-21.96
"573","lim",1,6,-25.7,-25.7,1,-21.96,-21.96
"574","lim",1,6,-19.7,-19.7,1,-21.96,-21.96
"575","lim",1,6,-21.1,-21.1,1,-21.96,-21.96
"576","lim",1,6,-21.1,-21.1,1,-21.96,-21.96
"577","lim",1,6,-19.3,-19.3,1,-21.96,-21.96
"578","lim",1,6,-19.8,-19.8,1,-21.96,-21.96
"579","lim",1,6,-22.1,-22.1,1,-21.96,-21.96
"580","lim",1,6,-24.5,-24.5,1,-21.96,-21.96
"581","lim",3,6,-16,-16,1,-20.45,-20.45
"582","lim",3,6,-15.3,-15.3,1,-20.45,-20.45
"583","lim",3,6,-16.1,-16.1,1,-20.45,-20.45
"584","lim",3,6,-18.2,-18.2,1,-20.45,-20.45
"585","lim",3,6,-19.2,-19.2,1,-20.45,-20.45
"586","lim",3,6,-21.4,-21.4,1,-20.45,-20.45
"587","lim",3,6,-21.8,-21.8,1,-20.45,-20.45
"588","lim",3,6,-23.6,-23.6,1,-20.45,-20.45
"589","lim",3,6,-25.1,-25.1,1,-20.45,-20.45
"590","lim",3,6,-27.8,-27.8,1,-20.45,-20.45
"591","lim",5,6,-17.7,-17.7,1,-16.65,-16.65
"592","lim",5,6,-15.6,-15.6,1,-16.65,-16.65
"593","lim",5,6,-17.5,-17.5,1,-16.65,-16.65
"594","lim",5,6,-20.3,-20.3,1,-16.65,-16.65
"595","lim",5,6,-15.9,-15.9,1,-16.65,-16.65
"596","lim",5,6,-15.2,-15.2,1,-16.65,-16.65
"597","lim",5,6,-15.1,-15.1,1,-16.65,-16.65
"598","lim",5,6,-15.8,-15.8,1,-16.65,-16.65
"599","lim",5,6,-16.6,-16.6,1,-16.65,-16.65
"600","lim",5,6,-16.8,-16.8,1,-16.65,-16.65
"601","con",0,7,-24.3,-24.3,1,-31.26,-31.26
"602","con",0,7,-30.3,-30.3,1,-31.26,-31.26
"603","con",0,7,-34.5,-34.5,1,-31.26,-31.26
"604","con",0,7,-28.9,-28.9,1,-31.26,-31.26
"605","con",0,7,-29.4,-29.4,1,-31.26,-31.26
"606","con",0,7,-30.6,-30.6,1,-31.26,-31.26
"607","con",0,7,-37.6,-37.6,1,-31.26,-31.26
"608","con",0,7,-36.9,-36.9,1,-31.26,-31.26
"609","con",0,7,-29.8,-29.8,1,-31.26,-31.26
"610","con",0,7,-30.3,-30.3,1,-31.26,-31.26
"611","dec",1,7,-29,-29,1,-30.76,-30.76
"612","dec",1,7,-30,-30,1,-30.76,-30.76
"613","dec",1,7,-31.3,-31.3,1,-30.76,-30.76
"614","dec",1,7,-23.8,-23.8,1,-30.76,-30.76
"615","dec",1,7,-24.8,-24.8,1,-30.76,-30.76
"616","dec",1,7,-29.2,-29.2,1,-30.76,-30.76
"617","dec",1,7,-30.6,-30.6,1,-30.76,-30.76
"618","dec",1,7,-31.9,-31.9,1,-30.76,-30.76
"619","dec",1,7,-37.9,-37.9,1,-30.76,-30.76
"620","dec",1,7,-39.1,-39.1,1,-30.76,-30.76
"621","dec",3,7,-35.7,-35.7,1,-32.88,-32.88
"622","dec",3,7,-34.8,-34.8,1,-32.88,-32.88
"623","dec",3,7,-35.6,-35.6,1,-32.88,-32.88
"624","dec",3,7,-34.2,-34.2,1,-32.88,-32.88
"625","dec",3,7,-33.4,-33.4,1,-32.88,-32.88
"626","dec",3,7,-28.3,-28.3,1,-32.88,-32.88
"627","dec",3,7,-30.9,-30.9,1,-32.88,-32.88
"628","dec",3,7,-30.8,-30.8,1,-32.88,-32.88
"629","dec",3,7,-31.8,-31.8,1,-32.88,-32.88
"630","dec",3,7,-33.3,-33.3,1,-32.88,-32.88
"631","dec",5,7,-35,-35,1,-37.54,-37.54
"632","dec",5,7,-36.9,-36.9,1,-37.54,-37.54
"633","dec",5,7,-36,-36,1,-37.54,-37.54
"634","dec",5,7,-39,-39,1,-37.54,-37.54
"635","dec",5,7,-37.3,-37.3,1,-37.54,-37.54
"636","dec",5,7,-36,-36,1,-37.54,-37.54
"637","dec",5,7,-36.1,-36.1,1,-37.54,-37.54
"638","dec",5,7,-36.7,-36.7,1,-37.54,-37.54
"639","dec",5,7,-41.1,-41.1,1,-37.54,-37.54
"640","dec",5,7,-41.3,-41.3,1,-37.54,-37.54
"641","dol",1,7,-21.4,-21.4,1,-26.12,-26.12
"642","dol",1,7,-22.1,-22.1,1,-26.12,-26.12
"643","dol",1,7,-22.9,-22.9,1,-26.12,-26.12
"644","dol",1,7,-22.9,-22.9,1,-26.12,-26.12
"645","dol",1,7,-24.4,-24.4,1,-26.12,-26.12
"646","dol",1,7,-26.6,-26.6,1,-26.12,-26.12
"647","dol",1,7,-27.9,-27.9,1,-26.12,-26.12
"648","dol",1,7,-27.8,-27.8,1,-26.12,-26.12
"649","dol",1,7,-32.3,-32.3,1,-26.12,-26.12
"650","dol",1,7,-32.9,-32.9,1,-26.12,-26.12
"651","dol",3,7,-32.4,-32.4,1,-32.75,-32.75
"652","dol",3,7,-32.9,-32.9,1,-32.75,-32.75
"653","dol",3,7,-31.6,-31.6,1,-32.75,-32.75
"654","dol",3,7,-31.2,-31.2,1,-32.75,-32.75
"655","dol",3,7,-32.3,-32.3,1,-32.75,-32.75
"656","dol",3,7,-31.2,-31.2,1,-32.75,-32.75
"657","dol",3,7,-32,-32,1,-32.75,-32.75
"658","dol",3,7,-32.4,-32.4,1,-32.75,-32.75
"659","dol",3,7,-36.2,-36.2,1,-32.75,-32.75
"660","dol",3,7,-35.3,-35.3,1,-32.75,-32.75
"661","dol",5,7,-43.4,-43.4,1,-46.52,-46.52
"662","dol",5,7,-41,-41,1,-46.52,-46.52
"663","dol",5,7,-42.3,-42.3,1,-46.52,-46.52
"664","dol",5,7,-50.7,-50.7,1,-46.52,-46.52
"665","dol",5,7,-46.9,-46.9,1,-46.52,-46.52
"666","dol",5,7,-45.4,-45.4,1,-46.52,-46.52
"667","dol",5,7,-49.8,-49.8,1,-46.52,-46.52
"668","dol",5,7,-47.6,-47.6,1,-46.52,-46.52
"669","dol",5,7,-51.5,-51.5,1,-46.52,-46.52
"670","dol",5,7,-46.6,-46.6,1,-46.52,-46.52
"671","lim",1,7,-21.9,-21.9,1,-24.9,-24.9
"672","lim",1,7,-23.4,-23.4,1,-24.9,-24.9
"673","lim",1,7,-22.6,-22.6,1,-24.9,-24.9
"674","lim",1,7,-24,-24,1,-24.9,-24.9
"675","lim",1,7,-25.2,-25.2,1,-24.9,-24.9
"676","lim",1,7,-26.2,-26.2,1,-24.9,-24.9
"677","lim",1,7,-25.7,-25.7,1,-24.9,-24.9
"678","lim",1,7,-27.5,-27.5,1,-24.9,-24.9
"679","lim",1,7,-26,-26,1,-24.9,-24.9
"680","lim",1,7,-26.5,-26.5,1,-24.9,-24.9
"681","lim",3,7,-22.4,-22.4,1,-22.2,-22.2
"682","lim",3,7,-23.9,-23.9,1,-22.2,-22.2
"683","lim",3,7,-23.2,-23.2,1,-22.2,-22.2
"684","lim",3,7,-23.1,-23.1,1,-22.2,-22.2
"685","lim",3,7,-25,-25,1,-22.2,-22.2
"686","lim",3,7,-19.3,-19.3,1,-22.2,-22.2
"687","lim",3,7,-19.9,-19.9,1,-22.2,-22.2
"688","lim",3,7,-18.8,-18.8,1,-22.2,-22.2
"689","lim",3,7,-22.9,-22.9,1,-22.2,-22.2
"690","lim",3,7,-23.5,-23.5,1,-22.2,-22.2
"691","lim",5,7,-28.2,-28.2,1,-28.38,-28.38
"692","lim",5,7,-29,-29,1,-28.38,-28.38
"693","lim",5,7,-29.4,-29.4,1,-28.38,-28.38
"694","lim",5,7,-27.2,-27.2,1,-28.38,-28.38
"695","lim",5,7,-28.2,-28.2,1,-28.38,-28.38
"696","lim",5,7,-26.4,-26.4,1,-28.38,-28.38
"697","lim",5,7,-27.6,-27.6,1,-28.38,-28.38
"698","lim",5,7,-28.4,-28.4,1,-28.38,-28.38
"699","lim",5,7,-29,-29,1,-28.38,-28.38
"700","lim",5,7,-30.4,-30.4,1,-28.38,-28.38
Freddy Omar López Quintero
2015-May-14 23:03 UTC
[R-es] comportamiento de data.table al hacer calculos por grupos
Holap. Reproduje tu ejemplo y ciertamente las filas se duplican innecesariamente. No supe cómo corregir este comportamiento. Si no se pudiera arreglar esto, por favor considera la librería sqldf y obtendrás el resultado correcto:>sqldf("select sol, dia, con, avg(media) as ave from dd group by sol, dia,con") sol dia con ave 1 con 1 0 -29.37000 2 con 2 0 -31.65000 3 con 3 0 -28.25000 ... donde dd es tu conjunto de datos. Salud. 2015-05-14 20:00 GMT-03:00 eric <ericconchamunoz en gmail.com>:> Estimada comunidad tengo un problema del que no encuentro datos que me > ayuden mucho en la web. > > Estoy haciendo calculos por grupos con data,table. Tengo un archivo > (zp.res) con tres columnas que clasifican los datos (sol, con, dia) y una > columna de datos numericos (media), de la siguiente forma: > > sol con dia media > 1: con 0 1 -22.6 > 2: con 0 1 -36.6 > 3: con 0 1 -35.6 > > y quiero calcular el promedio de "media" (la col 4) agrupando con las > variables sol,con,dia. Lo hago de la siguiente forma: > > med <- zp.res[, mean(media), by="sol,dia,con"] > > cuando reviso "med" esta todo bien, se han calculado las medias y el > objeto tiene solo las filas que corresponden a los promedios con sus > respectivas columnas sol,con,dia que los identifican. Pero como notaran por > el codigo, la columna con el resultado no tiene un nombre asi es que R la > bautiza como V1. Quise ponele un nombre y agregue este al codigo de la > siguiente forma: > > med <- zp.res[, "ave":=mean(media), by="sol,dia,con"] > > el problema es que ahora el objeto med tiene el mismo numero de filas que > zp.res y repite el promedio para cada dato dentro del grupo obtenido con > by="sol,dia,con". Esto no me sirve pues tengo que graficar los promedios ... > > mmmmm puede que mi explicacion sea algo confusa, espero que me entiendan. > > Encontre que luego puedo renombrar la columna, pero no lo quiero hacer, > pues pienso que deberia ser estandar poder ponerle el nombre y que se > construya bien el objeto con los promedios inmediatamente. Ademas el > promedio es solo uno de los calculos que debo hacer y los otros tambien > quedan con el nombre V1 en la mismo data.table. > > Alguna idea de como hacerlo ? > > Adjunto archivo con datos. > > Saludos y gracias, eric. > > > > > > -- > Forest Engineer > Master in Environmental and Natural Resource Economics > Ph.D. student in Sciences of Natural Resources at La Frontera University > Member in AguaDeTemu2030, citizen movement for Temuco with green city > standards for living > > Nota: Las tildes se han omitido para asegurar compatibilidad con algunos > lectores de correo. > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- «No soy aquellas sombras tutelares que honré con versos que no olvida el tiempo.» JL Borges [[alternative HTML version deleted]]
Carlos Ortega
2015-May-14 23:28 UTC
[R-es] comportamiento de data.table al hacer calculos por grupos
Hola, La forma de hacerlo con data.table es esta:> library(data.table) > datIn <- fread("zp.res") > med <- datIn[, .(vmed=mean(media)), by=.(sol,dia,con)] > head(med)sol dia con vmed 1: con 1 0 -29.37 2: dec 1 1 -30.35 3: dec 1 3 -37.15 4: dec 1 5 -31.55 5: dol 1 1 -24.02 6: dol 1 3 -32.95 - La solución que has empleado actualiza el valor de una columna. Y aunque te calcula el valor medio por grupo, lo adjunta al data.table original. - En cambio esta otra forma se calcula únicamente el valor medio por grupo. Y bueno si quieres ordenar para que salga el mismo orden que en sqldf, puedes hacer esto otro:> setorder(med, sol, dia) > head(med)sol dia con vmed 1: con 1 0 -29.37 2: con 2 0 -31.65 3: con 3 0 -28.25 4: con 4 0 -26.17 5: con 5 0 -27.94 6: con 6 0 -28.68 Saludos, Carlos Ortega www.qualiytexcellence.es El 15 de mayo de 2015, 1:00, eric <ericconchamunoz en gmail.com> escribió:> Estimada comunidad tengo un problema del que no encuentro datos que me > ayuden mucho en la web. > > Estoy haciendo calculos por grupos con data,table. Tengo un archivo > (zp.res) con tres columnas que clasifican los datos (sol, con, dia) y una > columna de datos numericos (media), de la siguiente forma: > > sol con dia media > 1: con 0 1 -22.6 > 2: con 0 1 -36.6 > 3: con 0 1 -35.6 > > y quiero calcular el promedio de "media" (la col 4) agrupando con las > variables sol,con,dia. Lo hago de la siguiente forma: > > med <- zp.res[, mean(media), by="sol,dia,con"] > > cuando reviso "med" esta todo bien, se han calculado las medias y el > objeto tiene solo las filas que corresponden a los promedios con sus > respectivas columnas sol,con,dia que los identifican. Pero como notaran por > el codigo, la columna con el resultado no tiene un nombre asi es que R la > bautiza como V1. Quise ponele un nombre y agregue este al codigo de la > siguiente forma: > > med <- zp.res[, "ave":=mean(media), by="sol,dia,con"] > > el problema es que ahora el objeto med tiene el mismo numero de filas que > zp.res y repite el promedio para cada dato dentro del grupo obtenido con > by="sol,dia,con". Esto no me sirve pues tengo que graficar los promedios ... > > mmmmm puede que mi explicacion sea algo confusa, espero que me entiendan. > > Encontre que luego puedo renombrar la columna, pero no lo quiero hacer, > pues pienso que deberia ser estandar poder ponerle el nombre y que se > construya bien el objeto con los promedios inmediatamente. Ademas el > promedio es solo uno de los calculos que debo hacer y los otros tambien > quedan con el nombre V1 en la mismo data.table. > > Alguna idea de como hacerlo ? > > Adjunto archivo con datos. > > Saludos y gracias, eric. > > > > > > -- > Forest Engineer > Master in Environmental and Natural Resource Economics > Ph.D. student in Sciences of Natural Resources at La Frontera University > Member in AguaDeTemu2030, citizen movement for Temuco with green city > standards for living > > Nota: Las tildes se han omitido para asegurar compatibilidad con algunos > lectores de correo. > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
eric
2015-May-15 01:00 UTC
[R-es] comportamiento de data.table al hacer calculos por grupos
Muchas gracias Freddy y Carlos ... estuve intentando con .() y con list(), para calcular la media y el error estandar al mismo tiempo en dos columnas, pero me arrojaba un error que no supe interpretar. Ahora ya funciona como sugiere Carlos. Muchas gracias de nuevo. Saludos, Eric. On 14/05/15 19:28, Carlos Ortega wrote:> Hola, > > La forma de hacerlo con data.table es esta: > > > library(data.table) > > datIn <- fread("zp.res") > > med <- datIn[, .(vmed=mean(media)), by=.(sol,dia,con)] > > head(med) > sol dia con vmed > 1: con 1 0 -29.37 > 2: dec 1 1 -30.35 > 3: dec 1 3 -37.15 > 4: dec 1 5 -31.55 > 5: dol 1 1 -24.02 > 6: dol 1 3 -32.95 > > * La solución que has empleado actualiza el valor de una columna. Y > aunque te calcula el valor medio por grupo, lo adjunta al data.table > original. > * En cambio esta otra forma se calcula únicamente el valor medio por > grupo. > > Y bueno si quieres ordenar para que salga el mismo orden que en sqldf, > puedes hacer esto otro: > > > setorder(med, sol, dia) > > head(med) > sol dia con vmed > 1: con 1 0 -29.37 > 2: con 2 0 -31.65 > 3: con 3 0 -28.25 > 4: con 4 0 -26.17 > 5: con 5 0 -27.94 > 6: con 6 0 -28.68 > > > Saludos, > Carlos Ortega > www.qualiytexcellence.es <http://www.qualiytexcellence.es> > > El 15 de mayo de 2015, 1:00, eric <ericconchamunoz en gmail.com > <mailto:ericconchamunoz en gmail.com>> escribió: > > Estimada comunidad tengo un problema del que no encuentro datos que > me ayuden mucho en la web. > > Estoy haciendo calculos por grupos con data,table. Tengo un archivo > (zp.res) con tres columnas que clasifican los datos (sol, con, dia) > y una columna de datos numericos (media), de la siguiente forma: > > sol con dia media > 1: con 0 1 -22.6 > 2: con 0 1 -36.6 > 3: con 0 1 -35.6 > > y quiero calcular el promedio de "media" (la col 4) agrupando con > las variables sol,con,dia. Lo hago de la siguiente forma: > > med <- zp.res[, mean(media), by="sol,dia,con"] > > cuando reviso "med" esta todo bien, se han calculado las medias y el > objeto tiene solo las filas que corresponden a los promedios con sus > respectivas columnas sol,con,dia que los identifican. Pero como > notaran por el codigo, la columna con el resultado no tiene un > nombre asi es que R la bautiza como V1. Quise ponele un nombre y > agregue este al codigo de la siguiente forma: > > med <- zp.res[, "ave":=mean(media), by="sol,dia,con"] > > el problema es que ahora el objeto med tiene el mismo numero de > filas que zp.res y repite el promedio para cada dato dentro del > grupo obtenido con by="sol,dia,con". Esto no me sirve pues tengo que > graficar los promedios ... > > mmmmm puede que mi explicacion sea algo confusa, espero que me > entiendan. > > Encontre que luego puedo renombrar la columna, pero no lo quiero > hacer, pues pienso que deberia ser estandar poder ponerle el nombre > y que se construya bien el objeto con los promedios inmediatamente. > Ademas el promedio es solo uno de los calculos que debo hacer y los > otros tambien quedan con el nombre V1 en la mismo data.table. > > Alguna idea de como hacerlo ? > > Adjunto archivo con datos. > > Saludos y gracias, eric. > > > > > > -- > Forest Engineer > Master in Environmental and Natural Resource Economics > Ph.D. student in Sciences of Natural Resources at La Frontera University > Member in AguaDeTemu2030, citizen movement for Temuco with green > city standards for living > > Nota: Las tildes se han omitido para asegurar compatibilidad con > algunos lectores de correo. > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org <mailto:R-help-es en r-project.org> > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es <http://www.qualityexcellence.es>-- Forest Engineer Master in Environmental and Natural Resource Economics Ph.D. student in Sciences of Natural Resources at La Frontera University Member in AguaDeTemu2030, citizen movement for Temuco with green city standards for living Nota: Las tildes se han omitido para asegurar compatibilidad con algunos lectores de correo.