Rubén C. Veliz¹, Casimiro D. Gadanha Jr.² y Javier A. Vásquez-Castro³

¹ M.Sc. en Máquinas Agrícolas, consultor privado, Lima, Perú. rdcollantes@hotmail.com,

² Prof. Dr. Departamento de Engenharia Rural, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, SP 13418-900, Brazil. cdgadanh@esalq.usp.br,

³ Prof. Dr. Departamento de Entomología, Universidad Nacional Agraria La Molina, Av. La Universidad s/n, apartado 456, Lima100, Perú. jaque@lamolina.edu.pe Autor para correspondencia.

Recibido: 8-dic-2008 • Aceptado: 1-jul-2009

Mejoras en la uniformidad, eficiencia de deposición del plaguicida y reducción de los volúmenes de aplicación, son requisitos importantes para el éxito del control químico de plagas. Frecuentemente, la plaga se localiza o prefiere ciertas áreas sobre la planta, requiriéndose así una deposición precisa del plaguicida. La irregularidad en la deposición de plaguicidas ocasionará áreas con sub y sobredosis en los campos tratados, lo que podría provocar una presión de selección importante sobre la plaga, favoreciendo la evolución de la resistencia (Vásquez-Castro 2008). Por otro lado, existe tendencia a reducir el volumen de los caldos en las aplicaciones fitosanitarias, para la disminución de costos y el aumento de la eficiencia de pulverización, lo que es posible mediante el uso de sistemas que proporcionen una deposición y un espectro de gotas uniformes (Voll et al. 2006).

En el sistema de pulverización, la “punta” de pulverización es el componente más importante, pues ella es la responsable del caudal, generación y distribución de gotas que cargan el plaguicida hasta el blanco a ser controlado. Las puntas de pulverización son clasificadas básicamente por la forma de energía utilizada para la producción de gotas (energía gaseosa, centrífuga, térmica, eléctrica, hidráulica y combinados). Entre ellas, las puntas hidráulicas son las más importantes en la aplicación de plaguicidas. En ellas, un líquido bajo presión es forzado a través de una pequeña abertura de tal forma que el líquido se esparce, formando una lámina que posteriormente se desintegra en gotas de diferentes tamaños (Vásquez-Castro et al. 2007). Otro tipo de sistema utilizado principalmente para la pulverización en la aviación agrícola es el atomizador rotativo de discos, el cual utiliza la energía centrífuga para la formación de gotas. Diversos trabajos han demostrado las ventajas de estos atomizadores frente a las puntas hidráulicas entre las que se destacan la homogeneidad en el espectro de gotas y la reducción en el volumen de aplicación, características que conducen a una mayor eficiencia biológica y reducción de los costos de operación (Bals 1978; Derksen y Bode 1986).

En la literatura internacional se encuentran algunos trabajos científicos que demuestran la mayor deposición de plaguicidas cuando se disminuyen los volúmenes de aplicación (Hall et al. 1981; Carpenter et al. 1983; Bukovac 1985; Vanee et al. 1985; Raetano 1996; Ramos et al. 2004). Sin embargo, existen otros reportes demostrando lo contrario (Rigolin y Tersi 1999; Farooq y Salyani 2002). Por otro lado, la eficiencia en la deposición y penetración de las gotas en la copa de los árboles va a depender de la arquitectura de la planta (Travis et al. 1987) y del diámetro de las gotas (Salyani 1988). Los mayores depósitos de plaguicidas ocurren en la parte externa e inferior de la copa de los árboles las cuales son fácilmente alcanzadas por el chorro del pulverizador, al contrario de lo que ocurre con las partes interna y superior de la copa las cuales reciben muy poco o nada del producto químico (Morgan 1981; Palladini 1990; Juste et al. 1990). Esa característica operacional ocasionará problemas de residuo de plaguicidas ya que las frutas más expuestas recibirán sobredosis del producto químico y a pesar de cumplirse con el periodo de carencia ese tiempo será insuficiente para que el plaguicida se degrade a niveles por debajo del Limite Máximo de Residuo (LMR) permitido por la Ley, lo cual coloca en riesgo la salud del consumidor y la economía del productor (Vásquez-Castro 2006).

