Re: [R-es] Cómo buscar patrones de comportamiento en datos temporales?

[Juan Abasolo]

Muchas gracias, Rubén y Javier;
Con respecto a lo que decís, Javier compatriota, desde mi desconocimiento,
tu ve esa intuición e intenté encontrar algo por el lado de ecología. Pero
no fui capaz. Es un oceano demasiado grande distinguir en la investigación
de la ecología desde los pocos conocimientos de pedagogía y lingüística.
También se me ocurrió que en geología podría haber algo, pero no llegó a
buen puerto mi búsqueda. No creo que yo vaya saber quebrantar la ley en mis
búsquedas... y no te puedo preguntar cómo hacerlo, por pedir que los hagás
vos! Todo un dilema, Javier. Si ves, dentro de la legalidad, alguna opción
obvia, comentámela.

Como decís, Ruben, presuponemos que los datos tienen que tener dependencia
temporal. Hace un tiempo un compañero entusiasta de Rattle me enseño a
instalarlo y usarlo, así que voy a desempolvar conocimientos, me parece. De
momento, lo que más me desconcierta es cómo incluir el tiempo junto con
variables relacionadas. Si te resulta evidente, te agradeceré todavía más
alguna directriz. Mientras, voy a intentar conseguir el libro que comentás,
a ver qué llego a encontrar.

Que pasen y ayuden a pasar una Feliz Nochebuena!

Saludos.

Hau idatzi du Rubén Fernández Casal (rubenfca...@gmail.com) erabiltzaileak
(2018 abe. 22, lr. (12:45)):

> Hola a todos,
>
> Por lo que comentáis estáis en una etapa muy temprana y podéis tener todo
> tipo de datos. Desde mi punto de vista tener modelos adecuados para este
> tipo de casos puede ser complicado ya que en principio habría una
> dependencia temporal. Mi recomendación sería enfocarlo primero desde un
> punto de vista descriptivo (realmente siempre recomiendo esto...). Además
> de los métodos tradicionales, incluyendo análisis de correlaciones,
> componentes principales y análisis cluster,  los árboles de decisión pueden
> ser de utilidad para empezar. Yo suelo recomendar el entorno gráfico rattle
> (paquete de R con un libro asociado: Data Mining with Rattle and R) y el
> plugin de Rcommander FactoMineR para los que no tienen mucha experiencia,
> pero tienes muchas más opciones...
>
> Un saludo, Rubén.
>
> El vie., 21 dic. 2018 15:46, Javier Marcuzzi <
> javier.ruben.marcu...@gmail.com> escribió:
>
>> Estimado
>> En la parte de ecología hay librerías que le pueden ser útiles, el
>> problema
>> es en conocer primero la parte biológica, yo leí tesis doctoral donde
>> solicité las matrículas profesionales, porque informáticos obiaron puntos
>> que en mi país entran en delitos a la salud pública, el comportamiento
>> tiene causas que se estudian desde medicina, cuidado en ese punto.
>>
>> El vie., 21 de dic. de 2018 11:15 AM, Juan Abasolo 
>> escribió:
>>
>> > Buenas,
>> >
>> > Estoy intentando hacer un trabajo con una compañera, a la que le pedí
>> que
>> > me lo explicara en unas lineas. Esas lineas son las que siguen:
>> >
>> > Quisieramos analizar algunos comportamientos humanos a lo largo del
>> tiempo;
>> > esto es, encontrar patrones de en algunos comportamientos del habla y la
>> > comunicación:  vocalizaciones, gestos (divididos en categorías), mirada
>> (en
>> > categorías).
>> >
>> > Conocemos un programa para windows que se llama Theme, que permite
>> > encontrar patrones y determinar las posibilidades de que se dé ese
>> > comportamiento.
>> >
>> > Quiero ayudar a mi compañera a hacer esa búsqueda en R; conocen Uds.
>> algún
>> > paquete útil para encontrar patrones enter variables, que no sabemos
>> cuáles
>> > estan relacionadas y cuales no.
>> >
>> > Los datos los estamos etiquetando en ELAN primero y después en LAN;
>> esto es
>> > tenemos datos temporales sobre los que se identifican segmentos en los
>> que
>> > se dan comportamientos según categorías (tres variables con ters
>> > categorías, por simplificar).
>> >
>> > No aporto un ejemplo de datos, porque todavía estamos en el proceso de
>> > obtención y querríamos terminar de hacerlo en función de las técnicas de
>> > análisis.
>> >
>> > Muchas gracias, desde ya. Y que tengan buenas fiestas.
>> >
>> >
>> >
>> >
>> >
>> > --
>> > Juan Abasolo
>> >
>> > Hizkuntzaren eta Literaturaren Didaktika Saila | EUDIA ikerketa taldea
>> > Bilboko Hezkuntza Fakultatea
>> > Euskal Herriko Unibertsitatea
>> > UPV/EHU
>> >
>> > Sarriena auzoa z/g 48940 - Leioa (Bizkaia)
>> >
>> > T: (+34) 94 601 7567
>> > Telegram: @JuanAbasolo
>> > Skype: abasolo72
>> >
>> > Tutoretza ordutegia 
>> >
>> > [[alternative HTML version deleted]]
>> >
>> > ___
>> > R-help-es mailing list
>> > R-help-es@r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>> [[alternative HTML version deleted]]
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>

