similar to: "Help On optim"

Hello, I am facing a problem with optimization in R from 2-3 weeks. I have some Gaussian mixtures parameters and I want to find the maximum in that *Parameters are in the form * mean1 mean2 mean3 sigma1 sigma2 sigma3 c1 c2 c3 506.8644 672.8448 829.902 61.02859 9.149168 74.84682 0.1241933 0.6329082 0.2428986 I have used optima and optimx to find the

Hello, I am having some problem in retrieving the list elements Code Mean_of_mc is a list as described below mean_mc_1 mean_mc_2 mean_mc_3 mean_mc 0.9577246 2.000516 2.962521 mean_mc 7.957725 9.000516 9.962521 optim(par=c(lVal=2, mu=Mean_of_mc), fn=xnorm, method = "BFGS", control=list(fnscale=-1)) xnorm= function(par) { print(par)

Using R-2.12.1 and latticeExtra-0.6-14, I would like to understand why a rootogram displaying samples from the Poisson distribution looks like I expected it, whereas a rootogram using the normal distribution does not: library(latticeExtra) rootogram(~rpois(1000, lambda = 50), dfun = function(x) dpois(x, lambda = 50)) rootogram(~rnorm(1000), dfun = function(x) dnorm(x,mean(x),sd(x))) I

So, given what I've read in the documentation I would create a text file named document_database.txt that might have the following: remote 192.168.1.10:30000 chert /var/lib/xapian_database/segment1 remote 192.168.1.10:30000 chert /var/lib/xapian_database/segment2 remote 192.168.1.10:30000 chert /var/lib/xapian_database/segment3 etc. I would then in my PHP program open

Greetings, I would like to add 3D arrows (i.e. arrow-headed vectors linking X1Y1Z1 to X2,Y2,Z2) to a 3D plot; ideally the sort of plot that can be rotated interactively. Is this possible using plot3d, or another 3d plotter in R? While it is easy to draw segments in plot3d (e.g. below), I haven't figured out how to add arrow heads, or to create 3d arrows from scratch. ##two headless

There is a piece of source code.On selection of some part of source code and click on button say"Create Segment".Segment should be created. Structure for this is: /*Start Segment1*/ { ------------------------ ----------------------- ---------------------- } /*End Segment1*/ On clicking the create segment button, for the selected text start and end segment should be created.Please help

Hi everyone, I'm trying to fit a of piecewise regression model on a time series. The idea is to divide the series into segments and then to apply linear regression models on each segment but in a "global way" and considering heteroskedasticity between the segments. For example, I build a time series y with 3 segments: segment1=1:20+rnorm(20,0,2) segment2=20-2*1:30+rnorm(30,0,5)

Hi, I want to customise the segments on an xyplot. Below is a simple example of what I'm trying to do... #Example dataset x <- c(-0.25, 0.25, 0.8) y <- c(-0.5, 0, 0.75) gp <- c("A", "I", "C") my.data <- cbind.data.frame(x,y,gp) #setting up the parameters to customise the lines with. ltype <- c(1,2,3) env.col <-

Hi list, anyone know what is going on here? I don't get any sound from the out clip & I get the following when I dial in to asterisk (after which it just times out): -- Event [0=>[06] Ring] on vpb/1-7 -- Executing Wait("vpb/1-7", "2") in new stack Read_channel ## vpb/1-7: Setting record mode, bridge = 0 -- Executing Answer("vpb/1-7",

Hi, gfortran trunk and 9-branch now have an option to automatically generate C prototypes for old-style F77 procedures. I just did for a in *.f; do gfortran -fsyntax-only -fc-prototypes-external $a > ${a%.f}.h; done in the src/modules/lapack directory. This generates header files which contain prototypes like int ilaenv_ (int *ispec, char *name, char *opts, int *n1, int *n2, int *n3,

Hello list, I have a problem with the value of the result size ($msize in omegascript) when collapsing results. The index contains 151452 documents. I'm using Omega 1.0.18 on FreeBSD (I tried both the version in ports and the latest one from xapian.org). This is my indexscript: uniqueid: boolean=Q unique=Q field=uniqueid objectid: field=objectid boolean=XID value=0 objecttype: field=type

I have made many different attempts to run a R program in Batch mode. I run a DOS command like: C:\Progra~1\R\rw1061\bin\rterm.exe -q --no-restore --no-save myprogram.R out.txt What happens: 1. DOS Shell starts up well 2. R starts 3. I get the messages: "ARGUMENT 'myprogram.R' __ignored__" "ARGUMENT 'out.txt' __ignored__" 4. R Shell is ready to

Dear R users, I need to normalize a bunch of row vectors. At a certain point I need to divide a matrix by a vector of norms. I find that the behavior of Matrix objects differs from normal matrix objects. Example the following code examples differ only in xnormed changing from normal to Matrix object: x = matrix(1:12,3,4) x = as(x, "CsparseMatrix") xnorms = sqrt(colSums(x^2))

On Sat, May 04, 2019 at 06:42:47PM +0200, Thomas K?nig wrote: > > > - figure out Fortran2003 specification for C/Fortran interoperability > > -- this _sounds_ like the right solution, but I don't think many > > understand how to use it and what is implied (in particular, will > > it require making changes to LAPACK itself?) > > That would actually be fairly