similar to: How to estimate Variance Components

Hello, I would like to know if R provides the similar function to produce the 4-in-1 graph ("Normal Plot of Residuals", "I Chart of Residuals", "Histogram of Residuals" & "Residuals vs Fits") as MiniTab. If do not have, could anyone please tell me how to produce above-mentioned individual chart? Your help is greatly appreciated.

Hi, RWeb web site is down past couple days, I am insterested in this project and want to try it for my projects. Deos anyone have this installation and user guide? Thanks! Best Regards, WeiQiang Li

Hi ALL, I am trying to compute adjusted SS & estimated component variance in GLM with un-balanced data using R. Can anyone advise me? Thanks in advance. Best Regards, WeiQiang Li [[alternative HTML version deleted]]

Dear friends, I am trying to create a web application to produce some statistical result using R. In order to avoid high CPU usage of web server caused by R, I have to create an ASP.NET web service in another server to involve R. But I am facing the unauthorizedAccessException when I call the web service, even I have assigned access and launch permission to everyone using

Hello, I am a newbie in R project and trying to call prcomp(x) of R function using (D)COM server communicate with R in ASP, and encountering the error "Runtime error -2147221493(8004000b). Automation Error, Object is static, operation not allowed." Source code is shown as below: <% Set StatConn=Server.CreateObject("StatConnectorSrv.StatConnector")

Hi, I am facing the problem that I want to plot a histogram chart set freq to true and overlay with normal or weibull or exponential distribution curve. The sample code is shown as below: >samp<-c(-8.2262,-8.2262,-8.2262,-8.20209,-8.09294,-8.07321,-8.07321, -8.07321,-8.07175,-8.04948,-8.04948,-8.04948,-8.03848,-8.03848,

Hello, and please excuse this off-topic question, but I have not been able to find an answer elsewhere. Consider a value Z that is calculated using the product (or ratio) of two means X_mean and Y_mean: Z=X_mean*Y_mean. More generally, Z=f(X_mean, Y_mean). The standard error of Z will be a function of the standard errors of the means of X and Y. I want to calculate this se of Z. Can someone

Suppose I have x<-data.frame(v1=1:4, v2=c(2,4,NA,7), v3=rep(1,4), v4=LETTERS[1:4],v5=rep('Z',4)) or a much larger frame, and I wish to test for and remove the constant numeric columns. I made: is.constant<-function(x){identical(min(x),max(x))} and apply(x,2,is.constant) # Works for numerics x[,-which(apply(x,2,is.constant))] I'd really like to be able to

Hello, I am a newbie in R project and trying to call prcomp(x) of R function using (D)COM server communicate with R in ASP, and encountering the error "Runtime error -2147221493(8004000b). Automation Error, Object is static, operation not allowed." at line Result=StatConn.Evaluate("y<-prcomp(x)"). E ven I have added the line StatConn.EvaluateNoReturn("library

Hola! I ask here since I learnt from this list that the LaTeX package listings should be good for typesetting R code. I encountered one problem: \begin{lstlisting} X %*% V \end{lstlisting} in the output the * in %*% disappears! same with %/%, etc, the / disappears. Any ideas? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. --

Hello R'ers In previous versions of R (I now use 2.0.1) there was a package "bootstrap". I wrote some programs that depended heavily on this package but can unfortunately not find it on Cran today. Is there any way to find an older version of "bootstrap" and use it with the new version - 2.0.1? Sincerely Fredrik Lundgren Norrk??ping Sweden

Dear friends, I am trying to produce Contour Plot with R, but there are some NA in my data matrix. After I ran the following R script, I got the error message:"no proper `z' matrix specified". Does anybody know how to plot contour chart with R for the non-strict matrix? Thank you in advance!!!

Hi, Would you help me solve the following question? Thanks. Question: If I try to set the probability=0.05 and find the approximate x. (The answer should be somewhere between 2.1782 and 2.1783.) I write about this R program as follows but I don¡¦t know how to get the value of x which is between 2.1782 and 2.1783. library(mvtnorm) value<-array(1000) a<-array(1000)

Hi, I want to estimate a VAR-model and calculate the impulse response function and a variance decomposition. I am familiar with standard R-functions, like e.g. arima. But I have not found equivalent functions to estimate a VAR-modell. Are there any functions available or which R-package can I use to solve my problem? Thank you for your help. Marc Gronwald University of Hamburg Department of

How might you fit a generalized linear model (glm) with variance = mu+theta*mu^2 (where mu = mean of the exponential family random variable and theta is a parameter to be estimated)? This appears in Table 2.7 of Fahrmeir and Tutz (2001) Multivariate Statisticial Modeling Based on Generalized Linear Models, 2nd ed. (Springer, p. 60), where they compare "log-linear model fits to

On trying with very smll test packages to try drilling down to som problems earlier reported in this list, I detected the following: R CMD check on a package with an R subdirectory empty, reports: * checking R files for syntax errors ... ERROR Syntax error in file On removing the empty subdirectory, the error disappears. Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of

I found the need of converting a matrix into a list of its columns (for use with do.call), and was surprised there was no method as.list.matrix, it could easily be a part of as.list default I wrote my.as.list.matrix <- function(mat) { if(!is.matrix(mat))stop("Argument must be a matrix") n <- NCOL(mat) res <- vector(mode="list", length=n)

I am trying to install a local package and get this unexpected error: ---------- Making package UMSA ------------ adding build stamp to DESCRIPTION installing R files installing data files installing man source files installing indices Error: couldn't find function "na.omit" Execution halted na.omit of course is in package stats, and that is listed in the Depends field in

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil --

I am trying to run a simple nested manova with two levels of nesting, Sires, and Dams within Sires. The goal is to extract the among sires covariance matrix and secondarily, the among Dams within Sires covariance matrix. Both sires and dams are random effects. At present there are four dependent variables, but that may change. This is part of a larger bootstrapping project, so efficiency and