Displaying 6 results from an estimated 6 matches for "opcons".
Did you mean:
opcon
2012 Apr 05
4
Appropriate method for sharing data across functions
In trying to streamline various optimization functions, I would like to have a scratch pad
of working data that is shared across a number of functions. These can be called from
different levels within some wrapper functions for maximum likelihood and other such
computations. I'm sure there are other applications that could benefit from this.
Below are two approaches. One uses the <<-
2005 Jun 23
4
contrats hardcoded in aov()?
...<rjvbertin at gmail.com> wrote:
> Hello,
>
> I was just having a look at the aov function source code, and see that when the model used does not have an Error term, Helmert contrasts are imposed:
>
> if (is.null(indError)) {
> ...
> }
> else {
> opcons <- options("contrasts")
> options(contrasts = c("contr.helmert", "contr.poly"))
> on.exit(options(opcons))
> ...
>
>
> My reading of several contributed user guides' sections on ANOVA is that Helmert contrasts are not intuitiv...
2012 Mar 25
2
avoiding for loops
I have data that looks like this:
> df1
group id
1 red A
2 red B
3 red C
4 blue D
5 blue E
6 blue F
I want a list of the groups containing vectors with the ids. I am
avoiding subset(), as it is
only recommended for interactive use. Here's what I have so far:
df1 <- data.frame(group=c("red", "red", "red", "blue",
2004 Jan 30
0
Two apparent bugs in aov(y~ *** -1 + Error(***)), with suggested (PR#6510)
...t;))
c("maov", "aov", oldClass(fit)) else c("aov", oldClass(fit))
fit$call <- Call
return(fit)
} else {
## helmert contrasts can be helpful: do we want to force them?
## this version does for the Error model.
opcons <- options("contrasts")
options(contrasts=c("contr.helmert", "contr.poly"))
on.exit(options(opcons))
allTerms <- Terms
errorterm <- attr(Terms, "variables")[[1 + indError]]
eTerm <- deparse(errorterm[[2]],...
2004 Feb 02
0
Two apparent bugs in aov(y~ *** -1 + Error(***)), with (PR#6520)
...quot;, "aov", oldClass(fit)) else c("aov", oldClass(fit))
> fit$call <- Call
> return(fit)
> } else {
> ## helmert contrasts can be helpful: do we want to force them?
> ## this version does for the Error model.
> opcons <- options("contrasts")
> options(contrasts=c("contr.helmert", "contr.poly"))
> on.exit(options(opcons))
> allTerms <- Terms
> errorterm <- attr(Terms, "variables")[[1 + indError]]
> eTerm <-...
2001 Dec 17
1
environments again
In a previous message I was not clear enough in my querry.
I have the following program:
tst<- function() {
x <- c(32.7,32.3,31.5,32.1,29.7,29.1,35.7,35.9,33.1,
36.0,34.2,31.2,31.8,28.0,29.2,38.2,37.8,31.9,
32.5,31.1,29.7)
g <- rep(1:7,rep(3,7))
s <- rep(1:3,7)
cat(" Only x and g \n")
aov1(x,g)
cat("\n\n Now x, g and s \n")
aov1(x,g,s=s)
}