Displaying 6 results from an estimated 6 matches for "opcon".
Did you mean:
opcion
2012 Apr 05
4
Appropriate method for sharing data across functions
...mber 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 <<- assignment to a structure I call OPCON. The
other attempts to create an environment with this name, but fails. Though I have looked at
a number of references, I have so far not found an adequate description of how to specify
where the OPCON environment is located. (Both the green and blue books do not cover this
topic, at least not...
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 intuiti...
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)
}