similar to: custom loss function + nonlinear models

Hi, I want to impute the missing values in my data set multiple times, and then combine the results (like multiple imputation, but manually) to get a mean of the parameter(s) from the multiple imputations. Does anyone know how to do this? I have the following script: y1 <- rnorm(20,0,3) y2 <- rnorm(20,3,3) y3 <- rnorm(20,3,3) y4 <- rnorm(20,6,3) y <- c(y1,y2,y3,y4) x1 <-

Hi all, I'm looking for the algorithm to calculate maximum intrinsic and parameter effects curvature of Bates & Watts (1980). I have Bates & Watts (1980) original article, Bates et al (1983) article, Seber & Wild (1989) and Ratkowsky (1983) books. Those sources show steps and algorithm to get this measures but I don't know translate C code to R language and I've no success

Hello, I would like to run a script in which a loop is included. Since I'm new to R, I cannot manage the following problem. I really hope someone could help me out. Data in the variable Y should be removed from the simulated data set with probability 0.50 if the variable X has a value below zero, and with probability 0.10 if X has a value above zero (see script). However, the total number of

Hi, I wrote a script in order to simulate data, which I will use for evaluating missing data and imputation. However, I'm having trouble with the last part of my script, in which a dataframe is constructed without missing values. This is my script: y1 <- rnorm(10,0,3) y2 <- rnorm(10,3,3) y3 <- rnorm(10,3,3) y4 <- rnorm(10,6,3) y <- c(y1,y2,y3,y4) a1 <-rep(1,20) a2

I try to compile samba-2.07 in a separate directory and it fails with : Linking bin/smbd ld: 0711-934 SEVERE ERROR: Cannot create the output file: bin/smbd ld:open() A file or directory in the path name does not exist. collect2: ld returned 12 exit status gnumake: *** [bin/smbd] Error 1 because the directory `bin' does not exist in my build directory; the

Dear All,i applied "NbClust", to my data to find optimum number of clusters, and got following resultsNow, i don't know how to read these results. more precisely, i would like to know, how to see which method is more precise for my data considering these index values.your help is needed...thanks in advance Eliza Botto > dput(Eliza)structure(list(All.index = structure(c(2, 3, 4, 5,

Dear All, I am trying to specify the following fixed effects model for lme: y ~ constant1 - beta1*(x - beta2) where y is the response, x is the independent variable, and the operators above are real arithmetic operations of addition, subtraction, and multiplication. I realize that this model is just a reparameterization of y=beta0+beta1*x, but I am using this parameterization because I am

Buenos dias, desde hace algunos dias estoy realizando un trabajo,mi computadora es una DELL, windows 10 64 bits, 8G de RAM y disco de estado solido, estoy procesando 29000 filas y 23 columnas, mi codigo es este: nb <- NbClust(datos.scaled, distance = "euclidean", min.nc = 2, #max.nc = 10, method = "complete", index ="all") . y mi error es este: Error:

I am the maintainer of one of the CRAN mirrors in the U.S., cran.us.r-project.org. It appears that this site has been inaccessible since Friday. This has me puzzled because I can connect to the same machine through ssh and I can see that the http daemon (apache) is running. In fact I can make an http connection from this machine to this machine and from other machines on any of our subnets to

[This email is either empty or too large to be displayed at this time]

Ronaldo Your mathematical question is of interest to me too. But I am particularly interested in a translation into American English of the pun attached to you email. Could you try to explain O Flamengo e o Bangu terminaram o jogo em 0 x 0. Quem fez o gol? A Volkswagen My rough translation into English is: Two soccer teams ended play 0 to 0. Who got the goal? A Volkswagen. I need all the

Hi, I have this non-linear function: y=115.251 - 118.69 * exp(-0.123517*x) I try to discovery the x value for the assyntoptic value of y and the x value where the behavior of curve change, the inflection point. How to make this? Is poss?ble to make this on R? Thanks Ronaldo -- O Flamengo e o Bangu terminaram o jogo em 0 x 0. Quem fez o gol? A Volkswagen -- | // | \\

Hi all, The Rails book talks about running Rails under Apache, but is there a (relatively) easy way to deploy it to either Tomcat or JBoss? Is the CGI servlet the only option? Thanks, Ken -- Kenneth A. Kousen, Ph.D. President Kousen IT, Inc. <http://www.kousenit.com> http://www.kousenit.com <mailto:ken.kousen@kousenit.com> ken.kousen@kousenit.com

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

Many thanks to those took time replied to my question. They were very helpful and I solved my problem by reparameterization. With the help of optim() the fitting procedure is very robust and insensitive to initial starting value. Once again, many thanks. Kevin -------- Original post --------- >Dear R and S users, > >I've been working on fitting finite mixture of negative

Any and all attempts to connect to our samba server, packaged in the redhat RPM found on samba.org, have failed. The allow statement is set to localhost, x.y.z. where x.y.z is our class C(ended with a period, as per the online help in SWAT.) Local connections via smbclient work beautifully, but connections from any other machine fails. I've followed the troubleshooting procedure to the

Dear friends: I am a new R user, and not very good at statistics and maths I find it difficult to find the minimum of this item: ((p*(b-1)-1)*(p+1)^(b-1)+1)/p^2 It seems that the method optim can find the minimum. I tried several times , but R constanly tells that " failure to find b or p" Would you please give some suggestions? Thanks in advance. Best

Dear R-ers, I am having trouble understanding why I am getting an error using glmmPQL (library MASS). I am getting the following error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 The long story: I have data from an experiment on pairwise comparisons between 3 treatments (a, b, c). So a typical run of an experiment

Hi all, If presented with a singular design matrix, lm drops columns to make the design matrix non-singular. What algorithm is used to select which (and how many) column(s) to drop? Particularly, given a factor, how does lm choose levels of the factor to discard? Thanks for the help. Best, Anirban [[alternative HTML version deleted]]

When I SSH using protocol 1 from a Debian box running OpenSSH 2.9p2-4 to a sparc.sunos5 box running vanilla OpenSSH 2.9.1, after a little while (of inactivity?) I get the following message on the client terminal: Disconnecting: protocol error: rcvd type 98 Looking further, this message is actually caused by the SSH daemon. However, I'm at a loss to determine why sshd is doing this. I attach