similar to: Error in solve.default

R-Listers: I am doing a quasi-maximum likelihood estimation and I get a "subscript out of bound" error message, Typically I would think this means that a subscript used in the function is literally out of bounds however I don't think this is the case. All I change in the code is a constant, that is hard-wired in (not data dependent and not parameter dependent), furthermore,

I am replying to the email below my signature; it is nominally about dmp files, but apparently specific to problems with Bill DuMouchel's hblm program. I have been using this program on Splus for years, and ran into some of the same problems Jack did in trying to get it into R. I talked with Bill about it, and he didn't know that there would be any problems; he assumed that R was

---------- Forwarded message ---------- From: <r-help-bounces@r-project.org> Date: Mon, May 11, 2009 at 10:24 PM Subject: The results of your email commands To: das.moumita.online@gmail.com The results of your email command are provided below. Attached is your original message. - Results: Ignoring non-text/plain MIME parts - Unprocessed: What is the function for partial

Hi, I was very much hoping someone could help me with the following. I am trying to convert some SAS NLMIXED code to NLME in R (v.2.1), but I get an error message. Does anyone have any suggestions? I think my error is with the random effect "u" which seems to be parametrized differently in the SAS code. In case it's helpful, what I am essentially trying to do is estimate parameters

Hi, A error message arose while I was trying to fit a ordinal model with lrm() I am using R 2.8 with Design package. Here is a small set of mydata: RC RS Sex CovA CovB CovC CovD CovE 2 1 0 1 1 0 -0.005575280 2 2 1 0 1 0 1 -0.001959580 2 3 0 0 0 1 0 -0.004725880 2 0 0 0 1 0 0 -0.005504850 2 2 1 1 0 0 0 -0.003880170 1 2 1 0 0 1 0 -0.006074230 2 2 1 0 0 1 1 -0.003963920 2 2 1 0 0 1 0

Just noticed nobody has commented on this, so I took a look, and I can't see the problem (2.1.0 looks the same as 2.1.1 to me). What is supposed to be wrong? Duncan Murdoch On 5/23/2005 2:35 PM, pelzc at ohsu.edu wrote: > Full_Name: Carl Pelz > Version: 2.1.0 & 2.0.1 > OS: Windows 2000 > Submission from: (NULL) (137.53.64.158) > > > # Paste this all into an R

Hi WizaRds, I tried reading the S-PLUS file ftp://ftp.research.att.com/dist/bayes-meta/hblm.dmp into R using data.restore("hblm.dmp") but I got an error: Error in attributes(value) <- thelist[-match(c(".Data", ".Dim", ".Dimnames", : row names must be 'character' or 'integer', not 'double' In

Hi I'm encountering an error/warning when doing multiple comparisons with the package multcomp on a coxme model. My data: I'm looking at the removal of brood from the nest according to three treatments I applied on the brood. The brood and the workers caring about the brood in the nest, belonged to different colonies. Factor: treatment (3 levels: tx,uv,meta) Random effect 1: origin of

Dear All I am trying to estimate a mixture of regression and get the following error using the mixreg package: Error in y - yhat : non-conformable arrays The instruction I used were: x <- as.matrix(LRHUN) y <- as.matrix(LRINTER) TS <- list(list(beta=c(3.0,1.0),sigsq=1,lambda=0.4), list(beta=c(0.0,1.0),sigsq=1,lambda=0.6)) prova <- mixreg(x,y, ncomp=2, theta.start=TS)

Hello, struggling with the very basic needs... :( any help appreciated. #using the package doBY #who drinks how much beer per day and therefor cannot calculate rowise maxvals evaluation=data.frame(date=c(1,2,3,4,5,6,7,8,9), name=c("Michael","Steve","Bob", "Michael","Steve","Bob","Michael","Steve","Bob"),

Hi I need help in the recursive problem. this is my code #Generate two random Numbers minval=20 maxval=100 cutoffValue=50 optVal<- function(cutoffValue,minval,maxval) { x=runif(2) x=x*cutoffValue for( i in 1:2) { if(x[i] < 30 || x[i] >60) # checking it falls between the range { optVal(cutoffValue,minval,maxval) } } return(x) }

hi, I have been attempting to run this script and am getting some strange results. The script connects to a database and retrieves a series of tables, using sequential sql statements. I have tested all of the sql statements in the PostGreSQL terminal and they all return the desired results. I place each table into a list and run a FOR loop for 'i' in the list. The script generates

Thank you Brian Ripley, Gardar Johannesson, and Marcel Wolbers for your prompt and friendly help! I will share any further learnings as I move through these suggestions. -Bob Abugov Brian Ripley wrote: I just use optim() on the sum of squares with non-negativity constraints. That did not exist in 1999. Gardar Johannesson wrote: You can always just use the quadratic programing library in R

Dave Langer in this video https://www.youtube.com/watch?v=z8PRU46I3NY uses the titanic data as an example of using caret to create xgbTree models. The caret train() function has a tuneGrid parameter which takes a list set up like so: tune.grid <- expand.grid(eta = c(0.05, 0.075, 0.1), nrounds = c(50, 75, 100), max_depth = 6:8,

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

Hello, I am writing an extension that requires the user to pass a function into the C code. The function returns a covariance matrix. The user defines an arbitrary set of parameters as the function's arguments, which are used to construct this covariance matrix. Within the C code these parameters and hence covariance matrix are iteratively updated. The user's function might look like

On Thu, Sep 3, 2009 at 6:24 AM, Olivier Meurant<meurant.olivier at gmail.com> wrote: > Hi Tanya, > > I have tried the 2.6 pre-release on the following host : > > Windows XP pro SP2 with mingw/msys : > > uname -a > MINGW32_NT-5.1 OLIVE 1.0.10(0.46/3/2) 2004-03-15 07:17 i686 unknown > > gcc -v > Reading specs from d:/mingw/bin/../lib/gcc/mingw32/3.4.5/specs

Dear r-helpers, I have a question about multiple testing. Here an example that puzzles me: All matrixes and contrast vectors are presented in treatment contrasts. 1. example: library(multcomp) n<-60; sigma<-20 # n = sample size per group # sigma standard deviation of the residuals cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE, dimnames =

I have ported some R code to C to make it faster. I can perform .Call("foobar",....) once and it works fine. Absolutely correct answer. If I put a loop inside foobar and run the main code routine more than 100 times, it crashes R. Or if I call .Call("foobar"....) seperately more than two tims it crashes R. For the most part I am doing matirx multiplies using EXP

I am 100% certain that there is an easy way to do this, but after experimenting off and on for a couple of days, and searching everywhere I could think of, I haven't been able to find the trick. I have this piece of code: ... attach(d) if (ORDINATE == 'ds') { lo <- loess(percent ~ ncms * ds, d, control=loess.control(trace.hat = 'approximate')) grid <-