Dear R users,

I am analyzing data using a COZIGAM  model

The model is:

> res <- cozigam(Clupeidae~s(temperature,salinity), constraint =
"proportional", family = gaussian)


Clupeidae = abundance (from 132 observed 85 items=0 ), LDH = log of a
distance that I am interested in, and
SI3 = is a index of invasion severity, but these details are not important.

The result:

iteration = 2    norm = 0.9998562

iteration = 3    norm = 0.120129

iteration = 4    norm = 0.01653977

iteration = 5    norm = 0.000226035



==========================================

estimated alpha = -0.486186 ( 0.2076793 )

estimated delta = 0.01016082 ( 0.004780538 )


I am trying to plot for the effects of temperature and salinityindividually on
Clupeidae through the code below:


>summary(res)

Family: gaussian

Link function: identity

Formula:

Clupeidae ~ s(temperature, salinity)

Parametric coefficients:

                       Estimate     Std. Error    t value    Pr(>|t|)

(Intercept)    -20.665505  13.904231  -1.486   0.1401

alpha            -0.486186   0.207679    -2.341    0.0210 *

delta1            0.010161   0.004781     2.125    0.0358 *

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Approximate significance of smooth terms:

                                              Edf          Est.rank    F
p-value

s(temperature,salinity)   20.81       29 8.77      <2e-16 ***

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Scale est. = 270.56    n = 132


my problem is the interpret of the results ,is the result acceptable?

For checking the good fitness of result, which code is used? For example
:which code for  plotting of the residuals and the plot of residuals vs.
fitted Values is applied?


Thanks a bunch for any suggestions that anyone can offer.


-- 
Mahnaz

	[[alternative HTML version deleted]]

