For the record, i added EcoGenetics to the parcours (instead of spdep) for my real data (30 coordinates, distance range 0.5m to 11m, z is species abundance), and forced it to 15 distance bins (for ncf, i had to give a increment leading to 15 bins). The results were roughly comparable, if not for slightly different distance bin means. I find it curious that pgirmess doubled the amount of comparisons (as if a distance matrix wasnt split by diagonale):

um(unlist(df.eco$Z.size))

[1] 435

> sum(unlist(df.ncf$n))

[1] 435

> sum(unlist(df.pgir$n))

[1] 870

Anyway, here is the plot. Do you have an idea (if you even look at this post anymore^^), to force the three packages to use the exact same distances?

]]>How can we interpret the results?

]]>I´ve just found this article and I´m writing this email hoping you could shed some light on an analysis I´m performing.

I am trying to analyze some data about animal behaviour and would need some help or advice regarding which non-parametric test should I use.

The variables I have are:

-Response variable: a continuous one (both positive and negative)

-Explicatory variable: a factor with 6 levels

-Random effect variable: as the same animal performing some behavioural task was measured more than once.

As I have a random effect variable, I chose a GLM model. Then, when checking the normality and homoscedasticity assumptions, Shapiro-Wilks test showed there was no normality and QQplots revealed there weren´t patterns nor outliers in my data. So the question would be: which non-parametric test would be optimal in this case, knowing that I would like to perform certain a posteriori comparisons (and not all-against-all comparisons)?

My database has lots of zeros responses in some conditions, I´ve read that for t-students tests lacking of normality due to lots of zeros it´s OK to turn a blind eye on lack of normality (Srivastava, 1958; Sullivan & D'agostino, 1992) ... is there something similar with GLM?

Thank you so much in advance for any advice you could provide.

Kind regards,

Yair Barnatan

Ph.D. Student - Physiology and Molecular Biology Department

Faculty of Science

University of Buenos Aires

in google's search results, but there is a tool that can help you to

create 100% unique articles, search for: boorfe's tips unlimited content ]]>

Is it possible to extract the original variables from the "new" variables sorted by decreasing eigenvalues (V1 to V4) ?

Any guidance available ?

Thank you! ]]>

The help to the jagam() function states that you can use "s", "te", "ti" or "t2" splines. So that is probably why it does not work with "ns" splines.

One solution would be not to use JAGS, but STAN for that. There is a packages "brms" which allows you to fit GAMs with any kinds of splines in the Bayesian setting. It's great!

Cheers, Petr ]]>

Like:

jagam(y ~ x1 + ns(x2, df=7), family=possion) ]]>

Anyone knows any benchmarks difference? ]]>