-- 
Juan Abasolo

Hizkuntzaren eta Literaturaren Didaktika Saila | EUDIA ikerketa taldea
Bilboko Hezkuntza Fakultatea
Euskal Herriko 

Re: [R] Issues with R3.5.2

[Jeff Newmiller]

You could delete your 3.5 personal package library (using the File Explorer 
with Run as Admin if necessary) and re-install your packages without running as 
Admin. If that does not work try uninstalling R and re-installing 3.5.1.

On December 22, 2018 8:16:11 PM PST, Janh Anni  wrote:
>This issue only came up after I installed R3.5.2. Never had any
>problems
>with previous installations.  So it is likely a bug in the current
>version.  Any suggestions what to do now?
>
>On Sat, Dec 22, 2018 at 11:06 PM Jeff Newmiller
>
>wrote:
>
>> That normally only occurs if you have at some time used elevated
>> permissions, beyond which point you fall into a downward spiral of
>more
>> permissions trouble. You are apparently already in trouble, whether
>it was
>> of your own making or due to a bug in the installer.
>>
>> Also, never update the system R package library... always use a
>personal
>> library.
>>
>> On December 22, 2018 6:01:44 PM PST, Janh Anni 
>wrote:
>> >Hi Jeff,
>> >
>> >No, during the installation, there was not an option to Run as
>> >Administration.  But *after *installation, I found that if I
>selected
>> >Run
>> >as Administrator, then I could install packages using
>install.packages
>> >as
>> >usual without problems.
>> >
>> >Thanks
>> >Janh
>> >
>> >On Sat, Dec 22, 2018 at 8:26 PM Jeff Newmiller
>> >
>> >wrote:
>> >
>> >> Did you by any chance use Run As Administrator to install R? If so
>> >then
>> >> you need to uninstall it and delete all files created by it (e.g.
>> >> Documents/R/win-lib/3.5/) and re-install using UAC as prompted.
>> >>
>> >> On December 22, 2018 5:10:27 PM PST, Janh Anni
>
>> >wrote:
>> >> >Dear R Experts,
>> >> >
>> >> >I use Windows 10 and just installed the new R version, R3.5.2 but
>> >when
>> >> >I
>> >> >tried to load a data file using read.table, I got an error
>message
>> >like
>> >> >this:
>> >> >
>> >> >*Error in file(file, "rt") : cannot open the connection*
>> >> >*In addition: Warning message:*
>> >> >*In file(file, "rt") :*
>> >> >*  cannot open file 'StreamPCB.dat': No such file or directory*
>> >> >
>> >> >Also, I couldn't install packages using install.packages as
>usual,
>> >> >unless I
>> >> >run R as Administrator
>> >> >
>> >> >I wonder if anyone else had the same issues and any suggestions
>how
>> >to
>> >> >fix?
>> >> >
>> >> >Thanks a lot
>> >> >Janh
>> >> >
>> >> >   [[alternative HTML version deleted]]
>> >> >
>> >> >__
>> >> >R-help@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.
>> >>
>> >> --
>> >> Sent from my phone. Please excuse my brevity.
>> >>
>>
>> --
>> Sent from my phone. Please excuse my brevity.
>>

