similar to: Weighted Optimization

Hi All, I've put together some script which gives me the parameters of a modified NIG distribution. Now I'd like to include a weighted vector within the maximization function. The code below gives me parameter estimates which are implicitly equally weighted. nig.par.fit <- try(optim(vega, negllh, hessian = se, pdf = density.nig, tmp.data = data, transf =

Hi, #make example data dat <- data.frame(matrix(rnorm(15),ncol=5)) colnames(dat) <- c("ab","cd","ef","gh","ij") If I want to get a subset of the data for the middle 3 columns, and I know the names of the start column and the end column, I can do this: mysub <- subset(dat,select=c(cd:gh)) If I wanted to do this just on the column names,

Hi. I have a CentOS 6.2 virtual machine in a vmware ESXi 4.0 host 4G Ram 4 virtaul cpus and with about 4 TB of disk space formatted with XFS . I'm seeing a lot of : Mar 21 10:49:42 nzhmlfpr05 kernel: XFS (dm-4): xlog_space_left: head behind tail Mar 21 10:49:42 nzhmlfpr05 kernel: tail_cycle = 129, tail_bytes = 20163072 Mar 21 10:49:42 nzhmlfpr05 kernel: GH cycle = 129, GH bytes =

below is a code to compute hessian matrix , which i need to generate 29 number of different matrices for example first element in x1 and x2 is use to generate let say matrix (M1) and second element in x1 and x2 give matrix (M2) upto  matrix (M29) corresponding to the total number of observations and b1 and b2 are constant.  can some one guide me or help to implement this please. I did not

Sean, I add the empty file .nojekyll to https://github.com/Jonathan2251/lbd/tree/gh-pages and do "make gh-pages" again. It's the same. I check the uvbook has the gh-pages as mine as follows, https://github.com/nikhilm/uvbook/tree/gh-pages https://github.com/Jonathan2251/lbd/tree/gh-pages In addition to gh-pages, uvbook has the the other web site on

Dear R users, I'm trying to fit a simple random intercept model with a fixed intercept. Suppose I want to assign a weight w_i to the i-th contribute to the log-likelihood, i.e. w_i * logLik_i where logLik_i is the log-likelihood for the i-th subject. I want to maximize the likelihood for N subjects Sum_i {w_i * logLik_i} Here is a simple example to reproduce

Dear list, I have following data set, which I want to plot the "Scale" variable on the x-axis and "Mean"´on the y-axis for each Ageclass and for each sex. The Mean value of each Ageclass for each sex would be connected by a line. Totally, there should be 6 lines, from which three present the Mean values of each Ageclass for respective sex. Are there any easy ways to do

Dear all, Over the last few weeks with the help of GH folks I've been exploring the options of Bugzilla migration. I believe finally we came to the viable solution which is detailed below. It turned out that GitHub has an internal project rehydration tool that could be used to populate the empty repo contents from the simple serialized format. There is a big advantage of this approach as

Dear list, Is there a way to reorder the xaxis using lattice. Using the following data, the x-axis is ordered as BP GH MH PF RE RP SF VT but I would like the x-axis to be ordered as PF RP BP GH VT SF RE MH. Kön Skalor Tillfälle Medelvärde 1 Kvinnor BP 1-inskrivning 36.45283 2 Kvinnor GH 1-inskrivning 38.62255 3 Kvinnor MH 1-inskrivning 62.88889 4

Sean, It is not work even though I change the PDF link as you indicated. I think I can just use directory/link https://github.com/Jonathan2251/lbd/tree/gh-pages instead http://jonathan2251.github.com/lbd/WriteAnLLVMBackendTutorialForCpu0.pdf as a work around solution. User can click the link and click left top button "ZIP" and get the epub and pdf. Jonathan --- 12/12/4 (二)，Sean Silva

You may need a ".nojekyll" file. uvbook has one < https://github.com/nikhilm/uvbook/blob/gh-pages/.nojekyll>. Also see < https://help.github.com/articles/files-that-start-with-an-underscore-are-missing >. Let me know if that doesn't work. -- Sean Silva On Mon, Dec 3, 2012 at 9:15 PM, gamma_chen <gamma_chen at yahoo.com.tw> wrote: > Sean, > > I can

I am currently attempting to split a long list of strings (let's call it "string.list") that is of the format: "1234567.z3.abcdef-gh.12" I have gotten it to: "1234567" "z3" "abcdef-gh" "12" by use of the strsplit function. This leaves me with each element of "string.list" having a split string of the above format. What

R Experts: I am conducting a meta-analysis where the effect measures to be pooled are simple proportions. For example, consider this data from Fleiss/Levin/Paik's Statistical methods for rates and proportions (2003, p189) on smokers: Study N Event P(Event) 1 86 83 0.965 2 93 90 0.968 3 136 129 0.949 4 82 70 0.854 Total

Dear all, bind cannot load zone files, i see the following in the log for all domains. i have tried with 777 to all zone files but does not make difference. main directories are also owned by named:named May 21 15:45:23 nsfo1 named-sdb[2482]: zone dot.com.gh/IN: loading from master file master/db.gh.com.dot failed: permission denied help would be very much appreciated Thanks

On Mon, Dec 3, 2012 at 11:09 PM, gamma_chen <gamma_chen at yahoo.com.tw> wrote: > I add the empty file .nojekyll to https://github.com/Jonathan2251/lbd/tree/gh-pages and do "make gh-pages" again. It's the same. I check the uvbook has the gh-pages as mine as follows, The book looks really good at <http://jonathan2251.github.com/lbd/>! I think if you change the PDF

Dear R-users, I am very new to R, so maybe my question is very easy to answer. I have the following table: TAB1<-data.frame(Name,Number), "Name" and "Number" are all character strings, it looks like this: Name Number ab 2 ab 2 NA 15 NA 15 NA 15 cd 3 ef 1 NA 15 NA 15 gh 15 gh 15 I want to delete all the rows

Hello, The question looks like simple. It's probably even stupid. But I spent several hours searching Internet, downloaded tons of papers, where deviance is mentioned and... And haven't found an answer. Well, it is clear for me the using of entropy when I split some node of a classification tree. The sense is clear, because entropy is an old good measure of how uniform is distribution.

Greetings. 2.15.0 added this behavior http://developer.r-project.org/blosxom.cgi/R-devel/NEWS/2012/03/29#n2012-03-29 o Passing R objects other than atomic vectors, functions, lists and environments to .C() is now deprecated and will give a warning. Most cases (especially NULL) are actually coding errors. NULL will be disallowed in future. This seems to make sense,

Thank you - that makes sense. I've seen something about swapping and optimizing channels on the console, but I didn't realise "optimize" meant "not do what you wanted". OK, so here's why I'm dialling anything at all: The first dial is because I MUST limit the incoming call to less than 60 minutes. The second dial, which carries the gH option, is because I

Hi, I have a vector with 4 elements, e.g., tau_i=c(100,200,300,400), but potentially tau_i[0]=0. In a "for" loop, tau_i=c(100,200,300,400) m=4 tau_i[0]=0 # <------- ? P_i=1 for(i in 2:m) { P_i = P_i*(tau_i[i-1]-tau_i[i-2]) } Error in P_i = P_i * (tau_i[k - 1] - tau_i[k - 2]): replacement has length zero Unfortunately, I can add this potential element into