search for: coefmat

...on(object){ if (!is.null(object$call1)){ cat("Call: ") dput(object$call1) } est <- object$qbar se <- sqrt(diag(object$t)) tval <- est/se df <- object$df pval <- 2 * pt(abs(tval), df, lower.tail = FALSE) coefmat <- cbind(est, se, tval, pval) colnames(coefmat) <- c("Estimate", "Std. Error", "t value", "Pr(>|t|)") cat("\nCoefficients:\n") printCoefmat( coefmat, P.values=T, has.Pvalue=T, si...

Hello, I'm using the R package called mvpart, which is about the multivariate regression trees. The function I wrote is: mrt1<- mvpart(coefmat~sChip+sScreen+sMem,data=mixdata, xv="pick", plot.add=TRUE,uniform=TRUE,which=4,all=TRUE,xadj=2,yadj=2,rsq=TRUE,big.pts=TRUE,wgt.ave.pca=TRUE,legend=TRUE,bars=F, pca=TRUE) where "coefmat" is a matrix(of dimension N*K) to store the coefficients of my functional data object; and t...

...:25:06 -0000 @@ -24,6 +24,7 @@ } \value{ A character vector of the concatenated values. Thus will be of length +% objects are (....) unless ??? zero if all the objects are unless \code{collapse} is non-NULL, in which case it is a single empty string. } Index: src/library/base/man/print.coefmat.Rd =================================================================== RCS file: /home/cvsroot/R/src/library/base/man/print.coefmat.Rd,v retrieving revision 1.1.1.1 diff -u -r1.1.1.1 print.coefmat.Rd --- src/library/base/man/print.coefmat.Rd 7 Nov 2003 11:57:17 -0000 1.1.1.1 +++ src/library/base/ma...

Hello, I'm using the R package called mvpart, which is about the multivariate regression trees. The function I wrote is: mrt1<- mvpart(coefmat~sChip+sScreen+sMem,data=mixdata, xv="pick", plot.add=TRUE,uniform=TRUE,which=4,all=TRUE,xadj=2,yadj=2,rsq=TRUE,big.pts=TRUE,wgt.ave.pca=TRUE,legend=TRUE,bars=F, pca=TRUE) where "coefmat" is a matrix(of dimension N*K) to store the coefficients of my functional data object; and t...

R-helpers: I have a dataset that has 168 subjects and 12 variables. Some of the variables have missing data and I want to use the multiple imputation capabilities of the "mice" package to address the missing data. Given that mice only supports linear models and generalized linear models (via the lm.mids and glm.mids functions) and that I need to fit Cox models, I followed the previous

>>>>> Shu Fai Cheung >>>>> on Thu, 6 Jul 2023 17:14:27 +0800 writes: > Hi All, > I would like to ask two questions about printCoefmat(). Good... this function, originally named print.coefmat(), is 25 years old (in R) now: -------------------------------------------------------------------- r1902 | maechler | 1998-08-14 19:19:05 +0200 (Fri, 14 Aug 1998) | Changed paths: M R-0-62-patches/CHANGES M R-0-62-patches/s...

...n be given to students if you are worried about their minds TL> being contaminated by stars. I'm going to do this on the TL> installation here. TL> TL> Perhaps we should have an option to leave p-values off summary.{lm,glm}? KH> I thought we had this. print.coefmat() uses .Options$signif.stars ... >> PD> Note: *p-values*, not stars. Not a good idea, in my opinion. nor in mine. TL> I agree it's a bad idea to leave p-values off, but I still think TL> it's reasonable to have *an option* to leave them off. Some people...

Something for the wish list (not really a bug): I was reminded of what I see as a problem with summary.glm last week when some of my students fell into a trap in a homework exercise, defining a logit model in which coefficients were aliased. When this happens in lm, summary.lm prints a message ["Coefficients: (1 not defined because of singularities)"], but summary.glm is silent. In

...3 1474.65 491.55 8.7852 0.000376 *** Residuals 25 1398.80 55.95 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 This is what I expected. Under RW 1.0.0 I get Analysis of Variance Table Response: kaeufe Error in print.coefmat(x, digits = digits, signif.stars = signif.stars, : binary operator applied to invalid types Am I doing something wrong or is this a real bug? Regards, *** D.Trenkler *** ==================================================================== Dietrich Trenkler (trenkler...

Hello R-community, It's been a week now that I am struggling with the implementation of a cox model in R. I have 80 cancer patients, so 80 time measurements and 80 relapse or no measurements (respective to censor, 1 if relapsed over the examined period, 0 if not). My microarray data contain around 18000 genes. So I have the expressions of 18000 genes in each of the 80 tumors (matrix

Hi All, I would like to ask two questions about printCoefmat(). First, I found a behavior of printCoefmat() that looks strange to me, but I am not sure whether this is an intended behavior: ``` r set.seed(5689417) n <- 10000 x1 <- rnorm(n) x2 <- rnorm(n) y <- .5 * x1 + .6 * x2 + rnorm(n, -0.0002366, .2) dat <- data.frame(x1, x2, y) out <-...

...f the system Rprofile file is skipped BDR> and hence options show.signif.stars and show.coef.Pvalues are not set. >> but other options as well.. >> I think this shouldn't happen and *is* a bug. BDR> It looks to me as if a bit more protection is needed in print.coefmat, BDR> at least. which Brian has implemented in the mean time. .... >> For the 1.0.x patches, at least, I'll move the autoload()'s to the end of >> profile/Common.R BDR> I _think_ they should have try() around them, which is what S-PLUS BDR>...

...t; Yes, in that the rest of the system Rprofile file is skipped BDR> and hence options show.signif.stars and show.coef.Pvalues are not set. but other options as well.. I think this shouldn't happen and *is* a bug. BDR> It looks to me as if a bit more protection is needed in print.coefmat, BDR> at least. In extrapolation, wouldn't all R functions using non-internal *) options() {i.e. getOption() nowadays} need protection against these options being wrong or non existent ? I feel that this is exaggerated.. *) Currently, the following options are "internal" (i...

...ot; "format.char" "format.default" "format.pval" [13] "help.search" "hist.default" "legend" "ls.print" [17] "print.aov" "print.aovlist" "print.coefmat" "print.dummy.coef" [21] "print.glm" "print.glm.null" "print.htest" "print.lm" [25] "print.mtable" "print.summary.glm" "print.summary.lm" "print.summary.lm.null" [29] &quot...

...of the system Rprofile file is skipped > BDR> and hence options show.signif.stars and show.coef.Pvalues are not set. > but other options as well.. > I think this shouldn't happen and *is* a bug. > > BDR> It looks to me as if a bit more protection is needed in print.coefmat, > BDR> at least. > In extrapolation, wouldn't all R functions using non-internal *) > options() {i.e. getOption() nowadays} need protection against these options > being wrong or non existent ? > I feel that this is exaggerated.. Yes, they should, and when I write them I...

...hen adds row names o codes() now distinguishes between ordered and unordered factors o codes() had *opposite* semantics of Splus. Now it's the same. o replicating factors now yields factors (again) o print.summary.xxx functions more consistent, using new function print.coefmat(). o Changed many `T' to `TRUE' and `F' to `FALSE' in the base package. o binary operation on 1x1 matrix lost dimension o fix anova.glm for null model o glm.fit.null: ensure df.residual == df.null for a null model o summary.glm: correlations in s...

...hen adds row names o codes() now distinguishes between ordered and unordered factors o codes() had *opposite* semantics of Splus. Now it's the same. o replicating factors now yields factors (again) o print.summary.xxx functions more consistent, using new function print.coefmat(). o Changed many `T' to `TRUE' and `F' to `FALSE' in the base package. o binary operation on 1x1 matrix lost dimension o fix anova.glm for null model o glm.fit.null: ensure df.residual == df.null for a null model o summary.glm: correlations in s...

I think I've discovered what went wrong. My workspace included a function wilcox.test formerly copied from ctest. Now ctest ist part of the distribution and because of that I always got the message: [Previously saved workspace restored] Error in autoload("wilcox.test", "ctest") : Object already exists I didn't take care of this message. Now I removed

In RW0.651 and RW0.901, I discovered some unexpected behaviour when I used as.matrix.data.frame() > matrix('"', 2, 2) [,1] [,2] [1,] "\"" "\"" [2,] "\"" "\"" > unclass(as.data.frame(matrix('"', 2, 2))) $V1 [1] "\"" "\"" $V2 [1] "\"" "\""

Dear R Folks, I'm a big fan of R, but there are a couple of things that repeatedly annoy me, and I wondered if anyone has neat ways to deal with them. (a) When using "apply" row-wise to a matrix, it returns the results column-wise, and to preserve the original orientation, I've to do a transpose. E.g. I've to keep doing a transpose, which I consider to be quite