search for: overall

...legant way of creating a summary count (other than extracing it from ftable). I''d like to create a new data.frame containing counts of spell by loc ie have three columns showing spell,loc,count. Below the data.frame... Any help appreciated Thanks Herry spell loc 101 Parts 1 102 Overall 2 105 Parts 1 106 None 1 111 None 1 116 Parts 1 118 None 1 119 Overall 4 123 Overall 1 125 Parts 1 126 Overall 1 127 Parts 1 128 Overall 1 134 Overall 1 138 Overall 2 139 Overall 1 142 Overall 1 191 Parts 1 192 Parts 2 193 Parts 2 204 P...

Dear all, I have a 3x2 plot and in addition to the title of the individual plots I would like to have an overall title for each row. I managed to get an overall title for the whole plot matrix with mtext: par(mfrow=(c(3,2)), mar=c(6.4,4.5,4.2, 1.8), oma=c(0,0,3,0)) for (i in 1:6) barplot(sample(1:10,5), main=paste("Plot #",i)) mtext("Overall Title Row 1", outer=TRUE) but I cannot put...

...sibinomial))[1]) # save warnings wrn1 <- warnings(p) ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1]) ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2]) ## sample counts n <- unwtd.count(~grp, grp1)[1] ## combine into table overall <- data.frame(n, p, ci_l, ci_u) colnames(overall) <- c("counts", "Group1", "LL", "UL") ## add any warnings ind <- length(wrn1) ind if (ind == 0) { msg <- "No warnings" } if (ind > 0) {msg <- names(warnings()) } overall[1,5] <-...

Hello fellow R users, I am having a problem finding the estimates for some overall treatment effects for my mixed models using 'lme' (package nlme). I hope someone can help. Firstly then, the model: The data: Plant biomass (log transformed) Fixed Factors: Treatment(x3 Dry, Wet, Control) Year(x8 2002-2009) Random Factors: 5 plots per treatment, 5 quadrats per plot (N=...

...;- warnings(p) > > ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1]) > ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2]) > > ## sample counts > n <- unwtd.count(~grp, grp1)[1] > > ## combine into table > overall <- data.frame(n, p, ci_l, ci_u) > colnames(overall) <- c("counts", "Group1", "LL", "UL") > > ## add any warnings > ind <- length(wrn1) > ind > > if (ind == 0) { msg <- "No warnings" } > if (ind > 0) {msg <-...

...;- warnings(p) > > ci_l <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1]) > ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2]) > > ## sample counts > n <- unwtd.count(~grp, grp1)[1] > > ## combine into table > overall <- data.frame(n, p, ci_l, ci_u) > colnames(overall) <- c("counts", "Group1", "LL", "UL") > > ## add any warnings > ind <- length(wrn1) > ind > > if (ind == 0) { msg <- "No warnings" } > if (ind > 0) {msg <-...

I'm making some plots on the same page and would like to include an overall title instead of individual main titles as they are similar and their x and y axis labels are sufficient to distinguish them. Is there a way to assign an overall "main" to this page of plots? Mark -- Mark W. Kimpel MD ** Neuroinformatics ** Dept. of Psychiatry Indiana University Schoo...

...(confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[1]) > > ci_u <- (confint(svyciprop(~grp, grp1, family=quasibinomial), 'ci')[2]) > > > > ## sample counts > > n <- unwtd.count(~grp, grp1)[1] > > > > ## combine into table > > overall <- data.frame(n, p, ci_l, ci_u) > > colnames(overall) <- c("counts", "Group1", "LL", "UL") > > > > ## add any warnings > > ind <- length(wrn1) > > ind > > > > if (ind == 0) { msg <- "No warnings"...

Hello, I looking for a solution to display overall rsync progress on an LCD display as a bargraph. I have found 2 parameters: --progress This option tells rsync to print information showing the progress of the transfer. This gives a bored user something to watch. Implies --verbose if it...

I entered the following: formula<-nst~age+soc+inc+reg+imp pnstlm<-lm(formula,nst) summary(pnstlm) imp and soc are ordered categorical variables but the summary does not give an output of the overall p-values, just individual comparisons. I can't find help for this in the manual. Is there a command or option change in the summary to find out the overall p-value? thank you in advance [[alternative HTML version deleted]]

