similar to: printCoefmat with more p-values?

Dear useRs, I need to compute zero of polynomial function fitted by lm. For example if I fit cubic equation by fit=lm(y~x+I(x^2)+i(x^3)) I can do it simply by polyroot(fit$coefficients). But, if I fit polynomial of higher order and optimize it by stepAIC, I get of course some coefficients removed. Then, if i have model y ~ I(x^2) + I(x^4) i cannot call polyroot in such way, because there is

Dear useRs, Is there any more elegant way to convert dataframe to a vector of all its values than as.vector(as.matrix(x)) ? I did not have to do such conversion yet, so I am not sure... (of course as.vector() alone does not work). Regards, -- Lukasz Komsta Department of Medicinal Chemistry Medical University of Lublin Jaczewskiego 4, 20-090 Lublin, Poland Fax +48 81 7425165

Dear useRs, I was searching CRAN for implementation of kurtosis and skewness tests, and found that there is some kind of lack on it. So, I have written three functions: 1. Anscombe-Glynn test for kurtosis 2. Bonett-Seier test based on Geary's kurtosis (which is not widely known, but I was inspired by original paper describing it, found coincidentally in Elsevier database) 3.

Hi, I'm using a fairly simple HP Compaq desktop PC running Windows 2K. When running a large process in R, the process "RGUI.exe" will never exceed 50% of the CPU usage. The program used to be able to use more of the computer, but does not now. I don't believe this is a multiple processor machine. Can anyone give any advice on how to solve the problem? Thanks, Michael

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 <- lm(y ~ x1 + x2, dat) out_summary <-

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

Dear useRs, I have a curve which is a mixture of Gaussian curves (for example UV emission or absorption spectrum). Do you have any suggestions how to implement searching for optimal set of Gaussian peaks to fit the curve? I know that it is very complex problem, but maybe it is a possibility to do it? First supposement is to use a nls() with very large functions, and compare AIC value, but it is

Hi List, My goal is to force R not to print in scientific notation in the sixth column (rel_diff - for the p-value) of my data frame (not a matrix). I have used the format.pval () and printCoefmat () functions on the data frame. The R script is appended below. This issue is that use of the format.pval () and printCoefmat () functions on the data frame gives me the desired results, but coerces

printCoefmat(signif.legend =FALSE) does not work properly. The option "signif.legend = FALSE" is ignored as shown in the example below. cmat <- cbind(rnorm(3, 10), sqrt(rchisq(3, 12))) cmat <- cbind(cmat, cmat[,1]/cmat[,2]) cmat <- cbind(cmat, 2*pnorm(-cmat[,3])) colnames(cmat) <- c("Estimate", "Std.Err", "Z value", "Pr(>z)") #

On 29/10/2023 3:48 a.m., Shu Fai Cheung wrote: > Hi all, > > Just a minor issue that I am not sure whether this is considered a > "bug." It is about the help page. > > In the help page of printCoefmat(), for the argument 'eps.Pvalue', the > description is as below: > > number, .. > > I have to read the source to figure out that this argument

Dear All, I have small samples of data (between 6 and 15) for numerious time series points. I am assuming the data for each time point is normally distributed. The problem is that the data arrvies sporadically and I would like to detect the number of outliers after I have six data points for any time period. Essentially, I would like to detect the number of outliers when I have 6 data points then

Hi all, Just a minor issue that I am not sure whether this is considered a "bug." It is about the help page. In the help page of printCoefmat(), for the argument 'eps.Pvalue', the description is as below: number, .. I have to read the source to figure out that this argument is to be used by format.pval(). Maybe the description of 'eps.Pvalue' can be revised to refer

I want to use something like printCoefmat() in a print.summary method to print a more nicely formatted version of the result from a summary method, but where the estimates may be cross-classified by one or more factors. However, printCoefmat() assumes that the labels for the parameters are the rownames of the object, and prints factors as integers. With one factor, I could just assign that as

I want to print the coefficient estimates of a model in a way as consistent with other output in R as possible. stats provides the printCoefmat function for doing this, but there is one problem. I have an additional piece of textual information I want to put on the line with the other info on each coefficient. The documentation for printCoefmat says the first argument must be numeric, which

I'm trying to understand how R prints summary.lm objects and trying to change it slightly for a summary function that calculates standard errors using an alternative method. I've found that I can modify a summary.lm object and then it prints the modified way but I want to change a few things in the print method that I think I might just be able to do. One is that I want the coefficients

Dear R experts, I am posting this question on behalf of a Japanese R user who wants to know how to change the siginificant codes default. As you know, R's default significant codes are: Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 But he says that it is usual in economics to give codes such as `***' for 0.01, `**' for 0.05 and `*' for

I would like to extend Dixon's values beyond 30. I've read over Rorabacher article but didn't understand the equations well enough to convert them to Excel and then "drag" the cells out extending the n. Rorabacher,_1991.pdf <http://r.789695.n4.nabble.com/file/n4649819/Rorabacher%2C_1991.pdf> Dean_&_Dixon,_1951.pdf

Hello, I'd like to add custom labels to my pair() plot. These labels include math expression but they aren't correctly displayed... Further, I want that the boxes for the text.panel (diagonal) have an other background color (grey80). Is that generally possible? If yes how do I have to set it? What I've so far is: panel.cor <- function(x, y, digits=2, prefix="",

Hello everyone, I'm trying to estimate the parameters of the returns series attached using the GARCH code below, but I get the following error message: Error in solve.default(Hessian) : system is computationally singular: reciprocal condition number = 0 Error in diag(solve(Hessian)) : error in evaluating the argument 'x' in selecting a method for function 'diag' Can

I've defined the following for objects of a class called jml summary.jml <- function(object, ...){ tab <- cbind(Estimate = coef(object), StdError = object$se, Infit = object$Infit, Outfit = object$Outfit) res <- list(call = object$call, coefficients = tab, N = nrow(object$Data), iter = object$Iterations) class(res) <- "summary.jml" res }