similar to: Cutting vectors

Hello, I have an unbalanced data set with some factors and several variables. The missing values are spread unequally over the variables. Because of the unbalanced structure have i tried NLME, but I cannot understand how I should set up the na.omit. I would not like to use na.omit before setting up the model, as that would exclude too many observations (on the whole structure). Lets say I have

Hello, are there any good methods in R that will test if two correlation matrices (obtained in different ways) are equal? Something better than the Mantel test would be preferable. Regards, Patrik Waldmann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

The on-line documentation for par() refers to a 'gamma' argument. However, if 'par()' is invoked, there is no gamma argument. -- Art Salwin salwin@mitretek.org (202) 863-2985 (202) 863-2988 (fax) --please do not edit the information below-- Version: platform = Windows arch = x86 os = Win32 system = x86, Win32 status = major = 1 minor = 0.0 year = 2000 month = February

The on-line documentation for the 'ann' argument under par(), states "If set to 'FALSE'...". However, this parameter shows up as a numeric on my machine, not as a logical. -- Art Salwin salwin@mitretek.org (202) 863-2985 (202) 863-2988 (fax) --please do not edit the information below-- Version: platform = Windows arch = x86 os = Win32 system = x86, Win32

Hello, I want to fit some non linear models with the Levenberg-Marquardt algorithm. It doesn''t seem to exist any function to do this in R ( well, maybe one does, but I''m a new user, and the only documentation I have is "An introduction to R"). I''d like to know if this function exists, maybe throught an additionnal package. I''d also like to know if if

The documentation for 'order()' has the following sentence: `sort.list' is the same, using only one argument put allowing partial sorting. Should "put" be "but" or possibly "'partial'"? -- Art Salwin salwin@mitretek.org (202) 863-2985 (202) 863-2988 (fax) --please do not edit the information below-- Version: platform = Windows arch =

Is there an R function that inputs a value and a vector, and returns the percentile corresponding to the value? Something like pvalue<-function(x,vector){mean(x>=vector)} but with interpolation. Given the above function, is there a way to do the following without the loop: for(i in 1:10) y[i] <- pvalue(x[i],vector) where x and y are numeric vectors of length 10? -- Art Salwin

"An Introduction to R" notes that the function tapply() can be used to index a vector by multiple categories. But they don't show how to do this. How does one implement, for example, the example in the text in which tax accountants are indexed by both state and sex? Is there a way to use tapply() to put the vector into bins defined by the levels instead of having to exactly match

Dear Thomas & Dear List, the fitting function `coxph.fit' called by `coxph' may fail to estimate the regression coefficients when some values of the design matrix are very large. For example library(survival) ### load example data load(url("http://www.imbe.med.uni-erlangen.de/~hothorn/coxph_fit.Rda")) method <- "efron" ### copied from `coxph.fit' coxfit

Background: I administer a cluster of RedHat EWS 3U4 Linux workstations at a university. I built R 2.1.0 from source: ./configure \ --prefix=/sscc/opt/R-2.1.0 \ --with-blas=no \ 2>&1 \ | tee NUInstall.configure R is now configured for i686-pc-linux-gnu Source directory: . Installation directory: /sscc/opt/R-2.1.0 C compiler:

A number of users have spotted a terminal problem with the 'debug' package under R2.4.0, along the lines of > mtrace(x) > x() Error in attr(value, "row.names") <- rlabs : row names must be 'character' or 'integer', not 'double' This arose from a bug in 'rbind.data.frame' in R2.4.0 itself. The bug is fixed in R2.4.0 patched, so the

When attempting to merge 3 data frames, one of which has fewer columns than the others, rbind.data.frame correctly refuses to perform the bind. However, the error message given is a bit obscure due to a logical bug in the match.names() internal function to rbind.data.frame. Illustration: ## Three data frames with same column variable names: > foo <- data.frame(v1 = c('a',

Hi all, how can I remove objects from the workspace which starts with a certain pattern, e.g. lm (lm1, lm2, lm3 etc). Wildcards won?t work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

R CMD check tells me * checking R code for possible problems ... NOTE agexact.fit.rds: ... may be used in an incorrect context: ?optim(init, agfitfn, ...)? Warning: <anonymous>: ... may be used in an incorrect context: ?optim(init, agfitfn, ...)? Can anyone tell me what this message means? My searches haven't turned up anything useful. This is with R 2.7 and 2.9. The message

Dear All, I am trying to solve the optimization problem below, but I am always getting the following error: Error in optim(rep(20, nvar), f, gr, method = "L-BFGS-B", lower = rep(0, : L-BFGS-B needs finite values of 'fn' Any ideas? Thanks in advance, Paul ----------------------------------------------- k <- 10000 b <- 0.3 f <- function(x) { n <- length(x)

Hi, in R-3.0 I get the following error when calling lsfit with more observations than variables, which seems to come from an error in the formatting of the error message (note that this was not happening in 2.15.3): > nobs <- 5; nvar <- 6; lsfit(matrix(runif(nobs*nvar), ncol=nvar), runif(nobs), intercept=FALSE) Error in sprintf(ngettext(nry, "%d response", "%d

Dear R community members and R experts I am stuck at a point and I tried with my colleagues and did not get it out. Sorry, I need your help. Here my data (just created to show the example): # generating a dataset just to show how my dataset look like, here I have x variables # x1 .........to X1000 plus ind and y ind <- c(1:100) y <- rnorm(100, 10,2) set.seed(201) P <-

I'm hoping someone with some experience generating tables in Sweave will be able to solve this problem for me. I'm experiencing some inconsistency with the way floating point numbers are displayed -- for some tables, they will be in standard notation, e.g. 11±26 while in others, it uses exponential notation, e.g. 1.1e+01±2.6e+01 I am using code such as the following; note I have

Hi, I'm currently using the R package e1071 to train naive bayes classifiers and came across a bug: When the posterior probabilities of all classes are small, the result from the predict.naiveBayes function become NaNs. This is an issue with the treatment of the log-transformed probabilities inside the predict.naiveBayes function. Here is an example to demonstrate the problem (you might need

hi all i know that one should try and limit the amount of looping in R programs. i have supplied some code below. i am interested in seeing how the code cold be rewritten if we dont use the loops. a brief overview of what is done in the code. ============================================== ============================================== ============================================== 1. the input