...ion ##LH2 dimLH2 <- dim(y.modwt$LH2) dimLH2x <- dimLH2[1] dimLH2y <- dimLH2[2] varLH2xlist <- c(rep(0, dimLH2x)) varLH2ylist <- c(rep(0, dimLH2y)) ##Loop to get variance from x axis for(i in seq(dimLH2x)){ varLH2xlist[i] <- var(y.modwt$LH2[i,]) } ##Get the variance from the overall x variance varLH2x <- var(varLH2xlist) ##Loop to get variance from y axis for(i in seq(dimLH2y)){ varLH2ylist[i] <- var(y.modwt$LH2[,i]) } ##Get the variance from the overall y variance varLH2y <- var(varLH2ylist) #-----------------------------------------------------------------# #...

I have a 2x2 plot set up using: par(mfrow=c(2,2)) I'd like to put an overall title on the page, but I cannot figure out how. Any ideas? [[alternative HTML version deleted]]

Hi, I plotted 12 graphs on a page and output to a png file. I wanted to have an overall title for all 12 graphs. What command can I use to do this? Below is the code that plotted the 12 graphs in one page. # FM10 by Month/ Export the plot to Wash2005FM10.png png(file="Wash2005FM10.png",bg="white") par(mfrow = c(3,4)) # Plot 12 Month of OFM10, FFM10 for...

Dear R-help, I'm applying the surrogate Poisson glm, by following Venables & Ripley (7.3 pp238-42). >overall_cbind(expand.grid(treatment=c("Pema","control"),age=c("young","adult","old"),repair=c("excellent","good","poor")),Fr=c(8,0,7,1,2,0,2,7,1,4,7,1, 0,3,2,5,1,9)) >overall$age_ordered(overall$age,levels=c("young&quot...

Hi all, I know the overall display precision can be changed in R... but what about overall calculation precision? Thank you! [[alternative HTML version deleted]]

...ension ##LH2 dimLH2 <- dim(y.modwt$LH2) dimLH2x <- dimLH2[1] dimLH2y <- dimLH2[2] varLH2xlist <- c(rep(0, dimLH2x)) varLH2ylist <- c(rep(0, dimLH2y)) ##Loop to get variance from x axis for(i in seq(dimLH2x)){ varLH2xlist[i] <- var(y.modwt$LH2[i,]) } ##Get the variance from the overall x variance varLH2x <- var(varLH2xlist) ##Loop to get variance from y axis for(i in seq(dimLH2y)){ varLH2ylist[i] <- var(y.modwt$LH2[,i]) } ##Get the variance from the overall y variance varLH2y <- var(varLH2ylist) #-----------------------------------------------------------------# ##Ge...

Hi everyone, I'm trying to implement calculation of overall progress for rsync. Currently it shows progress only on per file basis. So I need for the feature size of all files to be synced - this is obvious and already exists in the code. I can calculate size of transferred files, but I can't find the piece of code that is responsible for calculation...

...0 0 0 0 10 1 -1 -1 0 1 0 0 0 -1 -1 0 0 11 1 -1 -1 0 0 1 0 0 0 0 -1 -1 12 1 -1 -1 -1 -1 -1 1 1 1 1 1 1 ... I have two questions: (1) I assume the 1st column (under intercept) is the overall mean, the 2rd column (under a1) is the difference between the 1st level of factor a and the overall mean, the 4th column (under b1) is the difference between the 1st level of factor b and the overall mean. Is this interpretation correct? (2) I'm not so sure about those interaction columns. For...

...gories and 1.0 penalty for misclassifying by 3 categories. I wanted to use this sort of penalty - but with equal distances between the 4 categories (0, 1/3, 2/3, 1). But, I found that if I simply re-scale the penalty matrix, while maintaining equal distances between categories, then the estimate of overall accuracy increases. In effect I can achieve any value for accuracy - including unity - by re-scaling the penalty matrix. So, I'd like to ask what, if any are the contraints on the scaling process? Here is a working code that illustrates my difficulty:- set.seed(1); gldstd=round(runif(200)*4);...

hey guys I want two plots in one window with an overall title and with individual titles for each plots. my code: par(mfrow=c(2,1)) bp_dirverq1=boxplot(dirverq1, col="orange",horizontal=TRUE, main="Q1 2012",cex.main=0.7) bp_dirverq2=boxplot(dirverq2, col="orange",horizontal=TRUE, main="Q2 2012",cex.main=0.7) tit...