Las pulverizaciones de plaguicidas son generalmente clasificadas en función del tamaño de la gota con particular referencia al diámetro mediano volumétrico (DMV), definido como el diámetro de la gota que divide el volumen de pulverización en dos partes iguales. Así, la mitad del volumen pulverizado formará gotas con diámetros mayores al DMV y la otra mitad con diámetros menores al DMV (Matthews 2000). La variación entre los diámetros de las gotas de una pulverización determina el espectro de gotas, característica definida como amplitud relativa (AR). Si las gotas presentan una diferencia pequeña entre sus diámetros, la pulverización tendrá un espectro de gotas homogéneo. De lo contrario el espectro será heterogéneo. Otra característica importante que debe ser evaluada en las pulverizaciones es la densidad de gotas (DG), definida como la cantidad de gotas por unidad de superficie (Monteiro 2006). A mayor densidad de gotas mayor eficiencia biológica ya que habrá mayor probabilidad que la plaga entre en contacto con la superficie tratada.

De esa forma, la aplicación de plaguicidas se torna una operación bastante compleja debido a muchos factores que intervienen durante su ejecución. Por lo anterior, el objetivo del presente trabajo fue evaluar la eficiencia de dos sistemas de aplicación de plaguicidas en árboles de cítricos.

Los experimentos se realizaron en un huerto de naranjas de la variedad “Valencia” de 10 años de edad, con árboles distanciados a 6 m entre hileras y 3,5 m entre plantas, y con altura promedio de 4 m. Se evaluaron cuatro árboles para cada sistema de aplicación, los cuales fueron seleccionados por sus semejanzas en la arquitectura de la copa a fin de estandarizar ese factor. Las aplicaciones para ambos sistemas se realizaron con un pulverizador hidroneumático de tipo chorro transportado Arbus 2000, con deflector unilateral de tipo “camarón” con ventilador de 0,85 m de diámetro, el cual genera un flujo de aire de 12,45 m3.s-1 según indicaciones del fabricante, tanque de fibra de vidrio de 2000 L de capacidad, bomba modelo JP-150 de 150 L.min-1, comando convencional y filtro modelo FVS-200. Para accionar el pulverizador se usó un tractor agrícola Massey Ferguson, modelo 5285 con motor Perkins, modelo 4001 de 85 CV de potencia en el motor.

Para el sistema de aplicación equipado con puntas hidráulicas se utilizó solamente el ramal inferior del deflector, provisto de 13 puntas de chorro cónico vacío modelo AD3 con difusor AC 45. La presión del sistema fue de 827 kPa, según recomendación del fabricante. Para el sistema de aplicación equipado con atomizadores rotativos se montó una estructura de fibra de vidrio con soportes para colocar cinco atomizadores rotativos de discos, modelo TT-88B (Fig. 1) con las siguientes características técnicas: cuatro discos de plástico con periferia aserrada y diámetro de 0,088 m accionado por motor eléctrico, el cual recibe energía de la batería del tractor, disco de orifício D-4, velocidad de operación de 1,4 a 8,5 m.s-1, velocidad de rotación de 2000 a 6000 rpm con un peso aproximado de 1,67 kg. El control de las rotaciones de los discos del atomizador se realizó mediante un controlador de rotación marca Motron, modelo VV 8702, el cual recibía energía de la batería del tractor. La presión del sistema fue de 276 kPa.

Para la obtención del caudal real correspondiente para cada sistema de aplicación en las presiones de trabajo indicadas se recolectó durante un minuto el volumen en probeta de 2 L de capacidad. La velocidad de operación para el sistema de aplicación equipado con puntas hidráulicas fue de 0,83 m.s^-1 correspondiente a la marcha del tractor 2ª baja y para los atomizadores rotativos fue de 1,66 m.s^-1 correspondiente a la marcha del tractor 3ª baja.

El volumen de caldo a utilizar en cada sistema de pulverización y la cantidad de insecticida (producto comercial) por hectárea se calcularon a través de las fórmulas 1 y 2 citadas por Ramos y Pio (2003).

Donde: Q es el volumen de aplicación (L.ha^-1), q es el caudal colectado de los atomizadores rotativos o puntas hidráulicas (L.min^-1), V es la velocidad del conjunto tractor-pulverizador (km.h^-1) y e es la distancia entre las hileras de los árboles (m).

