similar to: Maximum Likelihood using optim()

Dear Mailing List Members, the problem I've been grappling with für quite some time now is the following: I have a 100 rows x 200 columns matrix. data.set <- matrix(rnorm(20000, 100, 200)) Now I would like to get a vector of length 100 which collects the values from the following procedure: Take the sum of the minima of the two values from each row of columns 1 and 101, and divide it

Dear mailing list, how can I identify all those rows of matrix B which fulfill some condition based on another matrix A? More precisely, A <- matrix(c(1, 1, 0, 1, -9, 1, -9, 1, 0, 0, 1, 1, 0, -9, 1, 0, -9, 0, 1, 1, 1, -9, 1, 1, 1), ncol = 5, byrow = TRUE) B <- matrix(c(0,0,1,0,0,0,1,1,0,0,0,0,1,1,0,1,0,1,1,0,0,1,1,1,0,1,1,1,1,0,0,0,1,0,1,1,

Hi, I have the following problem: I am working on assessing the accuracy of diagnostic tests. I have a log likelihood and 3 unknown parameters. Now I want to apply the EM algorithm to find the maximum likelihood estimates of these parameters. But I don't know anaything about how to code that in R. Can somebody help me?? Greats, Elena -- View this message in context:

Dear All, I have just uploaded a new version of the QCA (Qualitative Comparative Analysis) on CRAN, and it will be propagated in a couple of days. This is version 1.0-0 ("Easter edition") of the package, straight from the previous version 0.6-5, and it represent a major re-write of the package in order to accomodate fuzzy-sets. I am pleased to welcome Alrik Thiem as a co-author of this

Dear All, I have just uploaded a new version of the QCA (Qualitative Comparative Analysis) on CRAN, and it will be propagated in a couple of days. This is version 1.0-0 ("Easter edition") of the package, straight from the previous version 0.6-5, and it represent a major re-write of the package in order to accomodate fuzzy-sets. I am pleased to welcome Alrik Thiem as a co-author of this

Hi everyone, I am looking to do some manual maximum likelihood estimation in R. I have done a lot of work in Stata and so I have been using output comparisons to get a handle on what is happening. I estimated a simple linear model in R with lm() and also my own maximum likelihood program. I then compared the output with Stata. Two things jumped out at me. Firstly, in Stata my coefficient

Dear R community, I tried to optimized my R code by using Rprof. In my R code I'm using MySQL database connections intensively. After a bunch of queries R fails with the following error message: Error in .Call("RS_MySQL_newConnection", drvId, con.params, groups, PACKAGE = .MySQLPkgName) : RS-DBI driver: (could not connect mylogin@mydatabase on dbname "myDB"

is there a way to use the actual index value for plotting zoo objects this is the way that the index is set up and a sample range of what I would like to plot 01/01/06 00:00:00 - 01/01/06 23:45:00 { library(zoo) # chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) }} x <- structure(c(15.57, 15.5,

Hello, according to documentation: https://wiki.samba.org/index.php/Setting_up_a_Share_Using_Windows_ACLs User is required to setup the file share using ms windows (Computer Management->Connect to another computer). That I am trying without success. When I connect with Windows 7 client I get error message: "Event Viewer cannot connect to computer 'SAMDOM.EXAMPLE.COM. The error ported

The current llvm/llvm-gcc-4.2 2.6 branch passes all of the Polyhedron 2005 benchmarks built with its gfortran. The results compare as follows... Compile Command : gfortran -ffast-math -funroll-loops -msse3 -O3 %n.f90 -o %n benchmark gcc-4.2.4 llvm-gcc-svn llvm-gcc-2.6 llvm-gcc-2.6 at -m32 20081031 -m32 at -m32 at -m64 ac 18.30

I checked my entire setup again and again. Now I can see the my share. The most relevant change I did (afaik - ofcourse): - since I am using ad backend I created dedicated unix admin groups (as recommended in the docs) - I removed the $ from the share section/defnition ([projects$] -> [projects]) - I changed browseable switch from 'no' to 'yes' Iam still getting the "event

Hello everybody, if I'm using par(xpd=NA) gridlines will plotted outside the plotting region. This is "new" (and, in my opinion, unaesthetic) in version 1.4 - compared to version 1.3. Is this a bug or a feature? Example: ------------------ par(xpd=NA) plot(1:11) grid() ------------------ platform sparc-sun-solaris2.6 arch sparc os solaris2.6

Dear R users, I have issues regarding latex() from Hmisc and numSummary() from Rcmdr. Here's an example: > library(Rcmdr) > data(Angell, package="car") > numSummary(Angell[,"hetero"], statistics=c("mean", "sd", "quantiles"), quantiles=c( 0,.25,.5,.75,1 )) > .numSummary <- popOutput() > latex(as.table(.numSummary$table),

Hello everybody, could anybody give me a hint, who I can use rbind on a list of data.frames, please? I have a list with a large number of data.frames of the same structure, like: LIST <- list(X1=data.frame(a=1,b=2), X2=data.frame(a=3,b=4), X3=data.frame(a=5,b=6), ...., XN=data.frame(a=i,b=k)) I would like to bind all data.frames very fast to a single data.frame, something like that: DF <-

Hello everybody, if I'm running a R-script on a command line R-session which uses the function "try" and I'd like to interrupt the execution of my R-script pressing "CTRL-C" than R ignores the first "CTRL-C" key stroke and after I pressed "CTRL-C" second time R by itself is interrupted and I fall back to the (UNIX-) command prompt. Without using

Hello everybody, if I try to (r)bind a number of large dataframes I run out of memory because R wastes memory and seems to "forget" to release memory. For example I have 10 files. Each file contains a large dataframe "ds" (3500 cols by 800 rows) which needs ~20 MB RAM if it is loaded as the only object. Now I try to bind all data frames to a large one and need more than

Hello everybody, I found that queries (send by "dbExecStatement" ) with more than 4000 characters length produces an error in ROracle (ver 0.3-3). Maybe there is a limitation of 4kB.... Is this a bug? If yes, is this problem solved in the latest version of ROracle (ver 0.5-0)? My system information: platform sparc-sun-solaris2.8 arch sparc os solaris2.8

Hello everybody, Could anybody give me a hint how I can use "layout" and "filled.contour" (or "image" plot with a color legend) together, please? What I want to do is something like the following example (Two or more plots with a legend for each at one page): data(volcano) layout(matrix(1:2, 1, 2, byrow = TRUE)) for (i in 1:2) { filled.contour(i*volcano, color =

Hello everybody, is there a possibility to install and use multiple versions of a R library at the same time? A few words regarding the background of this question: If I change functions of an already developed package, I have to verify complex functions in a lot of R applications. So it would be useful to use different versions of a function/library and switch to the newer function/library

Hello everybody, I have following string type and try to convert it to float: > a<-"20020104.123456" > str(as.numeric(a)) num 20020104 > str(as.real(a)) num 20020104 > str(as.double(a)) num 20020104 How can I get an float instead of an integer? How can I set the desired precision? Thanks in advance Lutz Thieme Lutz Thieme AMD Saxony Manfacturing GmbH Product