similar to: R 2.1.0 RH Linux Built from Source Segmentation Fault

Hi, I tried to use rmnlIndepMetrop (bayesm package) for my MNL model with 4 choice alternatives, 5 independent variables, 69 observations, dim(X) [1] 276 5, nu=6. So I run such code: if(nchar(Sys.getenv("LONG_TEST")) != 0) {R=2000} else {R=10} set.seed(66) df=read.table("X_metrop.dat",header=TRUE) inp=as.matrix(df) y=as.numeric(inp[,1]) n=length(y) p=4

Hi, I am running the following code: mfg0 <- ggplot(aes(x=Grade,y=Math,colour=RiskStatic45678),data=math.f) mfg1 <- mfg0 + geom_smooth(method="lm", formula=y ~ ns(x,2),size=1) + geom_smooth(aes(y=nalt.math,color="NALT"),size=1,data=nalt) + scale_colour_brewer("Risk Status", pal="Set1") + coord_cartesian(ylim = c(175, 245)) mfg1 And when I specify

Hi, I have the following data about courses (504) in a university, two attributes about the proportion of resources used (#resources_used / #resources_available), namely the average and the standard deviation. Thus I have: [1] n=504 rows [2] 1 id column and 2 attributes Here's a sample of the data: courseid,average,std 12741,1,0 17161,1,0 12514,1,0

Error in `[.data.frame`(x, order(x, na.last = na.last, decreasing = decreasing)) : undefined columns selected during the execution of following r sequence of commands X<-subset(data,select=c(V1,V2,V3,V4,V5,V6,V7,V8,V9)) y<-subset(data,selcet=10) Data = list(y=y, X=X, p=.75) Prior = list(betabar=c(rep(0,ncol(X))),A=.01*diag(ncol(X))) Mcmc = list(R=100000, keep=10, step=.2) out <-

Hi R-users, I have the following problem: I am trying to learn something about bayes methodology and started paying around bayesm package, but could not replicate the Conjunctive model's estimates as they appear in Rossi et al "Bayesian Statistics and Marketing", 2005, JWS, pages 264-265, Table CS4.4. I have downloaded in my working directory the documents from

Hello, I have built a multi-variate probit model using the package "bayesm", which requires that the X data is constructed using the function "CreateX". I've gone through the documentation and run my model, but wanted to be sure about my interpretation of the results for the coefficients - beta. Steps: 1) I have 5 choices for the dependent variable Y, so p=5 2) I have 8

Hello, I have built a multi-variate probit model using the package "bayesm", which requires that the X data is constructed using the function "CreateX". I've gone through the documentation and run my model, but wanted to be sure about my interpretation of the results for the coefficients - beta. Steps: 1) I have 5 choices for the dependent variable Y, so p=5 2) I have 8

We are currently working to trunk from a Nortel 81C to an Asterisk Server 1.4 running on Red Hat Linux. We have two PRI trunks which work with the exception of the clock slips, which is causing the Nortel to reset the PRIs once a hour. Thanks for any suggestions. 81C MSDL Asterisk Digium TE110P REQ prt TYPE adan dch 10

Hi, I am looking for a way in which I can check if rows in 1 dataframe are present in another data frame in a unique way. A row in dataframe should be super set of any row in another dataframe. I can write a for loop for it, however, that will be inefficient. So, I am looking for an efficient way to do this in R. I have explained it with an example below: I want to check if a row in dataframe

Hi, I'm new to this mailing list so apologies if this is too basic. I have confocal images 512x512 from which I have extracted x,y positions of the coordiates of labelled cells exported from ImageJ as a.csv file. I also have images that define an underlying pattern in the tissue defined as areas of different pixel values 0 or 255 (also 512x512) I've exported these images as .txt

hi all, as i heard, i need winbind to connect linux-clients on samba. i installed this - without pam (in the hope this would work...) but this constellation seems to be needed, as i still cannot connect, even if winbind logs my domain as accepted. Do i really have to do the pam-configuration on the samba-server-side too?? i just thought, it was enough to authentificate against smbpasswd. is

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

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

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

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

Hi, I would like to connect my GHA-system to R and to deliver the library to users of R. could you give me a comprehensive example containing the step by step procedure to: 1. create a callable library with R binaries (not having a standalone version of R with its own main()). If possible, I would like to run the system on unix,alpha,linux computers as well as on IBM's parallel

Hi, I sent the below message to r-devel but received no answer. I would like to contribute my system but I need som advice concerning the connections. Is my proposal of any interest to you? regards Ralf ?stermark ----- Forwarded message from Ralf ?stermark <rosterma at abo.fi> ----- Date: Mon, 21 Feb 2005 11:35:21 +0200 From: Ralf ?stermark <rosterma at abo.fi> Reply-To: Ralf