similar to: minimizing multivariate functions

Hello, I want to use optim with simulated data. My program generates a sample and then calls optim to optimise a function that depends on the generated sample. I want to do this with 100 samples. But optim does not 'work' for the 100 samples, sometimes I get a message saying that sigma matrix was not invertible and therefore the optimisation was not carried out, and then my program

Hi all, Hope someone could help me. I am trying to run automatically the conversion of an Rwn file to a tex file. I am using windows 7, and cygwin. I tried to run automatically the Sweave.sh script, in its the most recent version available at R webpage: http://cran.r-project.org/contrib/extra/scripts/Sweave.sh Unfortunately, I got this error message: =========================== Raquel at

Hi all, I would like to compute ergodic mean using MCMC output from WinBUGS. I tried using CODA package, but it seems that it is not implemented yet. Could anyone help me to compute this? Attached to this email are my output and index files. Kind regards, Raquel -- Raquel Rangel de Meireles Guimar?es Doutoranda em Demografia raquel at cedeplar.ufmg.br

Hi all, I have two basic questions, hope you should help me: 1. How do I save a log file in R with the results? For example, in Stata it can be done by using "log using c:\...\test.txt" 2. How do I display the execution time of one function like the one below: CalculaCorrelacao <- function(construto, n) { library(polycor) for (i in 2:n){ for (j in i:n+1){ x1

Hi, I'm running a tutorial ("Meta-analyses of data from two (or more) microarray data sets"), which use wgcna package. I have an error in the function modulePreservation (it is below). I'm using R2.13 Can you help me? Do you know, what is happens? Thanks Raquel multiExpr = list(A = list(data=t(badea)),B = list(data=t(mayo))) # two independent datasets (dim = 13447 x 36) mp =

hi list, i have to ask you again, having tried and searched for several days... i want to do a TukeyHSD after an Anova, and want to get the adjusted p-values after the Tukey Correction. i found the p.adjust function, but it can only correct for "holm", "hochberg", bonferroni", but not "Tukey". Is it not possbile to get adjusted p-values after

While I don't know anything about Box's M test, I googled around and found a Matlab M-file that computes it. Below is my straight-forward translation of the code, without knowing Matlab or the formula (and done in a few minutes). I hope this demonstrates one of Prof. Ripley's point: If you really want to shoot yourself in the foot, you can probably program R to do that for you. [BTW:

The EuroLLVM'18 Developers’ Meeting relies on volunteers to keep things running smoothly. Moderators are critical to this as they keep speakers on track and facilitate Q&A after the talk. I’m looking for community members who would be attending specific talks anyway, to volunteer to moderate the sessions. If you are interested in volunteering, *please respond to this email with your first

I''ve run across something that would save me days of trouble. Situation, the contents of one ZFS file system needs to be moved to another ZFS file system. The destination can be the same Zpool, even a brand new ZFS file system. A command to move the data from one ZFS file system to another, WITHOUT COPYING, would be nice. At present, the data is almost 1TB. Ideally a "zmv" or

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello: My name is Ruben and I am from Elche (Spain) We are trying to setup a Samba server in our network. Our server is Suse Linux 7.3 Professional, with this version of samba: at:/ # smbd -V Version 2.2.1a at:/ # nmbd -V Version 2.2.1a I think that smbd is running OK becouse we can access to the server using its IP from the outside. But Server

Hello, I am getting a GPG error with the ubuntu repository at the bristol UK mirror. When my source.list has this line: deb http://www.stats.bris.ac.uk/R/bin/linux/ubuntu/ intrepid/ On an "apt-get update" you get this: W: GPG error: http://www.stats.bris.ac.uk intrepid/ Release: The following signatures were invalid: BADSIG D67FC6EAE2A11821 Vincent Goulet <vincent.goulet at

Hi all, I am very new to R and this might be a simple question but I have looked everywhere you suggest before writing to you. I am trying to change font type from san-serif to a serif (Times New Romans) on all labels and axis of my boxplot. I have used this function in other plots before, e.g.: plot(residuals~lnlifespan, data=mydata, pch=psymb, font=6, xlab="ln reproductive

Dear List, I'm having a curious problem with lapply(). I've used it before to convert a subset of columns in my dataframe, to factors, and its worked. But now, on re-running the identical code as before it just doesn't convert the columns into factors at all. As far as I can see I've done nothing different, and its strange that it shouldn't do the action. Has anybody

Dear R users, I am trying to fit an lmer model with only random effects which is giving me the following error: Error : length(f1) == length(f2) is not TRUE In addition: Warning messages: 1: In P1L55:family : numerical expression has 390 elements: only the first used 2: In P1L55:family : numerical expression has 390 elements: only the first used I am trying to extract variance components

Dear list, I have been using multivariate multiple regression (MMR) in the form lm(Y~X) where Y and X are matrices of response and predictor variables. I know that summary(mlm.object) would give the usual lm statistics for each response variable separately and that anova.mlm(mlm.object) will give the analysis of variance table of the mlm object. However, anova.mlm (also manova(mlm.object))

hello. I have a query about the Matrix package for R. I wrote some code a while ago using the Matrix package version 1.6.2 with an early version of R, to do some linear least squares for regression: xn [,1] [,2] [,3] [,4] [1,] 1 0.7429352 0.5519528 0.4100652 [2,] 1 0.7443713 0.5540886 0.4124477 [3,] 1 0.7447385 0.5546355 0.4130584 [4,] 1 0.7459597 0.5564558

Hi all, There are quite a few of us working on LLVM in the Bristol area. I invite you all to the first Bristol LLVM social: Where: The Hope & Anchor 38 Jacobs Wells Road Bristol BS8 1DR http://www.hopeandanchor.net When: Thursday 20 September from 7pm Look forward to meeting you all there. Jeremy -- Tel:

Hello, <http://tolstoy.newcastle.edu.au/R/help/04/07/0312.html#0313qlink1> Could anybody tell me how to generate discrete Laplacian distribution? I need to sample uma discretised Laplacian density like this: J( g -> g´) ~ exp (-lambda | g´ - g |) g in {0,…, gmax} Thanks, Nicolette [[alternative HTML version deleted]]

Hi Jay/Samba peeps, Emailing in reference to http://lists.samba.org/archive/samba/2011-April/162277.html I have seen a very similar issue with a similar setup. Users fail to be verified with: getent passwd username Entry in the log at same time is: [2012/01/23 16:58:53.159761, 3] winbindd/winbindd_misc.c:352(winbindd_interface_version) [18510]: request interface version [2012/01/23

Dear All I am trying to do a repeated measures analysis using lmer and have a number of issues. I have non-orthogonal, unbalanced data. Count data was obtained over 10 months for three treatments, which were arranged into 6 blocks. Treatment is not nested in Block but crossed, as I originally designed an orthogonal, balanced experiment but subsequently lost a treatment from 2 blocks. My