search for: dfr1

...form: resmat %*% cfs = target # 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,dfr...

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 <-

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

Dear all, I have a problem that may be someone of you can help. I am a newbie and do not find how to do it in manuals. I have two arrays, for example: ar1 <- array(data=c(1:16),dim=c(4,4)) ar2 <- array(data=c(1:16),dim=c(4,4)) colnames(ar1)<-c("A","B","D","E") colnames(ar2)<-c("C","A","E","B") > ar1

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