Description Usage Arguments Details Value Examples

estimates linear contrasts of the elements of c, c_s, c_t, or c_st from a `predPref`

object

1 2 |

`x` |
a predPref object as fit by the eponymous function |

`b` |
a vector to linearly transform c_st |

`mu` |
a number to test the linear contrast against in the null |

`alternative` |
string to specify alternative hypothesis |

`conf.level` |
confidence level of the interval |

`sig.level` |
determines null/alternative hypothesis value of c_st from predPref |

The input vector b performs the linear transformation t(b) %*% matrix(c_st), so that c_st becomes a column vector by indexing t first and then s. Hence there is no requirement of a linear contrast, any linear transformation such that t(b) %*% matrix(1, nrow=length(b)) != 0 is allowed.

Of the two estimated hypotheses in the underlying call
to `predPref`

, the linear transformation b is applied to the
hypothesis that is determined by the choice of `sig.level`

.

A list with class '"htest"' containing the following components:

statistic: the value of the t-statistic.

parameter: the degrees of freedom for the t-statistic.

p.value: the p-value for the test.

conf.int: a confidence interval for the mean appropriate to the specified alternative hypothesis.

estimate: the estimated mean or difference in means depending on whether it was a one-sample test or a two-sample test.

null.value: the specified hypothesized value of the mean or mean difference depending on whether it was a one-sample test or a two-sample test.

alternative: a character string describing the alternative hypothesis.

method: a character string indicating what type of t-test was performed.

data.name: a character string giving the names of the data.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ```
# set parameters
Predators <- Traps <- 100
PreySpecies <- 2
Times <- 5
g <- matrix(sqrt(2), nrow=Times, ncol=PreySpecies) # gamma
l <- matrix(seq(0.4,1.8,length.out=5)*sqrt(2), nrow=Times, ncol=PreySpecies) # ct
# fit model and contrast
## Not run:
set.seed(0)
fdata <- simPref(PreySpecies, Times, Predators, Traps, l, g, EM=FALSE) # p-value=0.305
pref <- predPref(fdata$eaten, fdata$caught, hypotheses=c('ct', 'cst'))
testC(pref, b = c(0,1, -1, 0, 0)) # p-value > sig.level => ct is used, not cst
## End(Not run)
``` |

spiders documentation built on May 30, 2017, 3:31 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.