similar to: R-help

Hi I need help figure out how to fix my code. When I call into R >optimize(llik,init.params=F) I get this error message ####Error in optimize(llik, init.params = F) : element 1 is empty; the part of the args list of 'min' being evaluated was: (interval)#### My data and my code looks like below. R_j R_m 0.002 0.026567296 0.01 0.003194435 . . . . . . . . 0.0006

Hi May you help me correct my loop function. I want optim to estimates al_j; au_j; sigma_j; b_j by looking at 0 to 20, 21 to 40, 41 to 60 data points. The final result should have 4 columns of each of the estimates AND 4 rows of each of 0 to 20, 21 to 40, 41 to 60. ###MY code is n=20 runs=4 out=matrix(0,nrow=runs) llik = function(x) { al_j=x[1]; au_j=x[2]; sigma_j=x[3]; b_j=x[4]

Hello, the more general thing I'd like to learn here is how to compute Function of Data on the basis of grouping determiend by n variables. In terms of the reason why I am interested in this, I need to compute the average of my data based on the value of the month and day across years. I have come up withy the code below which, as far as I can see, does what I need but getting either a more

Hi I have developed the code below. I am worried that the parameters I want to be estimated are "not being found" when I ran my code. Is there a way I can code them so that R recognize that they should be estimated. This is the error I am getting. > out1=optim(llik,par=start.par) Error in pnorm(au_j, mean = b_j * R_m, sd = sigma_j) : object 'au_j' not found #Yet

Dear Sir/Madam, I am using R version 1.8.1. I am doing following tast: First generate 100 Gaussion(3,1) numbers, then write the likelihood function to estimate the parameters of Gaussian distribution by direct maximizing the likelihood function. My likelihood function is: >fn<-function(x) >(-50*log((sd(x))^2))-50*log(sqrt(2*pi))-(1/2*((mean(x))^2))*(sum((x-(mean(x))^2)) After I

Everyone, I have installed Centos 7.0 on my homework machine in order to take a test drive with it, and am low on the learning curve with it at this point. I have a small Gateway SX2855-UB12P. I have a critical hurdle in that when I try a reboot or when I do a 'shutdown now -r' command the system will start a reboot process but hangs right after the os choices are presented. My

On 11/23/2014 12:02 PM, Edward M wrote: > > On 11/23/2014 10:58 AM, Gregory P. Ennis wrote: >> Everyone, >> >> I have installed Centos 7.0 on my homework machine in order to take a >> test drive with it, and am low on the learning curve with it at this >> point. I have a small Gateway SX2855-UB12P. >> >> I have a critical hurdle in that when I try

I''m a ruby newbie and have what is probably a relatively simple problem I''m trying to solve. Say that I created a hash: TASKS = { "homework" => "Do Your Homework", "chores" => "Your Have Chores", "exercise" => "Don''t Forget to Exercise" } and then I had an Assignment model with: def

Dear Sir, I am using R to estimate two parameters in Normal distribution. I generated 100 normal distributed numbers, on which to estimate the parameter. The syntax is: >fn<-function(x)-50*log((y)^2)+50*log(2*pi)-(1/2*(z^2))*(sum((x-y)^2)) >out<-nlm(fn, x, hessian=TRUE) but it does not work. Could you please help me to compose the syntax for the purpose that find maximum

*** Travel and Accommodation Support *** Funds from the U.S. National Science Foundation may be available to provide partial support for travel and accommodation for some graduate students and junior faculty at U.S. post-secondary institutions to attend DSC 2009 and useR 2009. If you wish to apply for this support, please send an application to luke at stat.uiowa.edu. Your application should

*** Travel and Accommodation Support *** Funds from the U.S. National Science Foundation may be available to provide partial support for travel and accommodation for some graduate students and junior faculty at U.S. post-secondary institutions to attend DSC 2009 and useR 2009. If you wish to apply for this support, please send an application to luke at stat.uiowa.edu. Your application should

R version: 2.7.0 Running on: WinXP I am trying to model damage from fire losses (given that the loss occurred). Since I have the individual insured amounts, rather than sampling dollar damage from a continuous distribution ranging from 0 to infinity, I want to sample from a percent damage distribution from 0-100%. One obvious solution is to use runif(n, min=0, max=1), but this does not seem to be

Evening, I'm interested in finding out the current status of EVE Online under wine and creating a central point for all work to support the game (maybe in wiki form). Could any one with information please reply with locations of patches and what is happening at the moment. Also I'm interested as to how the progression of the appdb is going and if any one else likes the idea of setting

Hi all, I'd like to fit a multivariate regression with the variance of the error term porportional to the predictors, like the WLS in the univariate case. y_1~x_1+x_2 y_2~x_1+x_2 var(y_1)=x_1*sigma_1^2 var(y_2)=x_2*sigma_2^2 cov(y_1,y_2)=sqrt(x_1*x_2)*sigma_12^2 How can I specify this in R? Is there a corresponding function to the univariate specification lm(y~x,weights=x)??

Hi R users: I will like to know if somebody works on insurance statistics (actuarial problems) and had use TRICAST, and can tell me if with all the R tools it can be build a solution like TRICAST or similar. In a word: Do you think that R has all the statistical tools (I mean modeling tools) to make a job similar to TRICAST? Does TRICAST has modeling tools that are not implemented on R yet?

Hello R Power Users, There is a rather large introduction you may skip to go directly to my question about S4 classes. I am working on some toy code to help me get in through S4 classes. Doing this homework, I have some questions about S4 classes. I have carefully read available help and some additional material such as Robert Gentleman's slides and special attention to "Programming

Has anybody written an actuarial (life) survival procedure, this does not appear to be an option in the survival package? This approach is common in orthopaedic surgery to demonstrate the survival of prostheses. I need to apply the "modified" lower conf.int because of the censoring over time. I want a life table which I can then easily plot. Many thanks

I have an NI2 PRI connection to a 5ESS switch. I am explicitly told to use NI2 as the protocol, not 5ESS. I am running Asterisk 1.2.4. When making calls to the PSTN, Asterisk sends the calling name as part of the display information element in the SETUP message. I need to be able to send the calling name as part of the facility information element in the SETUP message. I've seen that work

I'll throw this thought out here so it doesn't get lost: when it came time for me to build a Windows release, I always used a quarantined Windows box that had the minimum stuff installed and had never been on a network, to avoid malware getting into the binaries. The last thing I ever wanted to hear was some Windows user blaming FLAC because a bad build infected him.? It was bad enough

Luke, Thank you for sharing the benchmark results. The improvement is very substantial, I am looking forward to the release of the byte compiler! The arithmetic shows that x[i]<- is still the bottleneck. I suspect that this is due to a very involved dispatching/search for the appropriate function on the C level. There might be significant gain if loops somehow cached the result of the initial