similar to: Questions on Random Forest

Hi all, I appreciate the help this list has given me before. I have a question which has been perplexing me. I have been working on doing a Bayesian calculating inserting studies sequentially after using a non-informative prior to get a meta-analysis type result. I created a function using three iterations of this, my code is below. I insert prior mean and precision (I add precision manually

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

Hi, erveryone, I show much thanks to Andy and Matthew on former questions. I now sample only a small segment of a image can segment the image into several classes by RandomForest successfully. Now I have some confusion on it: 1. What is the internal component classifier in RandomForest? Are they the CART implemented in the rpart package? 2. I use training samples to predict new samples. But

Hello, I am trying to do a fit a loglikelihood function with Multivariate distribution via copulas with fitMdvc. The problem is that it doesn't recognize that my beta is a vector of km parameter and when I try to run it it say that the length of my initial values is not the same as the parameter. Can somebody guide me where my mistake is. Thanks, Elisa. #################################

Hi Now a more theoretical question. I have two correlation matrices - one of a set of variables under a particular condition, the other of the same set of variables under a different condition. Is there a statistical test I can use to see if these correlation matrices are "different"? Thanks Mick

I''m trying to create a multi-device filesystem on top of regular files (not actual disks), and mount that to a loopback device. For a filesystem created on a single file, it works fine, but for a filesystem across multiple files, it doesn''t. dd if=/dev/zero of=img1 bs=4096 count=65536 dd if=/dev/zero of=img2 bs=4096 count=65536 dd if=/dev/zero of=img3 bs=4096 count=65536 dd

Hello all, I would like to read the specific line number row of a direct access file (which is stored as a n_row*n_col matrix of elements kind=p) without reading all the preceding lines (i.e 1,2,..,row-1). Is there a function in R that can perform this task? To solve my issue, I tried without to call Fortran from R by doing the following steps: I) I wrote a subroutine in fortran called

Moved to R-help because there's no obvious financial content. Michael On Sat, Apr 6, 2013 at 10:56 AM, Stat Tistician <statisticiangermany at gmail.com> wrote: > Hi, > I want to implement the EM algorithm manually, with my own loops and so. > Afterwards, I want to compare it to the normalmixEM output of mixtools > package. > > Since the notation is very advanced, I

hi im saleh please help me I WANT SOURSE CODE SIMILAR THIS CODE package asl; import java.awt.*; import java.awt.event.* ; import java.applet.*; public class animat1 extends Applet implements Runnable{ int c = 290, d = 300; int x = 20, k = 20, o = 0, m = 10, n = 0, y = 0; Thread t = new Thread(this); Image img1, img2, img3; // AudioClip a1; public void init() {

May I ask why you're doing 512x512 instead of 1024x1024? These are already scaled up coordinates, so 1024x1024 should work no? Or is it because of the seams on the edges? Do those not also appear with 512x512 or does it sample outside of the box? Separately, why not use this approach on nv40 as well? I can't imagine the blitter would be faster... does this result in lower quality? On

Hi, I am trying to fit a function of the form: y = A0 + A1 * exp( -0.5* ( (X - Mu1) / Sigma1 )^2 ) - A2 * exp ( -0.5* ( (X-Mu2)/Sigma2 )^2 ) i.e. a mean term (A0) + a difference between two gaussians. The constraints are A1,A2 >0, Sigma1,Sigma2>0, and usually Sigma2>Sigma1. The plot looks like a "Mexican Hat". I had trouble (poor fits) fitting this function to toy data

Maybe the following article is of interest for you (co-authored by P.J. Rousseeuw, like the book you mentioned): http://www.jstatsoft.org/v01/i04/paper/clus.pdf -Heinrich. > -----Urspr?ngliche Nachricht----- > Von: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] Im Auftrag von Fucang Jia > Gesendet: Donnerstag, 11. Dezember 2003 14:37 > An:

Hi, I would like to simulate data with a binary outcome and a set of predictors that are correlated. I want to be able to fix the number of event (Y=1) vs. non-event (Y=0). Thus, I fix this and then simulate the predictors. I have 2 questions: 1. When the predictors are continuous, I can use mvrnorm(). However, if I have continuous, ordinal and binary predictors, I'm not sure how to simulate

I fitted a gaussian mixture to my financial data. The data can be found here: http://uploadeasy.net/upload/32xzq.rar I look at the density with plot(density(dat),col="red",lwd=2) this has a skew of library(e1071) skewness(dat) -0.1284311 Now, I fit a gaussian mixture according to: f(l)=πϕ(l;μ1,σ21)+(1−π)ϕ(l;μ2,σ22) with:

I need use different parameters of distribution for different case to generate random value, but I use ifelse, the generated value is fixed without change. Here is example data1 y x 1 1 2 2 2 1 3 3 2 4 4 3 5 5 3 6 6 1 7 7 2 8 8 1 9 9 1 10 10 3 11 11 3 12 12 2 ifelse(data1$x==1,rnorm(1,2,1),ifelse(data1$x==2,rnorm(1,-2,1),rnorm(1,110,1))) [1] -1.8042172 0.8478681

Hello, Ok, I will post them in a folder on Google Drive and will give the link on my blog page, -because among other things, on my blog page, I can not display 2 bitmap 512x512 images side by side-... Hope it is ok! Cheers, Raphael 2018-07-15 1:41 GMT+02:00 J.B. Nicholson <jbn at forestfield.org>: > Raphael Canut wrote: > >> I downloaded 20 images from the Internet with

Hi again, I have improved precision of my codec (on the encoder and decoder).I have also improved the -h1 quality setting.-Source code and binaries at http://nhwcodec.blogspot.com/-. I am still trying to improve precision of my codec, with keeping my low-complexity (fast) approach. I do not totally use the reference (and impressive) block prediction with different modes + residual coding scheme,

Hello, Just a quick message to let you know that I get back to NHW Project currently. I downloaded 20 images from the Internet with rather good quality (mainly faces), I compressed them at high compression -l7 setting with the NHW Project, and on 19 images out of 20, I visually prefer the results of NHW compared to x265 (HEVC)!!! -I can make available these 20x3=60 512x512 24bit bitmap images

Hi, I have a bivariate longitudinal dataset. As an example say, i have the data frame with column names var1 var2 Unit time trt (trt represents the treatment) Now suppose I want to fit a joint model of the form for the *i* th unit var1jk = alpha1 + beta1*timejk + gamma1* trtjk + delta1* timejk:trtjk + error1jk var2 = alpha2 + beta2*timejk + gamma2* trtjk + delta2* timejk:trtjk +

I'm running some multi-level binomial models with lme4 and have a question regarding the estimated random effects. Suppose I have nested data e.g. clinic and then patient within clinic. The standard deviations of the random effects at each level are roughly equal in a model for real life data. Attention then turns to examining the individual random effects at each level. I'm extracting