Displaying 6 results from an estimated 6 matches for "dfr2".
Did you mean:
df2
2001 Jul 18
1
hypothesis testing in models
...get
# resmat: matrix/vector of restriction multipliers (i x k) ARRAY!!
# target: restriction value vector/scalar (k)
htest.lm <- function(object, resmat, target = 0, level = 0.95)
{
tsu <- summary(object, corr=FALSE)
cfs <- as.array(object$coefficients)
dfr1 <- nrow(resmat)
dfr2 <- object$df.residual
S <- aperm(resmat %*% aperm(cfs) - target) %*% solve(resmat %*% tsu$cov.unscaled %*% aperm(resmat)) %*% (resmat %*% aperm(cfs) - target) / (nrow(resmat) * tsu$sigma^2)
cat(" F(",dfr1,",",dfr2,") =",qf(level,dfr1,dfr2), " vs. S =&...
2009 May 07
2
Show name of dataset in graph
Hi,
I?ve written a script to run several multivariate statistical analysis
automatically.
As one result a biplot and screeplot is produced.
Now I?d like to display the name of the inputdatset as part of the title of
these graphics and I do not want to enter it each time I run the script.
How can I extract the name of a dataset?
An (shortened) extraction from the script:
>Test1 <-
2008 Dec 03
1
reduce limit number of arguments in methods:::cbind
...nstruction of the function
which generates nested loops.
A workaround could be to use the internal cbind function on blocks of
non S4 objects. The limitation would then be reduced to the number of
consecutive S4 objects.
##### R code #####
dfr <- data.frame(matrix(0, nrow = 1 , ncol = 1000))
dfr2 <- is.na(dfr)
mlist <- rep(list(matrix(0, 2, 1)), 400)
cb1 <- do.call("cbind", c(mlist, mlist))
methods:::bind_activation(TRUE)
dfr2 <- is.na(dfr) # fails
cb2 <- do.call("cbind", mlist) # ok
cb3 <- do.call("cbind", c(mlist, mlist)) # fails
# This...
2008 Nov 18
2
sequencially merge multiple files in a folder
Dear all,
If the question is too easy, please forgive me since I am only few weeks old in R.
I have worked on this question a few days and still cannot figure it out.
Here I have a folder with more than 50 tab-delimited files. Each file has a few hundreds of thousands rows/subjects, and the number of columns/variables of each file varies.The 1st row consists of all the variable names.
Now I
2009 Oct 21
3
Missing data and LME models and diagnostic plots
Hello
Running R2.9.2 on Windows XP
I am puzzled by the performance of LME in situations where there are missing data. As I understand it, one of the strengths of this sort of model is how well it deals with missing data, yet lme requires nonmissing data.
Thus,
m1.mod1 <- lme(fixed = math_1 ~ I(year-2007.5)*TFC_,
data = long,
random =
2005 Nov 03
3
newbie graphics question: Two density plots in same frame ?
I swear I've scoured the help files and several texts before posting
what feels like a dumb newbie question.
How can I draw two kernel density plots in the same frame ? I have
similar variables in two separate data frames, and I would like to show
their two histograms/densities in a single picture. Same units, scale,
range for both, so I'm simply trying to draw one and then add the