search for: yob

...R ? Thanks Arnaud ###################################################################### ## function 1 x1 <- 1:100 y1 <- 5.468 * x + 3 # + rnorm(100,0, 10) dfxy <- cbind(x1,y1) # Objective function optfunc <- function(x, dfxy){ a <- x[1] b <- x[2] xtest <- dfxy[,1] yobs <- dfxy[,2] ysim <- a*xtest + b sum((ysim - yobs)^2) } out<- optim(par=c(0.2,5), fn=function(x){optfunc(x, dfxy)}, method = "Nelder-Mead", hessian = F) ## function 2 x2 <- seq(0.01, 0.1, length=100) y2 <- exp(30*x2) dfxy2 <- cbind(x2,y2) # objective function op...

...> E[i]<- 1/( i+ sum(R[(m+1-i):m])) > i<-i+1 > } > > V<-W^E > > i<-1 > while (i<=m) { > > U[i]<- (1- prod( V[(m+1-i):m])) > i<-i+1 > } > > x<-( (1-U)^(-1/theta)-1 )/beta # simulation ends > > Yobs<-x > > em.lomax <- function(Y){ > r <- length(Yobs) > > ##initial value > > theta.t=theta > beta.t=beta > > # Define log-likelihood function > ll <- function(y, k, c){ > n*log(c*k)-(k+1)*( sum(log(1+c*Y))+sum( R*log(1+c*Y)+1/k ) ) > } > &g...

Hello, I am using logistic discriminant analysis to check whether a known classification Yobs can be predicted by few continuous variables X. What I do is to predict class probabilities with multinom() in nnet(), obtaining a predicted classification Ypred and then compute the percentage P(obs) of objects classified the same in Yobs and Ypred. My problem now is to figure out whether P(o...

...HE Submission from: (NULL) (84.9.38.207) I'm using the FEAR library to perform Data Envelopment analysis with a 36,000 obs dataset. The function dea is working fine but then when I try to use the boot.sw98 to come up with some sensitivity analysis I get the following error: boot.sw98(XOBS=x,YOBS=y,NREP=100,DHAT=dhat,RTS=3,ORIENTATION=2) Error in 0:L : NA/NaN argument In addition: Warning message: longer object length is not a multiple of shorter object length in: kappam * Gcounts What does this mean? Is this because I have too many observations? What can I do to fix it? Please l...

Dear R users, I am using the MICE package. Specifically, at some point in my code I have imp2=mice(PoptotalMICE,imputationMethod="logreg2") And R returns... iter imp variable 1 1 MICEYError in logitreg(xobs, yobs, intercept=F) : couldn't find function "ms" I have been looking for this ms function on the web, hoping it was just a matter of downloading a package. I didn't find anything. I have tried getAnywhere without success. What can it be? Where is the function ms? Thanks...

...m. for (i in 2:N){ xp[[i]]=C%*%xf[[i-1]] Pp[[i]]=C%*%Pf[[i-1]]%*%t(C)+Q siginv=A[[i]]%*%Pp[[i]]%*%t(A[[i]])+R sig[[i]]=(t(siginv)+siginv)/2 # make sure sig is symmetric siginv=solve(sig[[i]]) # now siginv is sig[[i]]^{-1} K=Pp[[i]]%*%t(A[[i]])%*%siginv innov[[i]]=as.matrix(yobs[i,])-A[[i]]%*%xp[[i]] xf[[i]]=xp[[i]]+K%*%innov[[i]] Pf[[i]]=Pp[[i]]-K%*%A[[i]]%*%Pp[[i]] like= like + log(det(sig[[i]])) + t(innov[[i]])%*%siginv%*%innov[[i]] } like=0.5*like list(xp=xp,Pp=Pp,xf=xf,Pf=Pf,like=like,innov=innov,sig=sig,Kn=K) } I tried to fit my problem and observe that I...

...I get mean of perr1 of 0.2336 and SD of 0.000139. When using a home-made cross validation, code below, I get mean of perr2 of 0.2338 and SD of 0.02184. The means are similar but SD are different. Questions are: (1) how the $delta is computed in the cv.glm? In the home-made version, I simply use ((Yobs - Ypred)^2)/n. The equation might be correct because the mean is similar. (2) in the cv.glm, I have the impression the system is using glm0.dmi that was generated using all the data points whereas in my homemade version I only use the "test" database. I am confused if the cv.glm generat...

...running fine, so i decided to run something more fun, starcraft. The installation program works, even with sound and all, exept when i decide to run the previews. The previews run rine, but whithout sound. When i return to the installation program the sound is gone. But anyway, it does it's yob, i.e. it installes starcraft. When i start the game i get the following error: WineDbg starting on pid f Loaded debug information from ELF 'wine' ((nil)) No debug information in 32bit DLL 'C:\Program Files\Starcraft\StarCraft.exe' (0x400000) No debug information in 32bit DLL '...

I. most of ssh manual and all sshd manual present server and client as one machine, called host. All files mentioned are placed on one machine. This is incorrect, and makes the explanation unclear. For example, man sshd SSH_KNOWN_HOSTS FILE FORMAT suggests to copy keys from /etc/ssh/ssh_host_key.pub into /etc/ssh/ssh_known_hosts, as if those files are on the same machine. II. a general

I have data that looks like this: > df1 group id 1 red A 2 red B 3 red C 4 blue D 5 blue E 6 blue F I want a list of the groups containing vectors with the ids. I am avoiding subset(), as it is only recommended for interactive use. Here's what I have so far: df1 <- data.frame(group=c("red", "red", "red", "blue",

...221UI1uPgoAy796QQinElDqmn0ogWTioECCovJ06wKBTWxh1XHaO66Iu/FUIGEJgdBFGl5Hl dX+ROrPVbuh/zKWitka4WVp72VkvzYrCbafaso9UgHIGatQXR4Jb6vZgKR1oyzqSt57n9V8r6k+W 8EMK2PRHB3Bm3RhZ1ubJD52lrUt5X+rZv7CczO3Lcil6H8kPobaxSRNlRZD2PvpUTqGLlchC1y7C WsubK0omV40MJcMSFQAchZmjL8Un/TH/ItDUCEg0eVTsEKMtTtRV3SrzsOYRPorMv+UchAyP+yOb Sku2ND36ZAwi9STAsJac72Zt7kbZVHJ6uRFvHragtco34ZfJyo9Wgz6pbQjMWDvwpraPPtfRINW5 qh6qpQ8tqNNTBheaxvqEkHk0r2ZvLj+ksbpQQ7dBQW1QJKS7nkcrkOfkKg9ojc5t8vxAnJZcW323 OI+aOHS4oCDg2MP8c9ZtESc7z/rmMdgXdYhTi/LIqlLDXfHBm67Y63OJnr9igTeVixGZJghD5VU8 kx5YY6kvmdmAN3dErw/0AaXqizk6+/etPtKo6zkUaxofjoiQn3VvZEi1XoVhJWc+h5Bl...