[R-sig-eco] cozigam

2012-06-11 Thread Mahnaz Rabbaniha
hi

i try to find regression between clupeidae,with temperature,salinity and
depth. the response variable is inclued many zero ( 86 from 133 observed)

therefore i used this code :

 res - cozigam(Clupeidae~s(temperature,salinity)+s(depth), constraint =
proportional, family = gaussian)


the result:


iteration = 2norm = 1.001743

iteration = 3norm = 0.3377464

iteration = 4norm = 9.172232e-05



==

estimated alpha = -0.5337883 ( 0.1789113 )

estimated delta = -0.0009891505 ( NaN )

==



Warning message:

In sqrt(V.theta[2, 2]) : NaNs produced


what is exactly meaning?


thanks for your help

-- 
Mahnaz

[[alternative HTML version deleted]]

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


[R-sig-eco] error message with function rda of package vegan

2012-06-11 Thread Duarte Viana
Hello all,

Apologies for cross postings.

I am running the function rda of the vegan library, like I did many
times without having troubles, and I get the following error message:
Error in if (rank) { : argument is not interpretable as logical.

I am just testing a response matrix against a single explanatory
matrix, which is the result of a model matrix with 4 dummy variables
(coding for a factor of 5 levels).

myrda-rda(resp, expl)

When I, for instance, change the response matrix, it goes well.

Does anyone know what might be happening?

Thanks,

Duarte

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


Re: [R-sig-eco] error message with function rda of package vegan

2012-06-11 Thread Sarah Goslee
Hi,

On Mon, Jun 11, 2012 at 12:35 PM, Duarte Viana viana.s...@gmail.com wrote:
 Hello all,

 Apologies for cross postings.

 I am running the function rda of the vegan library, like I did many
 times without having troubles, and I get the following error message:
 Error in if (rank) { : argument is not interpretable as logical.

 I am just testing a response matrix against a single explanatory
 matrix, which is the result of a model matrix with 4 dummy variables
 (coding for a factor of 5 levels).

 myrda-rda(resp, expl)

 When I, for instance, change the response matrix, it goes well.

 Does anyone know what might be happening?

Not without a reproducible example, I'm afraid.

Sarah

-- 
Sarah Goslee
http://www.functionaldiversity.org

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


[R-sig-eco] error message with function rda of package vegan

2012-06-11 Thread Duarte Viana
Hi,

Thank you for your reply. It's not easy to give you a reproducible
example, because this is a very specific case (something in my data ,
I believe). The only way is to send the data and try yourself. If
someone is willing to do that, I would be grateful.

I can send the response matrix and the explanatory matrix, and all one
has to do is load library vegan, load the data, and run the rda
function:

library(vegan)

resp-read.table(resp.txt, head=F, sep=\t)
expl-read.table(expl.txt, head=F, sep=\t)

myrda-rda(resp, expl)

Thanks,

Duarte



On Mon, Jun 11, 2012 at 6:42 PM, Sarah Goslee sarah.gos...@gmail.com wrote:
 Hi,

 On Mon, Jun 11, 2012 at 12:35 PM, Duarte Viana viana.s...@gmail.com wrote:
 Hello all,

 Apologies for cross postings.

 I am running the function rda of the vegan library, like I did many
 times without having troubles, and I get the following error message:
 Error in if (rank) { : argument is not interpretable as logical.

 I am just testing a response matrix against a single explanatory
 matrix, which is the result of a model matrix with 4 dummy variables
 (coding for a factor of 5 levels).

 myrda-rda(resp, expl)

 When I, for instance, change the response matrix, it goes well.

 Does anyone know what might be happening?

 Not without a reproducible example, I'm afraid.

 Sarah

 --
 Sarah Goslee
 http://www.functionaldiversity.org

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


Re: [R-sig-eco] error message with function rda of package vegan

2012-06-11 Thread Simpson, Gavin
Feel free to send your data to me, off list. I'll take a look at what might be 
going on inside vegan.

Gavin

Sent from my HTC

- Reply message -
From: Duarte Viana viana.s...@gmail.com
Date: Mon, Jun 11, 2012 20:47
Subject: [R-sig-eco] error message with function rda of package vegan
To: Sarah Goslee sarah.gos...@gmail.com
Cc: r-sig-ecology@r-project.org r-sig-ecology@r-project.org

Hi,

Thank you for your reply. It's not easy to give you a reproducible
example, because this is a very specific case (something in my data ,
I believe). The only way is to send the data and try yourself. If
someone is willing to do that, I would be grateful.

I can send the response matrix and the explanatory matrix, and all one
has to do is load library vegan, load the data, and run the rda
function:

library(vegan)

resp-read.table(resp.txt, head=F, sep=\t)
expl-read.table(expl.txt, head=F, sep=\t)

myrda-rda(resp, expl)

Thanks,

Duarte



On Mon, Jun 11, 2012 at 6:42 PM, Sarah Goslee sarah.gos...@gmail.com wrote:
 Hi,

 On Mon, Jun 11, 2012 at 12:35 PM, Duarte Viana viana.s...@gmail.com wrote:
 Hello all,

 Apologies for cross postings.

 I am running the function rda of the vegan library, like I did many
 times without having troubles, and I get the following error message:
 Error in if (rank) { : argument is not interpretable as logical.

 I am just testing a response matrix against a single explanatory
 matrix, which is the result of a model matrix with 4 dummy variables
 (coding for a factor of 5 levels).

 myrda-rda(resp, expl)

 When I, for instance, change the response matrix, it goes well.

 Does anyone know what might be happening?

 Not without a reproducible example, I'm afraid.

 Sarah

 --
 Sarah Goslee
 http://www.functionaldiversity.org

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


[[alternative HTML version deleted]]

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