Donde: Pc es la cantidad de producto comercial (L), D es la dosis del producto (L.ha^-1), Ct es la capacidad del tanque (L) y Q es el volumen de aplicación (L.ha^-1).

Para el estudio del espectro de gotas se utilizó el método propuesto por Derksen y Gray (1995) con algunas modificaciones. Los papeles hidrosensibles (0,076 m de largo, 0,026 m de ancho) se colocaron en diferentes lugares de la planta. La copa fue imaginariamente dividida en dos sectores, hilera de plantas y calle; cada sector fue dividido en tres tercios, inferior, medio y superior, y cada tercio fue dividido en tres profundidades, externa, media e interna, generando así 32 puntos de muestreo por árbol (Fig. 2). Para el tratamiento de los árboles se utilizó el insecticida Metidation (Supracid 400 CE^®) en la concentración de 0,05 g por 100 L de agua. Después de la pulverización, los papeles hidrosensibles se colectaron y analizaron a través de un sistema de análisis de imagen computarizado, e-Sprinkle EMBRAPA, Sao Paulo, Brasil. Durante los ensayos se registraron temperaturas entre 28 y 30°C, humedades relativas del aire entre 55 y 60% y velocidades de viento entre 2,9 y 4,3 km.h^-1.

Para el análisis de la deposición del insecticida sobre las hojas de los árboles se utilizó el método propuesto por Hoffmann y Salyani (1996) con algunas modificaciones. Cuatro hojas de cada punto de muestreo fueron colectadas y almacenadas en cámara fría (-20 ± 5°C) hasta la extracción y cuantificación de los depósitos. El método analítico fue adaptado a partir de Andersson y Palsheden (1998). Para la extracción de los residuos, se colocaron cuatro muestras circulares de hojas de 25 mm de diámetro en tubos de ensayo de 50 mL. Luego se adicionaron 15 mL de acetato de etilo y 2 g de sulfato de sodio, siendo posteriormente, el conjunto homogenizado en Vortex y Ultrason por espacio de 1 y 10 minutos, respectivamente. Se transfirieron alícuotas de 12 mL del sobrenadante a tubos de ensayo de 15 mL, adicionando inmediatamente 50 μL de dodecano. Los extractos se evaporaron en Turbo- Vap LV, en baño-maría a 30°C, con auxilio de aire circulante, previamente secado en filtro de silica-gel desecante azul. Los residuos se resuspendieron en 2 mL de una mezcla de ciclohexano/acetato de etilo (1:1, v/v), homogenizados en Vortex y Ultra-son, y se filtraron a través de un filtro de membrana, Millipore, FG, poro 0,2 μm montado en jeringa hipodérmica de plástico (5 mL). La limpieza de los extractos se realizó mediante técnica de Cromatografía de Permeación en Gel (GPC), utilizando una mezcla de ciclohexano/acetato de etilo (1:1, v/v). Las muestras se analizaron con Cromatografía en Fase Gaseosa, utilizándose cromatógrafo a gas, Hewlett-Packard, modelo 6890, equipado con detector selectivo de masas. Los residuos se calcularon usando el software ChemStation versión B.01.00. El método analítico utilizado para las hojas de cítricos se validó mediante fortificación de la matriz con los niveles de 0,05, 0,5 y 5,0 mg.kg^-1, con tres repeticiones para cada nivel y las recuperaciones entre 70- 120% se consideraron aceptables.

