similar to: odd behaviour of 'while'?

Dear R users, I want to conduct a small simulation study and I have to use the lme function in a loop to save the restricted log likelihood. However, for one simulated data set the lme function gives this error Error en lme.formula(yboot ~ X[, -1], data = data.fr, random = Z.block) : nlminb problem, convergence error code = 1 message = singular convergence (7) and then, the

I am using the 'glmmPQL function in the 'MASS' library to fit a mixed effects logistic regression model to simulated data. I am conducting a series of simulations, and with certain simulated datasets, estimation of the random effects logistic regression model unexpectedly terminates. I receive the following error message from R: Error in lme.formula(fixed=zz + arm.long,random=~1 |

Hello Helpers, I have some problems with fitting the model for my data... -->my Literatur says (crawley testbook)= Non-normality of errors-->I get a banana shape Q-Q plot with opening of banana downwards Structure of data: origin wt pes gender 1 wild 5.35 147.0 male 2 wild 5.90 148.0 male 3 wild 6.00 156.0 male 4 wild 7.50 157.0 male 5 wild 5.90

nlminb( obj = function(x) x, start=1, lower=-Inf, upper=Inf ) $par [1] 0 $objective [1] 0 $convergence [1] 0 $message [1] "absolute function convergence (6)" $iterations [1] 1 $evaluations function gradient 2 2 [[alternative HTML version deleted]]

Hello R Help, I am trying solve an MLE convergence problem: I would like to estimate four parameters, p1, p2, mu1, mu2, which relate to the probabilities, P1, P2, P3, of a multinomial (trinomial) distribution. I am using the mle2() function and feeding it a time series dataset composed of four columns: time point, number of successes in category 1, number of successes in category 2, and

I set up a simple BlindUp and BlindDown div that works seamlessy in Firefox but doesn''t work as smoothly in IE. When blinding down, the whole div flashes and then Blinds down. The same thing when blinding up. Is this something that is fixable? _______________________________________________ Rails-spinoffs mailing list Rails-spinoffs-1W37MKcQCpIf0INCOvqR/iCwEArCW2h5@public.gmane.org

Dear all, I need to read in 4 sets of tab delimited data in a loop. The 4 data sets are called "simu1.dat", "simu2.dat" and so on. I know what I need on the righthand side of the read.table expression but I can't the left hand side of it to work (see the line in bold below). Can you kindly help? Many thanks. simu1 <- matrix(0,30,3) simu2 <- matrix(0,30,3) simu3

Dear R-Users, Is it possible to get the covariance matrix from an lme model that did not converge ? I am doing a simulation which entails fitting linear mixed models, using a "for loop". Within each loop, i generate a new data set and analyze it using a mixed model. The loop stops When the "lme function" does not converge for a simulated dataset. I want to

Senior Quantitative Research Analyst Based in Irish Life Investment Managers, Lower Beresford Place, Dublin 1 Irish Life Investment Managers (ILIM), the investment management company within the Irish Life & Permanent Group manages assets of circa EUR30bn and provides fund management services to a large domestic and multi-national client base. We provide investment products to both

Dear R users, I am trying to use tseries' garch function in order to determine the volatility of a return series generated by quantmod. Here is the code that I am using: > library(quantmod) > getSymbols("AAPL") convert daily closing prices into continuous log returns > dret<-dailyReturn(AAPL,type='log') check to see that the autocorrelations decay >

Hi Guys I need help.. I want to making 1 controller for covert html to pdf.. But that html have amchart (swf). When I tried convert, the swf not loadded... Any body can help me Please???? -- Senior Rails Developer Anton Effendi - Wu You Duan -- You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk" group. To post to this group, send

Dear R-programmers, I am trying out certain methods in R, and the statistics require me to calculate n-(sample size) dimensional equations. They are not really very hard to solve - my home-brew implentation of Newton-Raphson in R succeeds most of time with simulated data. (Note that I am assured of a unique solution by theory). Problem comes in with real data, for which I should really implement

Hi. I'm a student in a graduate program at Simon Fraser University in Canada. I am trying to run a simple screening experiment with some simulated data. I simply want to do an ANOVA of an experiemnt with 5 factors (4 have 2 levels, the last has 3 levels) and 48 runs (ie, full factorial). The thing is that I have multiple observations for each level combination (run). So, 1) How do I do

Dear R-listers, I am writing a function that will generate a simple plot for the user . The size of the window (that contains the plot) may not be what the user wants. Supposing the user adjusts the window to a size of her/his choice, I would like to be able to save the new window size and specify the plot to be of this new window size when the user calls my function again. My question is

My app was fine before "rake freeze_edge" but after, I can no longer boot up webrick. I''m on Win Xp MySql backend and this is what I get => Booting WEBrick... ./script/../config/../vendor/rails/activesupport/lib/active_support/core_ext/ker nel/agnostics.rb:7:in ``'': Exec format error - rake tmp:create (Errno::ENOEXEC) from

Railers, I am running an application on Rails 2.3.2. For some reason within our organization, we wish to remain on this version of Rails. For BDD/Testing, we use Cucumber. So I installed Cucumber (0.8.5), Capybara (0.3.9) and Selenium-webdriver (0.0.17). But when I run Cucumber on any feature, I get an error saying *" wrong number of arguments (2 for 1) (ArgumentError)"*, (Check here

Hi all, I wrote a piece of code that generates simulated variables. after variable generation I use them in several analyzes. However, when I use a for to repeat the procedure 1000 times I get an erro message in one of the "for" steps, precisely at this time: gls.temp<- gls(y2 ~ x2,correlation=corExp(form=~coord2[,1]+coord2[,2])) # coord 2 are spatial coordinates and the error

Hello, I have a question regarding competing risk regression using cmprsk package (function crr()). I am using R2.9.1. How can I do to assess the effect of qualitative predictor (gg) with more than two categories (a,b,c) categorie c is the reference category. See above results, gg is considered like a ordered predictor ! Thank you for your help Jan > # simulated data to test > set.seed(10)

Hi all, I'm trying to fit a nonlinear (logistic-like) regression, and I'd like to get some recommendations for which package to use. The expression I want to fit is something like: y ~ A * exp(X * Beta1) / (1 + exp(-(x + X * Beta2 - xmid)/scal)) Basically, it's a logistic function, but I want to be able to modify the saturation amplitude by a few parameters (Beta1) and shift the

Hi Folks! Just started using the gamlss package and I tried a simple code example (see below). Why the negative sigma? John > y <- rt(100, df=1)> m1<-fitDist(y, type="realline")Warning messages:1: In MLE(ll3, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma, : possible convergence problem: optim gave code=1 false convergence (8)2: In MLE(ll4, start = list(eta.mu =