similar to: mixed ordinal logistic regression

Dear List, I'm trying to solve this simple optimization problem in R. The parameters are the exponents to the matrix mm. The constraints specify that each row of the parameter matrix should sum to 1 and their product to 0. I don't understand why the constraints are not satisfied at the solution. I must be misinterpreting how to specify the constrains somehow. library(alabama) ff <-

Hello everyone, I want to calculate type I error rate for modified F statistic for one way robust anova. I need to find the j group trimmed mean and winsorized sum of squared deviations. Here I attached my code for j=2 to make it simple. Originally I have j=4. Hope you can help. I need to run it for 1000 times My problem is: i) the value of F-test obtain from my simulation below is in negative

Dear List, I'm new in R. I'm trying to solve a simple constrained optimization problem. Essentially, let's say I have a matrix as in the object 'mm' inside the function below. My objective function should have a matrix of parameters, one parameter for each element 'mm' (4 in this case). The problem is to select the value of the parameters to maximize the function

1. The matrices are only 3x2. In your loop you'd be assigning rows 4 and 5 in the last two iterations. Are you sure that's what you want? 2. The reason it ``didn't work'' is because assign() takes the first argument as the name of the object to create, literally, instead of evaluating it. 3. You probably ought to be working with a list. E.g.: m <- replicate(5,

Dear All, I have a dataframe made up of individual beetles consisting of individual number, family number, mother's family number, father's family number, and sex of the beetle. I would like to pair up the individuals for breeding. I would, however, like to avoid breeding beetles of the same sex (obviously), the same family, and with the same mother's family or father's family,

Hello all i have postfix running with dovecot-sasl and mysql as a backend. It all runs good. I run into trouble as where outlook 2003 fails to authenticate when sending e-mail. I have thunderbird outlook2007 and 2003 clients. The tunderbird and 2007 clients are working OK, the outlook2003 client get the relay access denied message. In the postfix log i see it is not initiating sasl they all

I am runnig posfix on Centos 4.4 as a Mailgateway. It only accepts mails for domains and then forwards mails to Lotus domino Server. All clients sends outgoing mails to that Lotus domino Server. Then , That Lotus Domino Server sends mails to Postfix mailgateway. This postfix mailgateway sends mails to all the destinations. But, This Postfix mailgateway has about 150 messages in the mailq. Some

Dear all, I've sent this question 2 days ago and got response from Sarah. Thanks for that. But unfortunately, it did not really solve our problem. The main issue is that we want to use our own (manipulated) covariance matrix in the calculation of the mahalanobis distance. Does anyone know how to vectorize the below code instead of using a loop (which slows it down)? I'd really appreciate

>> Can I change my vote to http://samugliestdog.com/Sam162edited.jpg yikes! how''d that get past my firewall? +1 On 6/30/06, mongrel-users-request at rubyforge.org <mongrel-users-request at rubyforge.org> wrote: > Send Mongrel-users mailing list submissions to > mongrel-users at rubyforge.org > > To subscribe or unsubscribe via the World Wide Web, visit

Is there a way to listen on port 25 for repeated dictionary attacks to harvest email address and blacklist that Ip with shorewall? Thanks, Mike

Hi, As it has been brought to my attention (by several parties), my current MB does not support VT-d. Since I don''t feel like purchasing a brand new system, I have searched for replacement MB-s that would fit into my existing uATX case, support my C2D CPU, and support at least 4x2GB ram. I have narrowed it down to the following four models: SUPERMICRO MBD-C2SBM-Q-O (Q35 chipset)

Hi everyone! I have some doubts about mixed effect models and I hope someone could help me. I´m trying to analyze a dataset coming from samples of dung beetles in the same forest fragments along 3 consecutive years (1994, 1995 and 1996) and 14 years after (2010). I sampled dung beetles in 18 different fragments with different sizes and different degrees of isolation. My aim is to determine

Dear R-Helpers, Can anyone please help me to interpret warning messages from zeroinfl (package pscl) while fitting a zero inflated negative binomial model? The console reports convergence and the parameters seam reasonable, but these <<Warning messages: 1: algorithm did not converge in: glm.fit(X, Y, family = poisson()) 2: fitted rates numerically 0 occurred in: glm.fit(X, Y, family =

Hello from a presently sunny Helsinki! I've bee trying to repeat an analysis I did about 18 months ago (the reviewers of the paper want something adding to it). I'm using R1.5.0, but I couldn't see anything in the list of changes between this and 1.5.1 to suggest it would act any differently. The data is observational, on the changes in the population status of carabid beetles, and

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to

I'm hoping to take advantage of authentication caching via dovecot-auth. I see that posfix can communicate with Dovecot-Auth, via SASL, but from what I can see, postfix only does this for authentication checks. Not 100% sure how postfix handles virtual_transport (I assume it needs to do a lookup via virtual_mailbox_maps), but I can't see anyway to get virtual_mailbox_maps to use

Hi I'm having a problem with model.frame, encapsulated in this example: y1 <- matrix(c(3,1,0,1,0,1,1,0,0,0,1,0,0,0,1,1,0,1,1,1), nrow = 5, byrow = TRUE) y1 <- as.data.frame(y1) rownames(y1) <- paste("site", 1:5, sep = "") colnames(y1) <- paste("spp", 1:4, sep = "") y1 model.frame(~ y1) Error in model.frame(formula, rownames,

Hello, I am using function lmodel2 to calculate RMA regression between life-history traits in ladybirds beetles. It works well but I am not able to plot an RMA regression line on the scatterplot of my data. I am of course unable to plot the confidence intervals. For ordinary least square regression abline(lm(x~y))works well but for RMA regression, abline(lmodel2(x~y))does not do it. If somebody

We have a bunch of lists hosted on a CentOS-4.9 server running sendmail and mailman that we are moving to a CentOS-6 vm running postfix and mailman. I am testing the setup and I am running into a problem no doubt caused by my own ignorance wrt postfix. Based on groking the web I discovered that the recommended way to configure postfix to work with mailman is a script called

Hi, I have a couple of samples that produce interesting artifacts when encoded with the CVS snapshot of 2001-02-13. Both are about a meg. ftp://slumber.dhs.org/tmp/4.wav.bz2 When encoded with oggenc -b 128, there is a sort of stereo separation in the drums, while in the original they are 'solidly' positioned in the stereo field. This also occurs with -b 160 and is barely audible with