[R] Trouble with SPI package
Jeff Newmiller
jdnewmil at dcn.davis.CA.us
Mon May 25 19:56:04 CEST 2015
Why are you overwriting your input file after you read it in? This seems likely to end up corrupting your input data if you make a mistake.
For instance, when you read a file into a data frame that has pure numeric values in the header line, the default behaviour is to convert those numeric values to valid labels, which means that they must start with a letter (in this case X).
I have never used the spi function, but its help file example for the filename argument suggests that the year numbers should not have letters in them.
You should be able to use a text editor to repair your input file, and remove the line that overwrites it from your code.
---------------------------------------------------------------------------
Jeff Newmiller The ..... ..... Go Live...
DCN:<jdnewmil at dcn.davis.ca.us> Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Embedded Controllers) .OO#. .OO#. rocks...1k
---------------------------------------------------------------------------
Sent from my phone. Please excuse my brevity.
On May 25, 2015 8:56:32 AM PDT, Miller Andres Ruiz Sanchez <mruiz at cenipalma.org> wrote:
>Hello,
>
>I write to ask you about an error that I have when I use the script
>below.
>I'm working with monthly precipitation data for the period between
>1990
>and 1998.
>
>I really thanks your help.
>
>_______________________________________________________________
>
>> require(SPEI)
>> require(spi)
>> require(sm)
>
>> dir()
>[1] "IndexScript.R" "PradoCorre.txt" "spi.txt"
>
>> Prado=read.table("PradoCorre.txt", head=TRUE, dec=".")
>
>> write.table(Prado,file="PradoCorre.txt",quote=FALSE,row.names=TRUE)
>
>> spi(3,"PradoCorre.txt",1990,1998)
>Error en data[i, ] : subíndice fuera de los límites
>
>
>> Prado
> Months X1990 X1991 X1992 X1993 X1994 X1995 X1996 X1997 X1998
>1 Jan 0.0 0.00 0.0 22.3 0.0 0.0 15.4 0.0 0.0
>2 Feb 0.0 11.00 0.0 0.0 0.0 0.0 0.0 2.5 0.0
>3 Mar 0.0 8.70 0.0 13.1 0.3 0.4 34.3 0.0 3.8
>4 Apr 52.0 32.20 96.8 70.0 61.4 251.0 21.0 31.0 18.0
>5 May 130.0 34.11 249.1 348.4 211.0 141.5 144.8 36.3 314.3
>6 Jun 102.0 142.20 188.5 70.6 24.1 116.9 90.6 159.4 126.9
>7 Jul 86.0 98.00 80.6 89.0 39.9 228.0 234.0 26.2 27.8
>8 Aug 135.0 82.78 76.3 173.5 245.9 370.9 44.2 51.8 162.4
>9 Sep 132.0 103.30 216.4 214.3 120.0 177.7 84.3 132.3 403.9
>10 Oct 432.0 245.60 180.6 253.2 101.2 356.1 241.9 92.8 117.8
>11 Nov 42.6 24.03 48.3 64.3 155.1 15.0 120.0 14.2 0.0
>12 Dec 109.3 1.60 0.0 26.0 0.0 9.3 0.0 0.0 0.0
>>
>
>
>________________________________________________________________
>
>
>
>
>Cordialmente
>
>*Miller Ruiz*
>Ingeniero Agrónomo
>Auxiliar de Investigación
>Área de Geomática
>Cenipalma - Zona Norte
>Cel: 3164807973
>Fundación - Magdalena
>
>--
>[image: Fedepalma]
><http://www.google.com/url?q=http%3A%2F%2Fwww.fedepalma.org&sa=D&sntz=1&usg=AFrqEzcV3Ok1p8dfn_mAY_mbeKqx9TcBUA>
>
><http://www.google.com/url?q=http%3A%2F%2Fwww.cenipalma.org%2F&sa=D&sntz=1&usg=AFrqEzfg7qal8zogbo9Aq6drmj-iPne9eQ>
>El presente mensaje de datos como sus archivos adjuntos se consideran
>información confidencial y de valor estratégico para la Federación;
>motivo
>por el cual sólo podrá ser empleada por su exclusivo destinatario,
>según
>las indicaciones impartidas por el remitente y dada la naturaleza de la
>
>misma. En consecuencia, cualquier uso, explotación, reproducción,
>modificación, distribución, puesta a disposición entre otras
>posibilidades
>diferentes a las autorizadas, se entenderán expresamente prohibidas. Si
>
>Usted No es el destinatario, deber eliminar completamente el mensaje
>y, en
>lo posible, notificar al remitente del mismo. Es responsabilidad del
>destinatario de este mensaje comprobar que el mensaje de datos y sus
>adjuntos no representan un riesgo informático para su propio sistema.
>
>La Federación ha tomado las medidas adecuadas para tratar de prevenir
>la
>transmisión de virus y programas malignos, no obstante, no se hace
>responsable por su eventual transmisión por este conducto. La
>Federación
> no acepta responsabilidad alguna por eventuales daños o alteraciones
>derivadas de la recepción o uso del presente mensaje.
>
>--
>
><http://www.google.com/url?q=http%3A%2F%2Fwww.cenipalma.org%2F&sa=D&sntz=1&usg=AFrqEzfg7qal8zogbo9Aq6drmj-iPne9eQ>
>El presente mensaje de datos como sus archivos adjuntos se consideran
>información confidencial y de valor estratégico para la Federación;
>motivo
>por el cual sólo podrá ser empleada por su exclusivo destinatario,
>según
>las indicaciones impartidas por el remitente y dada la naturaleza de la
>
>misma. En consecuencia, cualquier uso, explotación, reproducción,
>modificación, distribución, puesta a disposición entre otras
>posibilidades
>diferentes a las autorizadas, se entenderán expresamente prohibidas. Si
>
>Usted No es el destinatario, deber eliminar completamente el mensaje
>y, en
>lo posible, notificar al remitente del mismo. Es responsabilidad del
>destinatario de este mensaje comprobar que el mensaje de datos y sus
>adjuntos no representan un riesgo informático para su propio sistema.
>
>La Federación ha tomado las medidas adecuadas para tratar de prevenir
>la
>transmisión de virus y programas malignos, no obstante, no se hace
>responsable por su eventual transmisión por este conducto. La
>Federación
> no acepta responsabilidad alguna por eventuales daños o alteraciones
>derivadas de la recepción o uso del presente mensaje.
>
> [[alternative HTML version deleted]]
>
>______________________________________________
>R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>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.
More information about the R-help
mailing list