similar to: Using odesolve to produce non-negative solutions

Is there any R package addressing problems of constrained optimization ? I have the following "apparently" simple problem: Given a set V with fixed cardinality: nv Given a set S whose cardinality is a parameter: nHat Let the cardinality of the intersection S.and.V be: nHatv The problem consists of maximizing nHatv/nv subject to a penalty if nHat >

Dear r-helpers, I can run the examples on the mcsamp help page. For example: **************************************** > M1 <- lmer (y1 ~ x + (1|group)) > (M1.sim <- mcsamp (M1)) fit using lmer, 3 chains, each with 1000 iterations (first 500 discarded) n.sims = 1500 iterations saved mean sd 2.5% 25% 50% 75% 97.5% Rhat n.eff beta.(Intercept)

Hello Samba Team, I am new to Linux, but not to networking. However, I am having a bear of a time setting up Samba on a small test network - 3 Machines. I have read all your example setup files for Samba, config, lmhosts, hosts, resolve, etc., and have read back to back "Linux - Samba Server Administration" by Roderick W. Smith (Craig Hunt Linux Library) My test network (10-100

Hi all, 2 days ago, managed to install rhat and asterisk, starting with 2 sip phones, simple config. CLI reports : "playing 'vm-theperson' (language 'no')" when transferred to voicemail after timeout, or "playing 'vm-password' (language 'no')" when dialing into voicemail extension. But, no sound, nothing heard in phones. First time, first

Is anyone using Webmin with CentOS? I've always found it works very well with RHat Linux. I've not used RHEL so I'm not sure. I know on my SuSE server it doesn't seem to work well and suspect there are some location problems. -- <<JAV>>

This is more of a pam question then a samba question, but I thought I'd start here and see if I can get an answer. I've gotten pam_winbind.so working with gdm (on RHAT 8) using the following /etc/pam.d/gdm file. I've put + signs to show the lines I added I added to the stock RHAT 8 gdm pam def.: #%PAM-1.0 + auth sufficient /lib/security/pam_winbind.so + auth

Hi, I have a 4 column file (long/lat/elev/variable) and I tried to fit the values of my variable to the XYZ space using Tps and I keep getting the following message: Warning messages: 1: GCV search gives a minumum at the endpoints of the grid search in: Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, 2: GCV search gives a minumum at the endpoints of the grid search in:

Hi, we have installed ocfs2 on Rhat 4. It is being used as a general clustered file system. There are no oracle binaries or datafiles on the ocfs volume. So no need for RMAN or oarcle agents etc etc. The customer would like to backup the volume using arcserve. However CA are saying that they cannot backup from an ocfs volume. Has anyone out there got any words of wisdom about backing up

Hi All, I have a user who is receiving this error after running the following: library(rjags) library(R2jags) x = c(1.6907, 1.7242, 1.7552, 1.7842, 1.8113, 1.8369, 1.8610, 1.8839) n = c(59, 60, 62, 56, 63, 59, 62, 60) r = c(6, 13, 18, 28, 52, 53, 61, 60) N = 8 data=list("x","n","r","N") inits=function(){ list(alpha.star=0,beta=0) }

Greetings all. I need a transparent way of selectively controlling packet flow between two network devices such that for example I can "see" [e.g. tcpdump] a packet coming into eth0, queue it, and then manually release it out eth1, one packet at a time. See diagram below. (Reason for doing this is to watch the behavior of the device under test to certain packets but I cannot control

Dear all, I tried to run WinBUGS model in R by using package R2WinBUGS, but I failed because of the massage : cannot calculate DIC for model in log file of WinBUGS14. In fact, the model works in WinBUGS and I can get summaries of the model, such as density, stats..., except DIC , because DIC may not be appropriate in this kind of model due to some reasons(according to the manual of WinBUGS14),

Hi All: I had a basic question to ask. I am running R2WinBUGS so that I could automate the running of my model using 1000 simulated datasets. Below is the code I am using. The only problem I am having is the bugs output that comes out shows my parameters as nos with 1 decimal place after. I would want to have the parameters with 5 places after decimal. How would I specify that in my code for

Hi, I'm bit confused about ecdf (read the help files but still not sure about this). I have an analytical expression for the pdf, but want to get the empirical cdf. How do I use this analytical expression with ecdf? If this helps make it concrete, the pdf is: f(u) = \sum_{t = 1}^T 1/n_t \sum_{i = 1}^{n_t} 1/w K((u - u_{it})/w) where K = kernel density estimator, w = weights, and u_{it} =

Dear Users, I would like to simulate an AR(1) (y_t=ct1+y_t-1+e_t) model in R where the innovations are supposed to follow a t-GARCH(1,1) proccess. By t-GARCH I want to mean that: e_t=n_t*sqrt(h_t) and h_t=ct2+a*(e_t)^2+b*h_t-1. where n_t is a random variable with t-Student distribution. If someone could give some guidelines, I can going developing the model. I did it in matlab, but the loops

Hi, it would be nice if the VaR package could get expanded to work with more models. I wrote a function which calculates the VaR for a Garch(1,1) AR(2) model. Since I am fairly new to R I am sure it can be written much neater. If someone feels like having a go, here is what I have written library(fSeries) library(tseries) VaR.Garch <- function(data,stockId=1,p=0.01,dt=1) { r <-

I am just starting to learn Sweave (really neat tool). I am pretty early in the learning curve (I had to think a moment ago whether a # or % was the appropriate comment character). I have successfully incorporated simple graphics and outputs, but am having trouble getting a latex (xtable) table to function properly. Latex is seemingly treating the xtable code as input or verbatim text. That

Have any of you folks actually managed to get pam_mount working? A quick google shows a ton of messages saying "you can use pam_mount" to automatically mount a user's home directory on log in, but no messages saying "I use pam_mount" etc -- I'm suspecting it might not actually work. Or at least I might not be smart enough to make it work. First -- pam_mount 0.5.11

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints

starcraft is a game ive enjoyed for a long time and last time i played it was on windows. now im trying it on ubuntu feisty, with wine installed. starcraft installed amazingly fast, but when i started it up, its as if the computer was giving no visual output, and the monitor itself gave me an error message. basicly i could still hear the starcraft music in the background, but i couldnt see a

Hi - I've been using the option main=bquote(paste(mu==.(mu),", ",lambda==.(lambda),", ",truncation==.(truncation),", ",N[T]==.(n))) to produce a title when using the "plot" command - a title which includes variable names (two Greek) along with their values. The above option, however, does not work within the "boxplot" command. Any