search for: ynames

I have sort of an emergency question for the list. One of my professors for an S-Plus intensive class distributed a function to produce partial regression plots. I need to run it under R, because I'm doing the homework on my home computer with a modem; hence I don't have the speed required to emulate X-Windows and run S Plus off one of the campus servers. Bottom line: I'm using R.

...glm.fit <- 2 function (x, y, weights = rep(1, nobs), start = NULL, etastart = NULL, 3 mustart = NULL, offset = rep(0, nobs), family = gaussian(), 4 control = glm.control(), intercept = TRUE) 5 { 6 x <- as.matrix(x) 7 xnames <- dimnames(x)[[2]] 8 ynames <- names(y) 9 conv <- FALSE 10 nobs <- NROW(y) 11 nvars <- NCOL(x) 12 if (nvars == 0) { 13 cc <- match.call() 14 cc[[1]] <- as.name("glm.fit.null") 15 return(eval(cc, parent.frame())) 16 } 17 if (is.null(...

...As I didn't get any answer to my bug-report last week I have taken the effort and extracted a minimal data set from my data (see below) where the following bug occurs: > glm(SKR.ein.aus ~ ., family = binomial, data = bugdata, na.action = na.omit) Error in names<-.default(*tmp*, value = ynames) : names attribute must be the same length as the vector In addition: Warning messages: 1: fitted probabilities of 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, 2: fitted probabilities of 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.f...

Dear R-team There seems to be a bug in glm.fit - I got the following error message: > > > + Error in names<-.default(*tmp*, value = ynames) : names attribute must be the same length as the vector In addition: Warning messages: 1: fitted probabilities of 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, 2: fitted probabilities of 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.f...

Hello, I've been looking through ?phantom and ?expression and this forum for examples of how I might be able to manipulate some of the names that appear on the y-axis of the barplot below. For example, the "gw" in "ECgw" would appear as a subscript...or "qr" would be the theta symbol followed by subscript "r". My attempts haven't even come close

Is there a way to create a 'bubble plot' in R? For example, if we define the following data frame containing the level of y observed for 5 patients at three time points: time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5))

Dear R users, I am looking for the source code of the implementation of ppr (Projection Pursuit Regression) in R. It will be great if citations of the source papers on which the implementation is based, are also provided. Thank you, Arvind Iyer, Grad student, Deptt. of Biomedical Engineering Viterbi School of Engineering University of Southern California, Los Angeles [[alternative HTML

Hi, I want to test for a contrast from a regression where I am regressing the columns of a matrix. In short, the following. X <- matrix(rnorm(50),10,5) Y <- matrix(rnorm(50),10,5) lm(Y~X) Call: lm(formula = Y ~ X) Coefficients: [,1] [,2] [,3] [,4] [,5] (Intercept) 0.3350 -0.1989 -0.1932 0.7528 0.0727 X1 0.2007 -0.8505 0.0520

I'm not understanding something. I'm trying to add xlab & ylab to a balloon plot of a table object. From docs I thought following should work: require(gplots) # From balloonplot example: # Create an example using table xnames <- sample( letters[1:3], 50, replace=2) ynames <- sample( 1:5, 50, replace=2) tab <- table(xnames, ynames) balloonplot(tab) # Try xlab, ylab: balloonplot(tab, xlab = "MyX", ylab = "MyY") But second plot is no different from first. R.version.string: Version 2.3.1 (2006-06-01) gplots version: 2.3.0 on Win...

Hello, The e1071 function naiveBayes gives an error when called by do.call, although it doesn't give any error if I call it manually. Browse[1]> trainParams at classifier function (x, ...) UseMethod("naiveBayes") <environment: namespace:e1071> Browse[1]> trained <- do.call(trainParams at classifier, paramList) Error in names(dimnames(tables[[i]])) <- c(Yname,

Over the years there have been inquiries to these lists as to how to produce plots with the direction(s) of axes reversed. In response to some of these inquiries, I have provided a function call ``rplot'' which automated the procedure for such axis reversal. Recently Herberto Ghezzo of McGill University pointed out to me that there was a bug in my code. I have now fixed this bug (and a

...> X <- matrix(scan("miX_e.txt", 0), ncol=2, byrow=TRUE) Read 450 items > Y <- matrix(scan("g1_sr_e.txt", 0), ncol=1, byrow=TRUE) Read 225 items > g1.ppr <-ppr(X,Y, nterms=3, max.terms=5) Error in matrix(smod[q + 6 + p * ml + 1:(q * mu)], q, mu, dimnames = list(ynames, : length of dimnames[1] not equal to array extent Where "miX_e.txt" is a file containing 225 rows (225 cases for training) with 2 columns (2 explanatory variables), and "g1_sr_e.txt" is a file with one column (1 response) and 225 rows (one for each of the 225 trining...

...mlm() fails. This works well: > y <- matrix(rnorm(20),nrow=10) > x <- matrix(rnorm(10)) > mod1 <- lm(y~x) > summary(mod1) ... But this does not: > f <- y~x > mod2 <- lm(f) > summary(mod2) Error en object$call$formula[[2L]] <- object$terms[[2L]] <- as.name(ynames[i]) : objeto de tipo 'symbol' no es subconjunto I would say that the problem is in the following difference: > class(mod1$call$formula) [1] "call" > class(mod2$call$formula) [1] "name" As far as I understand, summary.mlm() creates a list of .lm objects from the...

...object, ...) { .local <- function (object, standardize = FALSE) { m <- ncol(object at X) n <- nrow(object at X) resl <- lapply(object at models, residuals, standardize = standardize) resm <- matrix(c(unlist(resl)), ncol = m, nrow = n) ynames <- paste("y", 1:2, sep = "") colnames(resm) <- ynames rownames(resm) <- rownames(object at X) return(resm) } .local(object, ...) } <environment: namespace:gogarch> Signatures: object target "GoGARCH" defined &...

I think this is another one of the same kind of bugs in ts: Version 0.64.1 (May 8, 1999) ... > z <- ts(matrix(1:20,10,2), start=c(1969,1), frequency=12) > (z > 5) | (z < 2) Error: invalid time series parameters specified > Paul -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey<-c(0,0,0,1,1,1) Wax<-c(1,0,1,0,1,0) f<-function(x,y){ table(x,y) } f(Monkey,Wax) so that the printout is not y x 0 1 0 1 2 1 2 1 but Wax Monkey 0 1 0 1 2 1 2 1 -- View this message in context:

R2.1.1 Win 2k I have a function, B, within a function, A. I would like to have B print the name of the argument passed to it (not the value of the arguments). i.e., A<-function() { B<-function(x,y) { fit1<-lm(y~x,data=jo) print(summary(fit1) I want B to print the string "age" and the string "height". } B(age,height) }

To the purpose of fitting a 2nd order polynomial (a + b*x + c*x^2) to the chunk of signal falling in a 17 consecutive samples window I wrote the following very crude script. Since I have no previous experience of using Least Square Fit with R I would appreciate your supervision and suggestion. I guess the returned coefficients of the oolynomial are: a = -1.3191398 b = 0.1233055 c = 0.9297401

Hi all, I have a data.frame with the following colnames pattern: x1 y11 x2 y21 y22 y23 x3 y31 y32 ... I.e. I have an x followed by a few y's. What I would like to do is turn this wide format into a tall format with two columns: "x", "y". The structure is that xi needs to be associated with yij (e.g. x1 should next to y11 and y12, x2 should be next to y21, y22, and

Hi Is it possible to plot values against strings, i.e. the x axis has people's names and the y axis as values, as shown below: 5 x 4 x x 3 x 2 x 1 0 Dave Tim Ian Steve Paul Sorry if this is a little primitive, but I think it gets the point across. regards Danny