El caudal promedio del atomizador rotativo de discos fue 2,15 L.min^-1. De esa forma, el caudal total del sistema (cinco atomizadores rotativos) fue 10,75 L.min^-1. La velocidad de operación del conjunto tractor-pulverizador fue 1,67 m.s^-1. Como se utilizó apenas un lado del pulverizador para la aplicación, la distancia entre las hileras de árboles fue divida en dos a fin de calcular el volumen del caldo, cuyo valor según la Fórmula 1 fue de 368 L.ha^-1. De acuerdo con la recomendación del fabricante del insecticida, el producto comercial se utilizó a la concentración de 0,125%. Por lo anterior, se utilizaron 500 L de agua y 0,625 L del producto comercial para preparar el caldo insecticida. Con el objetivo de utilizar la misma dosis del insecticida en ambos sistemas de aplicación se despejó la Fórmula 2, cuyo valor fue de 0,45 L.ha^-1 del producto comercial (180 gr de Metidation por hectárea). El caudal promedio de la punta hidráulica fue 1,7 L.min^-1. De esa forma, el caudal total del sistema (13 puntas hidráulicas) fue 22,1 L.min^-1. La velocidad de operación del conjunto tractor-pulverizador se estimó en 0,83 m.s^-1 (Fórmula 1) con el fin de aplicar el volumen de caldo recomendado por el fabricante del insecticida de 1500 L.ha^-1. Esos resultados muestran que el sistema de aplicación equipado con atomizadores rotativos de disco utilizó apenas un cuarto del volumen de caldo utilizado por los sistemas de aplicación tradicionales (puntas hidráulicas). Debido al comportamiento diferencial de las plagas donde algunas se localizan en la parte externa de la copa de los árboles y otras en la parte interna, y a la gran variabilidad en la arquitectura de las copas, las pulverizaciones generalmente se realizan con volúmenes de caldo elevados en la errónea creencia de que así, aumentará la eficacia del producto químico. Por otro lado, el alto costo de operación de las aplicaciones de plaguicidas asociado al requerimiento de las empresas agrícolas de reducir sus costos de producción se ha creado la imperiosa necesidad de mejorar la eficiencia de los sistemas de aplicación tradicionales.

En relación con el espectro de gotas, se observó que el efecto principal del sistema de aplicación solamente fue significativo para el diámetro mediano volumétrico (Tablas 1, 2 y 3). El sistema de aplicación equipado con puntas hidráulicas presentó un DMV (679 μm) significativamente superior al del sistema equipado con atomizadores rotativos (278 μm). En el proceso de formación de gotas, la energía hidráulica del líquido se transforma en energía cinética de las gotas (Lefebvre 1989). Así, las gotas grandes poseen mayor masa y, por lo tanto adquieren mayor energía cinética. De ese modo, gotas grandes tendrán mayor capacidad para vencer la resistencia del aire al desplazamiento horizontal, llegando a recorrer distancias superiores en relación a las gotas pequeñas (Vásquez- Castro et al. 2007). Por otro lado, las gotas pequeñas tienen mayor capacidad para penetrar la copa de los árboles, llegando a lugares donde las gotas grandes no pueden hacerlo. Sin embargo, las gotas muy pequeñas se evaporan fácilmente o son transportadas por el viento para áreas que no son objeto de tratamiento (deriva), ocasionando pérdidas del plaguicida. En ese sentido, el diámetro de gota más adecuado para el tratamiento de árboles de cítricos parece estar comprendido entre 240 y 340 μm (Salyani 1988) tal como se observó en el sistema de pulverización equipado con atomizadores rotativos de disco.

Los porcentajes de recuperación del insecticida en las hojas de cítricos fueron aceptables (70-120%), validándose el método analítico. En el testigo no se detectó el insecticida Metidation, indicando que las hojas estaban libres de este compuesto. En la Tabla 1 se observa que independientemente del sistema de aplicación utilizado, hubo mayor depósito de Metidation en el sector de la copa expuesta a la calle en comparación con la hilera. Así mismo, los valores de amplitud relativa y densidad de gotas fueron superiores en ese sector. A mayor valor de amplitud relativa, el espectro de gotas será más heterogéneo. Es decir, habrá una mayor variabilidad en el tamaño de gotas que componen el chorro, característica indeseable en las aplicaciones de plaguicidas. En la Tabla 2 se observa que independientemente del sistema de aplicación utilizado, hubo mayor depósito del insecticida en el tercio inferior de la copa frente a los tercios medio y superior. A su vez, los depósitos fueron superiores en la parte externa de la copa frente a las partes media e interna (Tabla 3). Esos resultados demuestran que la parte externa del tercio inferior de la copa del árbol de naranja es la región que recibe mayor cantidad de plaguicidas durante el tratamiento. Por otro lado, las ramas y hojas del árbol se comportan como un “filtro de gotas”, reteniendo las gotas grandes en la parte externa y dejando pasar las gotas pequeñas al interior de la copa. En la literatura se encuentran diversos trabajos que presentan resultados similares (Morgan 1981; Matuo 1988; Derksen y Gray 1995; Hoffmann y Salyani 1996; Cunningham y Harden 1998; Ramos et al. 2004). Finalmente, a pesar de que ambos sistemas de aplicación no presentaron diferencias significativas en la mayoría de los parámetros evaluados, el sistema de aplicación equipado con atomizadores rotativos de disco utilizó apenas un cuarto del volumen de caldo utilizado por el sistema de aplicación tradicional (puntas hidráulicas), lo cual representa beneficio económico para el productor agrícola. Es obvio que una reducción en el volumen de agua utilizado en las aplicaciones de plaguicidas reducirá el costo de operación y por ende habrá un beneficio económico para el productor. La cuantificación de ese beneficio económico es materia de otra investigación.

