similar to: London puppet masters

Hello everyone, I am new to R. I have data in the form of excel pivot table format and I want to cast it into a format which can make it compatible with computation. Since I already have the package in pivot format; I avoid melt function and use the cast directly. I inspect the dataread <- read.csv(".....", header=TRUE) Data in the format Id Region Country Industry

I'm having trouble figuring out how to format Date variables when used as axis labels in graphs. The particular case here is an attempt to re-create Nightingale's coxcomb graph with ggplot2, where I'd like the months to be labeled as "Mar 1885", "Apr 1885", using a date format of "%b %Y" applied to label the dates, or really anything other than

Hi everyone, I have already started documented and building this privately. However I refuse to believe that there is not a module for setting up a FreeIPA client server on RHEL to authenticate against an existing IPA authentication server. Has anyone actually already done this? If not I will probable post a link to this email with the puppet module I have built. Regards John -- You

Dear R-er, For a non-linear optimisation, I used optim() with BFGS method but it stopped regularly before to reach a true mimimum. It was not a problem with limit of iterations, just a local minimum. I was able sometimes to reach better minimum using several rounds of optim(). Then I moved to optimx() to do the different optim rounds automatically using "Nelder-Mead" and

Hi There, There are 300 boy students and 100 girl students in a class. One interesting question is whether boy is smarter than girl or not. first given the exam with a difficulty level 1, the number of the student who got A is below 31 for boy, 10 for girl. Then we increase the difficulty level of the exam to level 2, the number of the student who got A is below 32 for boy, 10 for girl. We

Hi, I have written an R translation of C.T. Kelley's Matlab version of the Nelder-Mead algorithm. This algorithm is discussed in detail in his book "Iterative methods for optimization" (SIAM 1999, Chapter 8). I have tested this relatively extensively on a number of smooth and non-smooth problems. It performs well, in general, and it almost always outperforms optim's

How does one implement a likelihood-ratio test, to test whether the variances of the random effects differ between two groups of subjects? Suppose your data consist of repeated measures on subjects belonging to two groups, say boys and girls, and you are fitting a linear mixed-effects model for the response as a function of time. The within-subject errors (residuals) have the same variance in

Hi r-users, I would like to solve system of nonlinear equation using BBsolve function and below is my code.  I have 4 parameters and I have 4 eqns. mgf_gammasum <- function(p) { t  <- rep(NA, length(p)) mn <- 142.36 vr <- 9335.69 sk <- 0.8139635 kur <- 3.252591 rh  <- 0.896 # cumulants k1 <- p[1]*(p[2]+p[3]) k2 <- p[1]*(2*p[2]*p[3]*p[4] +p[2]^2+p[3]^2) k3 <-

We are having a problem getting a NetApp filer to re-join a samba domain after a move to a new network. The filer worked fine with samba before the move. Apologies in advance for the long missive. I've tried the following: - re-running the CIFS setup program on the filer - removing the problem filer's samba account, replacing it, and re-running the setup program on the filer

Hi everyone, There seems to be a good level of interest in this; so we should settle on a time and location. It has been pointed out that previous socials were held at http://bbrovers.com/ - we seem to have a preference for the availability of alcoholic beverages, and otherwise no strong opinions have been expressed. Logistically, we should pick someplace likely to have flexibility regarding

Dear friends. Please see the program below and answer if it does simulate a population of 1.000.000 families, each with at max 20000 children (typical in Denmark, you know), constructed such that each family stops having children when more boys than girls are present ? Equal numbers of boys and girls are got in the population, according to the simulation, is that obvious ? ND <- NP <-

I have a Master's student working on a project which involves estimating parameters of a certain model via maximum likelihood, with the maximization being done via optim(). A phenomenon has occurred which I am at a loss to explain. If we use certain pairs of starting values for optim(), it simply returns those values as the ``optimal'' values, although they are definitely not

Let''s say I have three hypothetical MySQL tables: ? people, with columns id, gender, and source_id belongs_to :boys and :girls ? boys, with columns id and name has_many :people ? girls, with columns id and name has_many :people The gender column in people specifies which of the two source tables the source_id refers to. For example, if we have values: 1, boy, 1 in people,

..::: Drunk Young :::... http://210.101.95.254/drunk-young ^^^^ ^^^^ Click here to Join Drunk Young will take you to the new level of excitement. It is a guarantee of Your pleasure, because we choose models very seriously & carefully, depending on you offers and suggestions. 'Drunk Young' is regularly updates to bring you more and more sweet nymphets 100%

Watch Movies Online For Free Full Movie Downloads,stream movies online free,trailer, New Moon ,2012 ,Taylor Lautner,Twilight ,Precious ,Transformers: Revenge of the Fallen,Harry Potter and the Half-Blood Prince,Up,The Hangover,Star Trek,The Twilight Saga: New Moon,Monsters Vs. Aliens,Ice Age: Dawn of the Dinosaurs,X-Men Origins: Wolverine,Night at the Museum: Battle of the Smithsonian,The

Dear list,   I am looking for an R-only implementation of a Nelder-Mead process that can find local maxima of a spatially distributed variable, e.g. height, on a spatial grid, and outputs the coordinates of the new point during each evaluation. I have found two previous threads about this topic, and was wondering if something similar has been implemented since those messages were posted.   Thank

All, In looking at `optim', it doesn't appear that it is possible to impose nonlinear constraints on Nelder- Mead. I am sufficiently motivated to try to code something in C from scratch and try to call it from R.... Does anyone have some good references to barrier and/or penalization methods for Nelder-Mead? I would ideally like some papers with pseudocode for method(s) that are in

Hi! I'm trying to figure out how to use optim... I get some really strange results, so I guess I got something wrong. I defined the following function which should be minimized: errorFunction <- function(localShifts,globalShift,fileName,experimentalPI,lambda) { lambda <- 1/sqrt(147) # error <- abs(errHuber(localShifts,globalShift, #

Hi everyone, since I subscribed to this mailing list (about two days ago) I`m receiving the attached message about 5-10 times a day. Has anybody had the same problem? And how can the problem be solved?? Best regads, Sven P.S.: In case anybody wonders, the email address used to subscribe to the mailing list has never been used before. Begin forwarded message: > From: gee_two@ru.ru >

Dears useRs, I am using the mle function but this gives me the follow erros that I don't understand. Perhaps there is someone that can help me. thank you for you atention. Bernardo. > erizo <- read.csv("Datos_Stokes_1.csv", header = TRUE) > head(erizo) EDAD TALLA 1 0 7.7 2 1 14.5 3 1 16.9 4 1 13.2 5 1 24.4 6 1 22.5 > TAN <-