similar to: Determining sample size from power function

Hi, With the following Sweave minimal file: ---<--------------------cut here---------------start------------------->--- \documentclass{article} \usepackage{Sweave} \begin{document} <<binom-sim>>= thetas <- seq(0, 1, by=0.001) prior <- rep(1, length(thetas)) / length(thetas) lik <- dbinom(1, 1, thetas) lik.p <- prior * lik post <- lik.p / sum(lik.p)

Here's some code to make lissajous dance. I've attached a small sample GIF. Cheers, Rob Steele robsteele at yahoo dot com plot.lissajous = function(omega.x, omega.y, delta = 0, num.thetas = 200) { thetas = seq(0, 2 * pi, length = num.thetas) xs = sin(omega.x * thetas + delta) ys = cos(omega.y * thetas) plot(xs, ys, type = 'l', lwd = 3, ann = FALSE, axes = FALSE) }

Dear all, I am trying to optimize a logistic function using optim, inside the following functions: #Estimating a and b from thetas and outcomes by ML IRT.estimate.abFromThetaX <- function(t, X, inits, lw=c(-Inf,-Inf), up=rep(Inf,2)){ optRes <- optim(inits, method="L-BFGS-B", fn=IRT.llZetaLambdaCorrNan, gr=IRT.gradZL, lower=lw, upper=up, t=t, X=X)

Hello, I am trying to use optim() on a function involving a summation. My function basically is a thinned poisson likelihood. I have two parameters and in most cases optim() does a fine job of getting the minima. I am simulating my data based on pre specified parameters, so I know what I should be getting. However when my true parameters fall in a particular range, optim() gives

Hi all, I'm making an animation of a rotating persp plot of some data. It's easy, really: I just plot the persp in a loop over the theta parameter, saving each plot, and then chain them together with ImageMagick (as suggested on Paul Johnson's R-tips page). But... when you run the animation, because each graph is individually re-scaled, the image is not the same size in each

On Thu, 2003-04-03 at 12:59, Kaveh Goudarzi wrote: > Hello People, > > Are there any ongoing efforts to port speex enc/dec > to java? I know there are a couple of things lying about in > jcraft and other sites but they all seem a little half hearted. You might like to write a JNI Wrapper. Cheers David > I've used jspxdec.zip written by james.lawrence@hellonetwork.com

Hi all, I was wondering whether it is possible to use the lapply() function to alter the value of the input, something in the spirit of : a1<-runif(100) a2<-function(i){ a1[i]<-a1[i-1]*a1[i];a1[i] } a3<-lapply(2:100,a2) Something akin to a for() loop, but using the lapply() infrastructure. I haven't been able to get rapply() to do this. The reason is that the "real"

> Thanks, I'm having a look at it. I've also had a look at the > various files making up libspeex trying to see how hard the port to > java would be ... Having only had a cursory look it looks ok except > for filters_sse.h & lpt_sse.h which have some assembly. > > Jean-Marc can you throw some light on what these do ... Also > anyone interested in chiping in some

Dear All, I've got a weird problem with my asterisk box which has fxo interfaces (TDM400). Well, the problem is that the interface answers the call, but no caller id is being received. Also, sometimes this error happens: fsk_serie made mylen < 0 Any idea what is going on? Thanks, Kaveh __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has

Hello Mark: Count me in. I'm an excellent C programmer and I've done assembly coding. Although I've not done float to int conversions, I'll pick it up quickly ... send me an example, please. Also, I direct The Tech Startup Connection here in Silicon Valley, so I have a very large resource pool. Check out

Dear R-help, I am trying to estimate a Cox model with nested effects basing on the minimization of the overall AIC; I have two frailties terms, both gamma distributed. There is a error message (theta2 argument misses) and I don?t understand why. I would like to know what I have wrong. Thank you very much for your time. fitM7 <- coxph(Surv(lifespan,censured) ~ south + frailty(id,

Hi, here comes my problem, say I have the following functions (example case) #------------------------------------------------------------ function1 <- function (x, theta) {a <- theta[1] ( 1 - exp(-theta[2]) ) * theta[3] ) b <- x * theta[1] / theta[3]^2 return( list( a = a, b = b )) } #----------------------------------------------------------- function2<-function (x, theta) {P

Hi, I got my hands on the LCA 2003 CD today and I tried to listen to it on my powerbook running Linux. Instead of getting speech, I got an earfull of noise! Attempting the same experiment on an i386 with the same CD produced understandable speech. It seems that the output routine in speexdec converts the output data to little endian short format, which is incorrect on the ppc architecture. The

Dear All, I'm trying to implement the section of the manual pertaining to the gcc-ubsan test carried by CRAN on my local computer (ubuntu 14.04): http://www.stats.ox.ac.uk/pub/bdr/memtests/gcc-UBSAN/README.txt I was wondering whether someone could tell what the minimal version of the gcc tool chain needed to run the gcc-ASAN and gcc-UBSAN alignment tests on ones local

Dear List, I'd like to call pyhton function from within R. I tried installing the latest version of RSPython: wget http://www.omegahat.org/RSPython/RSPython_0.7-1.tar.gz R CMD INSTALL --clean RSPython_0.7-1.tar.gz I get a compile error (posted below). Did anyone else run against this ? Is there a solution ? checking for python... /usr/bin/python Python version 2.6 Using threads checking

This is not so much a question as a contribution, but comments are welcome. Comments: 1) thank you very much to Paul Smith in the post https://stat.ethz.ch/pipermail/r-help/2008-March/157249.html This is intended to build on that example with something more complex than a 2x2 set of constraints 2) "L-BFGS-B" does not appear to work in optimConst Problem: let's

Dear all, I was wondering if there is any way i could do a "Grid Search" on a parameter space using R (as SAS 6.12 and higher can do it) to start the Newton-Gauss Linearization least squares method when i have NO prior information about the parameter. W. N. Venables and B. D. Ripley (2002) "Modern Applied Statistics with S", 4 th ed., page 216-7 has a topic

Hi list, in the Rglpk_solve_LP function (::Rglpk), on line 26, the function calls a function as.glp_bounds() that i cannot access. i'm trying to alter the Rglpk_solve_LP function to add a line to retrieve column/row dual values. everytime i change the slightest line of code inside Rglpk_solve_LP() [to even add a print] i get a ': could not find function "as.glp_bounds"'

Dear List, I'm currently working on an R package that works great when compiled with GCC 4.8 under linux but which doesn't seem to work so well when compiled by the win builder (http://win-builder.r-project.org/). After asking around a bit, the primary culprit could be the older version of GCC used in Rtools. Now, for academic purposes, I would like to have a zip version of this

Hi all, were can i find the sources for the function sample() (the one from the base package) thanks in advance, -- View this message in context: http://r.789695.n4.nabble.com/source-for-function-sample-in-package-base-tp4447796p4447796.html Sent from the R help mailing list archive at Nabble.com.