A la Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES, Brasil por la concesión de la beca al primer autor y a los Drs. Gilberto Casadei de Baptista y Luiz Roberto Pimentel Trevizan por la orientación durante los análisis cromatográficos.

ANDERSSON, A.; PALSHEDEN, H. 1998. Multi-residue method for the analysis of. pesticides in fruit and vegetables using acetate extraction, GPC clean up and GC determination. En: National Food Administration. Pesticide Analytical in Sweden. Uppsala, p. 9-41.        [ Links ]

BALS, E.J. 1978. The reason for C.D.A. (controlled drop application). Proceedings British Crop Production. Weeds 2: 659-666.        [ Links ]

BUKOVAC, M.J. 1985. Spray application technology: shortcomings and opportunities with especial reference to tree fruits. In: Proceedings of the conference by agricultural research institute, Bethesda. p. 25-38.        [ Links ]

CARPENTER, T.G.; REICHARD, D.L.; KHAN, A.S. 1983. Spray deposition from a row-crop airblast sprayer. Transactions of the ASAE 26(2): 338-342.        [ Links ]

CUNNINGHAM, G.P.; HARDEN, J. 1998. Reducing spray volumes applied to mature citrus trees. Crop Protection 17(4): 289-292.        [ Links ]

DERKSEN, R.C.; BODE, L.E. 1986. Droplet size comparisons from rotary atomizers. Transactions of the American Society of Agricultural and Biological Engineers 29(5): 1204-1207.        [ Links ]

DERKSEN, R.C.; GRAY, L.E. 1995. Deposition and air speed patterns of air-carrier apple orchard sprayers. Transactions of the American Society of Agricultural and Biological Engineers 38(1): 5-11.        [ Links ]

FAROOQ, M.; SALYANI, M. 2002. Spray penetration into the citrus tree canopy from two air-carrier sprayers. Transactions of the American Society of Agricultural and Biological Engineers 45(5): 1287- 1293.        [ Links ]

HALL, F.R.; REICHARD, D.L.; KRUEGER, H.R. 1981. Effects of spray volume and nozzle pressure on orchard spray deposits. Journal of Economic Entomology 74(4): 461-465.        [ Links ]

HOFFMAN, W.C.; SALYANI, M. 1996. Spray deposition on citrus canopies under different meteorological conditions. Transactions of the American Society of Agricultural and Biological Engineers 39(1): 17-22.        [ Links ]

JUSTE, F.; SANCHEZ, S.; IBAÑEZ, R.; VAL, L.; GARCIA, C. 1990. Measurement of spray deposition and efficiency of pesticide application in citrus orchard. Journal of Agricultural Engineering Research 46: 187-196.        [ Links ]

LEFEBVRE, A.H. 1989. Atomization and sprays. Ed. Taylor & Francis, USA.        [ Links ]

MATTH EWS, G.A. 2000. Pesticide application methods. 3rd. ed. Blackwell science, London.        [ Links ]

MATUO , T. 1988. Desenvolvimento de um pulverizador intermitente operado fotoeletricamente para tratamento de pomares de citros. Tesis de libre-docencia, Facultad de Ciencias Agrarias y Veterinaria, Universidade Estadual Paulista, "Julio de Mesquita filho", Jaboticabal, Brasil. 167 p.        [ Links ]

MONTEIRO, M.V. 2006. Compêndio de aviação agrícola. Sorocaba, Brasil. 298 p.        [ Links ]

MORGAN, N.G. 1981. Minimizing pesticide waste in orchard spraying. Outlook on Agriculture 10: 342-344.        [ Links ]