-- 
Sent from my phone. Please excuse my brevity.

__
R-help@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.

Re: [R] Random seed

[Jeff Newmiller]

It doesn't matter. The whole point is to make the pseudo-random sequence 
repeatable... unless you have a specific reason to avoid repeatability.

On December 22, 2018 5:33:39 PM PST, Steven Yen  wrote:
>I have known from the old days to set a random seed of a LARGE ODD 
>NUMBER. Now I read instructions of set.seed and it requires ANY
>INTEGER. 
>Any idea? Or, does it matter. Thanks.

-- 
Sent from my phone. Please excuse my brevity.

__
R-help@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.

Re: [R] Issues with R3.5.2

[Janh Anni]

This issue only came up after I installed R3.5.2. Never had any problems
with previous installations.  So it is likely a bug in the current
version.  Any suggestions what to do now?

On Sat, Dec 22, 2018 at 11:06 PM Jeff Newmiller 
wrote:

> That normally only occurs if you have at some time used elevated
> permissions, beyond which point you fall into a downward spiral of more
> permissions trouble. You are apparently already in trouble, whether it was
> of your own making or due to a bug in the installer.
>
> Also, never update the system R package library... always use a personal
> library.
>
> On December 22, 2018 6:01:44 PM PST, Janh Anni  wrote:
> >Hi Jeff,
> >
> >No, during the installation, there was not an option to Run as
> >Administration.  But *after *installation, I found that if I selected
> >Run
> >as Administrator, then I could install packages using install.packages
> >as
> >usual without problems.
> >
> >Thanks
> >Janh
> >
> >On Sat, Dec 22, 2018 at 8:26 PM Jeff Newmiller
> >
> >wrote:
> >
> >> Did you by any chance use Run As Administrator to install R? If so
> >then
> >> you need to uninstall it and delete all files created by it (e.g.
> >> Documents/R/win-lib/3.5/) and re-install using UAC as prompted.
> >>
> >> On December 22, 2018 5:10:27 PM PST, Janh Anni 
> >wrote:
> >> >Dear R Experts,
> >> >
> >> >I use Windows 10 and just installed the new R version, R3.5.2 but
> >when
> >> >I
> >> >tried to load a data file using read.table, I got an error message
> >like
> >> >this:
> >> >
> >> >*Error in file(file, "rt") : cannot open the connection*
> >> >*In addition: Warning message:*
> >> >*In file(file, "rt") :*
> >> >*  cannot open file 'StreamPCB.dat': No such file or directory*
> >> >
> >> >Also, I couldn't install packages using install.packages as usual,
> >> >unless I
> >> >run R as Administrator
> >> >
> >> >I wonder if anyone else had the same issues and any suggestions how
> >to
> >> >fix?
> >> >
> >> >Thanks a lot
> >> >Janh
> >> >
> >> >   [[alternative HTML version deleted]]
> >> >
> >> >__
> >> >R-help@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.
> >>
> >> --
> >> Sent from my phone. Please excuse my brevity.
> >>
>
> --
> Sent from my phone. Please excuse my brevity.
>

