Displaying 2 results from an estimated 2 matches for "arcdi".
Did you mean:
ardi
2011 Feb 10
2
Rioja package, creating transfer function, WA, "Error in FUN"
...EST) (STEST is the
> species abundance file)
> Spe=sqlFetch(STEST, "Sheet8")
> odbcClose(STEST)
> Spe
(The species data contains the abundance of 32 species over 90 sites, set
out like this)
F1 AmpFlav AmpWri ArcCat ArcDis
1 N1_1 22.2929936 0.0000000 0.0000000 0.0000000
2 N1_2 30.9677419 0.0000000 0.0000000 3.2258065
> library(rioja)
> y <-as.matrix(Spe)
> x <-as.matrix(Env)
> WA(y, x, tolDW = FALSE, use.N2=TRUE, check.data=TRUE, lean=FALSE) (the
> command from the WA section...
2011 Feb 22
1
Transfer function observed vs predicted values graph problem
...t;30" "90.04269" "3.745020" "108.57278"
....
90 "GAL_15" "70" "94.07849" "3.777932" " 66.77673"
The species matrix (y) contains the abundance of 32 species over 90 sites,
set out like this
F1 AmpFlav AmpWri ArcCat ArcDis
1 1 22.2929936 0.0000000 0.0000000 0.0000000
2 2 30.9677419 0.0000000 0.0000000 3.2258065
fit <- WA(y, x, tolDW = FALSE, use.N2=TRUE, check.data=TRUE, lean=FALSE)
# plot predicted vs. observed
plot(fit)
plot(fit, resid=TRUE)
# Water-table reconstruction
pred <- predict(fit, y)
#plot the re...