PALLADINI, L.A. 1990. Efeito de condições operacionais de um turboatomizador na cobertura de folhas de citros. Tesis de maestría, Facultad de Ciencias Agrarias y Veterinaria, Universidade Estadual Paulista, "Julio de Mesquita Filho", Jaboticabal, Brasil. 93 p.        [ Links ]

RAETANO , C. 1996. Condições operacionais de turboatomizadores na distribuição e deposição da pulverização em citros. Tesis de doctorado, Escuela Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Piracicaba, Brasil. 93 p.        [ Links ]

RAMOS, H.; PIO, L.C. 2003. Tecnologia de aplicação de produtos fitossanitários. In: Zambolim L, Zuppi M, Da Conceição, Santiago T (Ed.). O que engenheiros agrônomos devem saber para orientar o uso de produtos fitossanitários, Viçosa, Brasil. 376 p.        [ Links ]

RAMOS, H.; KIYOSHI, Y.; CORRÊA, M.; BASSANEZI, R.B.; GARCIA, L.C. 2004. Efeito do volume de calda aplicado com turbopulverizador sobre a deposição e cobertura em folhas, ramos, e frutos de citros. In: Anales del simposio internacional de tecnología de aplicación de agrotóxicos, Botucatu-Brasil: Fepaf, p. 144-147.        [ Links ]

RIGOLIN, A.T.; TERSI, F.E.A. 1999. Redução nos custos de aplicação dos defensivos em citros. Laranja 20(1): 93-106.        [ Links ]

SAS institute. 1999. Sas/stat user´s guide 8.0. Cary: SAS Institute inc., 1 Cd-Rom.        [ Links ]

SALYANI, M. 1988. Droplet size effect on spray deposition efficiency of citrus leaves. Transactions of the American Society of Agricultural and Biological Engineers 31(6): 1680-1684.        [ Links ]

TRAVIS, J.W.; SKROCH, W.A.; SUTTON, T.B. 1987. Effect of canopy density on pesticide deposition and distribution in apple trees. Plant Diseases 71(7): 613-615.        [ Links ]

VANEE, G.R.; LEDEBUHR, R.L.; POTTER, H.S. 1985. Air curtain sprayer increases spraying efficiency. Journal of Agriculture Engineering Research 66(7): 15-17.        [ Links ]

VÁSQUEZ-CASTRO, J.A. 2006. Resíduos de fenitrotion e esfenvalerato em grãos de milho e trigo, em alguns de seus produtos processados e sua ação residual sobre Sitophilus oryzae (L., 1763), Sitophilus zeamais Motsch., 1855 (Coleoptera: Curculionidae) e Rhyzopertha dominica (Fabr., 1792) (Coleoptera: Bostrichidae). Tesis de doctorado, Escuela Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Piracicaba, Brasil. 213 p.        [ Links ]

VÁSQUEZ-CASTRO, J.A. 2008. Tecnología de aplicación de pesticidas. Revista Agronomía 139: 47-52.        [ Links ]

VÁSQUEZ-CASTRO, J.A.; BAPTISTA, G.C.; GADANHA Jr., C.D.; TREVIZAN, L.R.P. 2007. Influence of emulsifiable concentrate formulation on the physical properties of the fluid, spray characteristics, and insecticide deposits on stored grains. Journal of Agricultural and Food Chemistry 55: 3529-3534.        [ Links ]

VOLL, C.E.; VÁSQUEZ-CASTRO, J.A.; GADANHA Jr., C.D.; LINO, P.L. 2006. Distribución volumétrica de la boquilla de pulverización XR 11003. Revista Peruana de Entomología 45:95-99.         [ Links ] ]]> National Food Administration 1998 9-41 1978 2 659-666 Proceedings of the conference by agricultural research institute 1985 25-38 1983 26 2 2 338-342 1998 17 4 4 289-292 1986 29 5 5 1204-1207 1995 38 1 1 5-11 2002 45 5 5 1287- 1293 1981 74 4 4 461-465 1996 39 1 1 17-22 1990 46 187-196 1989 2000 3 1988 2006 298 1981 10 342-344 1990 93 1996 2003 376 2004 144-147 1999 20 1 1 93-106 SAS institute 1999 1988 31 6 6 1680-1684 1987 71 7 7 613-615 1985 66 7 7 15-17 2006 213 2008 139 47-52 2007 55 3529-3534 2006 45 95-99