[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Issues with R3.5.2

[Jeff Newmiller]

That normally only occurs if you have at some time used elevated permissions, 
beyond which point you fall into a downward spiral of more permissions trouble. 
You are apparently already in trouble, whether it was of your own making or due 
to a bug in the installer.

Also, never update the system R package library... always use a personal 
library.

On December 22, 2018 6:01:44 PM PST, Janh Anni  wrote:
>Hi Jeff,
>
>No, during the installation, there was not an option to Run as
>Administration.  But *after *installation, I found that if I selected
>Run
>as Administrator, then I could install packages using install.packages
>as
>usual without problems.
>
>Thanks
>Janh
>
>On Sat, Dec 22, 2018 at 8:26 PM Jeff Newmiller
>
>wrote:
>
>> Did you by any chance use Run As Administrator to install R? If so
>then
>> you need to uninstall it and delete all files created by it (e.g.
>> Documents/R/win-lib/3.5/) and re-install using UAC as prompted.
>>
>> On December 22, 2018 5:10:27 PM PST, Janh Anni 
>wrote:
>> >Dear R Experts,
>> >
>> >I use Windows 10 and just installed the new R version, R3.5.2 but
>when
>> >I
>> >tried to load a data file using read.table, I got an error message
>like
>> >this:
>> >
>> >*Error in file(file, "rt") : cannot open the connection*
>> >*In addition: Warning message:*
>> >*In file(file, "rt") :*
>> >*  cannot open file 'StreamPCB.dat': No such file or directory*
>> >
>> >Also, I couldn't install packages using install.packages as usual,
>> >unless I
>> >run R as Administrator
>> >
>> >I wonder if anyone else had the same issues and any suggestions how
>to
>> >fix?
>> >
>> >Thanks a lot
>> >Janh
>> >
>> >   [[alternative HTML version deleted]]
>> >
>> >__
>> >R-help@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.
>>
>> --
>> Sent from my phone. Please excuse my brevity.
>>

-- 
Sent from my phone. Please excuse my brevity.

__
R-help@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.

Re: [R] Issues with R3.5.2

[Janh Anni]

Hi Jeff,

No, during the installation, there was not an option to Run as
Administration.  But *after *installation, I found that if I selected Run
as Administrator, then I could install packages using install.packages as
usual without problems.

Thanks
Janh

On Sat, Dec 22, 2018 at 8:26 PM Jeff Newmiller 
wrote:

> Did you by any chance use Run As Administrator to install R? If so then
> you need to uninstall it and delete all files created by it (e.g.
> Documents/R/win-lib/3.5/) and re-install using UAC as prompted.
>
> On December 22, 2018 5:10:27 PM PST, Janh Anni  wrote:
> >Dear R Experts,
> >
> >I use Windows 10 and just installed the new R version, R3.5.2 but when
> >I
> >tried to load a data file using read.table, I got an error message like
> >this:
> >
> >*Error in file(file, "rt") : cannot open the connection*
> >*In addition: Warning message:*
> >*In file(file, "rt") :*
> >*  cannot open file 'StreamPCB.dat': No such file or directory*
> >
> >Also, I couldn't install packages using install.packages as usual,
> >unless I
> >run R as Administrator
> >
> >I wonder if anyone else had the same issues and any suggestions how to
> >fix?
> >
> >Thanks a lot
> >Janh
> >
> >   [[alternative HTML version deleted]]
> >
> >__
> >R-help@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.
>
> --
> Sent from my phone. Please excuse my brevity.
>

[[alternative HTML version deleted]]

__
R-help@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.

[R] Random seed

[Steven Yen]

I have known from the old days to set a random seed of a LARGE ODD 
NUMBER. Now I read instructions of set.seed and it requires ANY INTEGER. 
Any idea? Or, does it matter. Thanks.

-- 
st...@ntu.edu.tw (S.T. Yen)


[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Issues with R3.5.2

[Jeff Newmiller]

Did you by any chance use Run As Administrator to install R? If so then you 
need to uninstall it and delete all files created by it (e.g. 
Documents/R/win-lib/3.5/) and re-install using UAC as prompted.

On December 22, 2018 5:10:27 PM PST, Janh Anni  wrote:
>Dear R Experts,
>
>I use Windows 10 and just installed the new R version, R3.5.2 but when
>I
>tried to load a data file using read.table, I got an error message like
>this:
>
>*Error in file(file, "rt") : cannot open the connection*
>*In addition: Warning message:*
>*In file(file, "rt") :*
>*  cannot open file 'StreamPCB.dat': No such file or directory*
>
>Also, I couldn't install packages using install.packages as usual,
>unless I
>run R as Administrator
>
>I wonder if anyone else had the same issues and any suggestions how to
>fix?
>
>Thanks a lot
>Janh
>
>   [[alternative HTML version deleted]]
>
>__
>R-help@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.

-- 
Sent from my phone. Please excuse my brevity.

__
R-help@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.

[R] Issues with R3.5.2

[Janh Anni]

Dear R Experts,

I use Windows 10 and just installed the new R version, R3.5.2 but when I
tried to load a data file using read.table, I got an error message like
this:

*Error in file(file, "rt") : cannot open the connection*
*In addition: Warning message:*
*In file(file, "rt") :*
*  cannot open file 'StreamPCB.dat': No such file or directory*

Also, I couldn't install packages using install.packages as usual, unless I
run R as Administrator

I wonder if anyone else had the same issues and any suggestions how to fix?

Thanks a lot
Janh

[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Glitch in Kruskal-Wallis test?

[Bert Gunter]

... Moreover, you should not analyze proportions in this way, which treats
.5 = 2/4 or .5 = 2000/4000 identically. As David said, you need to work
with a statistician.

Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Sat, Dec 22, 2018 at 7:32 AM David L Carlson  wrote:

> You may need to spend some more time with the statistician who needs to
> see your data. It is not clear if you have a two sample test or a paired
> sample test. Kruskall-Wallis expects data for each observation, not grouped
> data. Without the observations, the test cannot compute the sample size and
> the degrees of freedom. You have run kruskal.test separately on each
> sample. The kruskal.test is designed for comparing two or more samples.
>
> 
> David L Carlson
> Department of Anthropology
> Texas A University
> College Station, TX 77843-4352
>
> -Original Message-
> From: R-help  On Behalf Of Jenny Liu
> Sent: Saturday, December 22, 2018 7:10 AM
> To: Michael Dewey 
> Cc: r-help@r-project.org
> Subject: Re: [R] Glitch in Kruskal-Wallis test?
>
> Hi Michael,
>
> Thank you for your reply! I'm testing the difference in proportions. Temp
> is temperature, and Prop is the proportion of insect pupae that survived at
> that temperature. I was told by a statistician that the K-W was appropriate
> for testing proportions, but perhaps you know of an alternative? I have
> already tested for heteroscedasticity using the Breusch-Pagan test.
>
> Thanks again,
> Jenny
>
>
>
> On Dec 22, 2018 7:38 AM, "Michael Dewey"  wrote:
>
> Dear Jenny
>
> What exactly do you think you are testing here? You are telling K-W you
> have seven groups each with a single value which is not the usual
> situation for K-W.
>
> Michael
>
>
> On 22/12/2018 04:58, Jenny Liu wrote:
> > Hi everyone,
> > I have been running a K-W test with the attached data, PupMort1. My code:
> > kruskal.test(Prop~Temp,data=PupMort1)
> > However, I found that I get the exact same result when I change the
> x-values, as
> > in the attached data PupMort2.
> > Test run with PupMort1Kruskal-Wallis rank sum testdata:  Prop by Temp
>
> > Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> > Test run with PupMort2Kruskal-Wallis rank sum testdata:  Prop by Temp
>
> > Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> > Does anybody know why this is happening?
> > Thank you!
> > Jenny
> >
> >
> > __
> > R-help@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.
> >
>
> --
> Michael
> http://www.dewey.myzen.co.uk/home.html
>
> [[alternative HTML version deleted]]
>
> __
> R-help@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.
>
> __
> R-help@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.
>

[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Glitch in Kruskal-Wallis test?

[David L Carlson]

You may need to spend some more time with the statistician who needs to see 
your data. It is not clear if you have a two sample test or a paired sample 
test. Kruskall-Wallis expects data for each observation, not grouped data. 
Without the observations, the test cannot compute the sample size and the 
degrees of freedom. You have run kruskal.test separately on each sample. The 
kruskal.test is designed for comparing two or more samples. 


David L Carlson
Department of Anthropology
Texas A University
College Station, TX 77843-4352

-Original Message-
From: R-help  On Behalf Of Jenny Liu
Sent: Saturday, December 22, 2018 7:10 AM
To: Michael Dewey 
Cc: r-help@r-project.org
Subject: Re: [R] Glitch in Kruskal-Wallis test?

Hi Michael,

Thank you for your reply! I'm testing the difference in proportions. Temp
is temperature, and Prop is the proportion of insect pupae that survived at
that temperature. I was told by a statistician that the K-W was appropriate
for testing proportions, but perhaps you know of an alternative? I have
already tested for heteroscedasticity using the Breusch-Pagan test.

Thanks again,
Jenny



On Dec 22, 2018 7:38 AM, "Michael Dewey"  wrote:

Dear Jenny

What exactly do you think you are testing here? You are telling K-W you
have seven groups each with a single value which is not the usual
situation for K-W.

Michael


On 22/12/2018 04:58, Jenny Liu wrote:
> Hi everyone,
> I have been running a K-W test with the attached data, PupMort1. My code:
> kruskal.test(Prop~Temp,data=PupMort1)
> However, I found that I get the exact same result when I change the
x-values, as
> in the attached data PupMort2.
> Test run with PupMort1Kruskal-Wallis rank sum testdata:  Prop by Temp

> Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> Test run with PupMort2Kruskal-Wallis rank sum testdata:  Prop by Temp

> Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> Does anybody know why this is happening?
> Thank you!
> Jenny
>
>
> __
> R-help@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.
>

-- 
Michael
http://www.dewey.myzen.co.uk/home.html

[[alternative HTML version deleted]]

__
R-help@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.

__
R-help@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.

Re: [R] Printing with cat in a procedure

[Steven Yen]

Thank you all - print works wonders.

On 12/22/2018 10:36 PM, Eric Berger wrote:
> Hi Steven,
> Here's one way, using print
>
> try5<-function(A,B){
>   C<-A+B
>   #cat("\nA =",A,"\nC = ",C)
>   cat("\nA = ")
>   print(A)
>   cat("\nC = ")
>   print(C)
>   structure(list(A=A,B=B,C=C))
> }
>
> HTH,
> Eric
>
>
> On Sat, Dec 22, 2018 at 4:32 PM Steven Yen  > wrote:
>
> How do I print a matrix running a procedure? In the code below, I
> print
> with the cat command and get a vector (from A and C).
>
> A<-matrix(rpois(16,lambda=5),nrow=4,byrow=T)
> B<-diag(4)
>
> try5<-function(A,B){
>   C<-A+B
>   cat("\nA =",A,"\nC = ",C)
> structure(list(A=A,B=B,C=C))
> }
>
> v<-try5(A,B)
> v$C
>
> -- 
> st...@ntu.edu.tw  (S.T. Yen)
>
>
> [[alternative HTML version deleted]]
>
> __
> R-help@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.
>

-- 
st...@ntu.edu.tw (S.T. Yen)


[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Printing with cat in a procedure

[Jeff Newmiller]

Try using print instead of cat [1], and please read about what the arguments 
are in the help file [2][3] for any function you are using before posting a 
question.

[1] 
https://stackoverflow.com/questions/31843662/what-is-the-difference-between-cat-and-print
[2] ?cat
[3] ?print

On December 22, 2018 6:31:52 AM PST, Steven Yen  wrote:
>How do I print a matrix running a procedure? In the code below, I print
>
>with the cat command and get a vector (from A and C).
>
>A<-matrix(rpois(16,lambda=5),nrow=4,byrow=T)
>B<-diag(4)
>
>try5<-function(A,B){
>  C<-A+B
>  cat("\nA =",A,"\nC = ",C)
>structure(list(A=A,B=B,C=C))
>}
>
>v<-try5(A,B)
>v$C

-- 
Sent from my phone. Please excuse my brevity.

__
R-help@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.

Re: [R] Printing with cat in a procedure

[Rui Barradas]

Hello,

Use print(A) and print(C). cat is meant for simpler objects.

Hope this helps,

Rui Barradas

Às 14:31 de 22/12/2018, Steven Yen escreveu:

How do I print a matrix running a procedure? In the code below, I print
with the cat command and get a vector (from A and C).

A<-matrix(rpois(16,lambda=5),nrow=4,byrow=T)
B<-diag(4)

try5<-function(A,B){
   C<-A+B
   cat("\nA =",A,"\nC = ",C)
structure(list(A=A,B=B,C=C))
}

v<-try5(A,B)
v$C



__
R-help@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.

Re: [R] Printing with cat in a procedure

[Eric Berger]

Hi Steven,
Here's one way, using print

try5<-function(A,B){
  C<-A+B
  #cat("\nA =",A,"\nC = ",C)
  cat("\nA = ")
  print(A)
  cat("\nC = ")
  print(C)
  structure(list(A=A,B=B,C=C))
}

HTH,
Eric


On Sat, Dec 22, 2018 at 4:32 PM Steven Yen  wrote:

> How do I print a matrix running a procedure? In the code below, I print
> with the cat command and get a vector (from A and C).
>
> A<-matrix(rpois(16,lambda=5),nrow=4,byrow=T)
> B<-diag(4)
>
> try5<-function(A,B){
>   C<-A+B
>   cat("\nA =",A,"\nC = ",C)
> structure(list(A=A,B=B,C=C))
> }
>
> v<-try5(A,B)
> v$C
>
> --
> st...@ntu.edu.tw (S.T. Yen)
>
>
> [[alternative HTML version deleted]]
>
> __
> R-help@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.
>

[[alternative HTML version deleted]]

__
R-help@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.

[R] Printing with cat in a procedure

[Steven Yen]

How do I print a matrix running a procedure? In the code below, I print 
with the cat command and get a vector (from A and C).

A<-matrix(rpois(16,lambda=5),nrow=4,byrow=T)
B<-diag(4)

try5<-function(A,B){
  C<-A+B
  cat("\nA =",A,"\nC = ",C)
structure(list(A=A,B=B,C=C))
}

v<-try5(A,B)
v$C

-- 
st...@ntu.edu.tw (S.T. Yen)


[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Glitch in Kruskal-Wallis test?

[Jenny Liu]

Hi Michael,

Thank you for your reply! I'm testing the difference in proportions. Temp
is temperature, and Prop is the proportion of insect pupae that survived at
that temperature. I was told by a statistician that the K-W was appropriate
for testing proportions, but perhaps you know of an alternative? I have
already tested for heteroscedasticity using the Breusch-Pagan test.

Thanks again,
Jenny



On Dec 22, 2018 7:38 AM, "Michael Dewey"  wrote:

Dear Jenny

What exactly do you think you are testing here? You are telling K-W you
have seven groups each with a single value which is not the usual
situation for K-W.

Michael


On 22/12/2018 04:58, Jenny Liu wrote:
> Hi everyone,
> I have been running a K-W test with the attached data, PupMort1. My code:
> kruskal.test(Prop~Temp,data=PupMort1)
> However, I found that I get the exact same result when I change the
x-values, as
> in the attached data PupMort2.
> Test run with PupMort1Kruskal-Wallis rank sum testdata:  Prop by Temp

> Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> Test run with PupMort2Kruskal-Wallis rank sum testdata:  Prop by Temp

> Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
> Does anybody know why this is happening?
> Thank you!
> Jenny
>
>
> __
> R-help@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.
>

-- 
Michael
http://www.dewey.myzen.co.uk/home.html

[[alternative HTML version deleted]]

__
R-help@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.

Re: [R] Glitch in Kruskal-Wallis test?

[Michael Dewey]

Dear Jenny

What exactly do you think you are testing here? You are telling K-W you 
have seven groups each with a single value which is not the usual 
situation for K-W.


Michael

On 22/12/2018 04:58, Jenny Liu wrote:

Hi everyone,
I have been running a K-W test with the attached data, PupMort1. My code:
kruskal.test(Prop~Temp,data=PupMort1)
However, I found that I get the exact same result when I change the x-values, as
in the attached data PupMort2.
Test run with PupMort1Kruskal-Wallis rank sum testdata:  Prop by Temp
Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
Test run with PupMort2Kruskal-Wallis rank sum testdata:  Prop by Temp
Kruskal-Wallis chi-squared = 6, df = 6, p-value = 0.4232
Does anybody know why this is happening?
Thank you!
Jenny


__
R-help@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.



--
Michael
http://www.dewey.myzen.co.uk/home.html

__
R-help@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.

Re: [R] Problem with Kruskal–Wallis test

[Michael Dewey]

Dear Giuseppe

If I understand you correctly you have a very large sample size so it is 
not surprising that you get very small p-values. Eevn a scientifically 
uninteresting difference can become statistically significant with large 
samples. You probably need to define a metric for meaningful differences 
between groups and calculate a confidence interval for it.


Michael

On 21/12/2018 15:37, Giuseppe Cillis wrote:

Dear all,
I am a beginner with R (and also with the statistics) for which I hope to
be clear.
I should do this non-parametric test on data I extracted from maps.
In practice I have a column that represents the landscape Dynamics of a
certain time period (there are 3 dynamics, each of them marked by the
number 1, 2 or 3) and the other column with the values of a topographic
variable (for example the slope) . In all, there are more than 90,000 pairs
of values.
Going to do the test in R, for all the dynamics and for all the variables,
I get out of the values of chi-square elevated (even in the order of
thousands) and a p-value always <2.2e-16  why? Where can the error be? in
the script or in the test approach?
Thanks in advance

[[alternative HTML version deleted]]

__
R-help@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.



--
Michael
http://www.dewey.myzen.co.uk/home.html

__
R-help@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.

Re: [R-es] Cómo buscar patrones de comportamiento en datos temporales?

[Rubén Fernández Casal]

Hola a todos,

Por lo que comentáis estáis en una etapa muy temprana y podéis tener todo
tipo de datos. Desde mi punto de vista tener modelos adecuados para este
tipo de casos puede ser complicado ya que en principio habría una
dependencia temporal. Mi recomendación sería enfocarlo primero desde un
punto de vista descriptivo (realmente siempre recomiendo esto...). Además
de los métodos tradicionales, incluyendo análisis de correlaciones,
componentes principales y análisis cluster,  los árboles de decisión pueden
ser de utilidad para empezar. Yo suelo recomendar el entorno gráfico rattle
(paquete de R con un libro asociado: Data Mining with Rattle and R) y el
plugin de Rcommander FactoMineR para los que no tienen mucha experiencia,
pero tienes muchas más opciones...

Un saludo, Rubén.

El vie., 21 dic. 2018 15:46, Javier Marcuzzi <
javier.ruben.marcu...@gmail.com> escribió:

> Estimado
> En la parte de ecología hay librerías que le pueden ser útiles, el problema
> es en conocer primero la parte biológica, yo leí tesis doctoral donde
> solicité las matrículas profesionales, porque informáticos obiaron puntos
> que en mi país entran en delitos a la salud pública, el comportamiento
> tiene causas que se estudian desde medicina, cuidado en ese punto.
>
> El vie., 21 de dic. de 2018 11:15 AM, Juan Abasolo 
> escribió:
>
> > Buenas,
> >
> > Estoy intentando hacer un trabajo con una compañera, a la que le pedí que
> > me lo explicara en unas lineas. Esas lineas son las que siguen:
> >
> > Quisieramos analizar algunos comportamientos humanos a lo largo del
> tiempo;
> > esto es, encontrar patrones de en algunos comportamientos del habla y la
> > comunicación:  vocalizaciones, gestos (divididos en categorías), mirada
> (en
> > categorías).
> >
> > Conocemos un programa para windows que se llama Theme, que permite
> > encontrar patrones y determinar las posibilidades de que se dé ese
> > comportamiento.
> >
> > Quiero ayudar a mi compañera a hacer esa búsqueda en R; conocen Uds.
> algún
> > paquete útil para encontrar patrones enter variables, que no sabemos
> cuáles
> > estan relacionadas y cuales no.
> >
> > Los datos los estamos etiquetando en ELAN primero y después en LAN; esto
> es
> > tenemos datos temporales sobre los que se identifican segmentos en los
> que
> > se dan comportamientos según categorías (tres variables con ters
> > categorías, por simplificar).
> >
> > No aporto un ejemplo de datos, porque todavía estamos en el proceso de
> > obtención y querríamos terminar de hacerlo en función de las técnicas de
> > análisis.
> >
> > Muchas gracias, desde ya. Y que tengan buenas fiestas.
> >
> >
> >
> >
> >
> > --
> > Juan Abasolo
> >
> > Hizkuntzaren eta Literaturaren Didaktika Saila | EUDIA ikerketa taldea
> > Bilboko Hezkuntza Fakultatea
> > Euskal Herriko Unibertsitatea
> > UPV/EHU
> >
> > Sarriena auzoa z/g 48940 - Leioa (Bizkaia)
> >
> > T: (+34) 94 601 7567
> > Telegram: @JuanAbasolo
> > Skype: abasolo72
> >
> > Tutoretza ordutegia 
> >
> > [[alternative HTML version deleted]]
> >
> > ___
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
>
> [[alternative HTML version deleted]]
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es