r - Graphical output of density for the function gammamixEM (package mixtools) -

-

i'm using function gammamixem package mixtools. how can return graphical output of density in function normalmixem (i.e., second plot in plot(...,which=2)) ?

update:

here reproducible example function gammamixem:

x <- c(rgamma(200, shape = 0.2, scale = 14), rgamma(200,       shape = 32, scale = 10), rgamma(200, shape = 5, scale = 6)) out <- gammamixem(x, lambda = c(1, 1, 1)/3, verb = true)

here reproducible example for function normalmixem:

data(faithful) attach(faithful) out <- normalmixem(waiting, arbvar = false, epsilon = 1e-03) plot(out, which=2)

enter image description here

i obtain graphical output of density function gammamixem.

here go.

out <- normalmixem(waiting, arbvar = false, epsilon = 1e-03)  x          <- out  whichplots <- 2 density = 2 %in% whichplots loglik  = 1 %in% whichplots  def.par    <- par(ask=(loglik + density > 1), "mar") # ask , mar changed mix.object <- x   k    <- ncol(mix.object$posterior) x    <- sort(mix.object$x)    <- hist(x, plot = false) maxy <- max(max(a$density), .3989*mix.object$lambda/mix.object$sigma)

i had dig source code of plot.mixem

so, gammamixem:

x <- c(rgamma(200, shape = 0.2, scale = 14), rgamma(200,                                                      shape = 32, scale = 10), rgamma(200, shape = 5, scale = 6)) gammamixem.out <- gammamixem(x, lambda = c(1, 1, 1)/3, verb = true)    mix.object <- gammamixem.out  k    <- ncol(mix.object$posterior) x    <- sort(mix.object$x)    <- hist(x, plot = false) maxy <- max(max(a$density), .3989*mix.object$lambda/mix.object$sigma)  main2 <- "density curves" xlab2 <- "data"  col2  <- 2:(k+1)   hist(x, prob = true, main = main2, xlab = xlab2,       ylim = c(0,maxy))   (i in 1:k) {   lines(x, mix.object$lambda[i] *            dnorm(x,                  sd = sd(x))) }

enter image description here

i believe should pretty straight forward continue example bit, if want add labels, smooth lines, etc. here's source of plot.mixem function.


Comments

Post a Comment

Popular posts from this blog

Spring Boot + JPA + Hibernate: Unable to locate persister -

-
i trying work entity inside request method of spring boot application. entity in project , pom.xml has dependencies it. why error message?: org.hibernate.unknownentitytypeexception: unable locate persister: com.ric.bill.model.bs.par package com.ric.web; import java.util.concurrent.atomic.atomiclong; import javax.persistence.entitymanager; import javax.persistence.persistencecontext; import net.sf.ehcache.cachemanager; import net.sf.ehcache.management.cache; import net.sf.ehcache.management.cachestatistics; import org.hibernate.stat.statistics; import org.springframework.beans.factory.annotation.autowired; import org.springframework.boot.autoconfigure.domain.entityscan; import org.springframework.cache.annotation.enablecaching; import org.springframework.transaction.annotation.enabletransactionmanagement; import org.springframework.transaction.annotation.propagation; import org.springframework.transaction.annotation.transactional; import org.springframework.web.bind.annotati
Read more

go - Golang: panic: runtime error: invalid memory address or nil pointer dereference using bufio.Scanner -

-
i implementing go program uses bufio.scanner , bufio.writer have packaged code follows package main import ( "fmt" "player/command" "strings" ) func main() { //enter code here. read input stdin. print output stdout commands.scanner.scan() { //scan new line , send comand variable check command exist or not input := strings.split(strings.trim(commands.scanner.text(), " "), " ") command := input[0] if command == "" { fmt.printf("$ %s:", commands.pwd) continue } if !commands.commands[command] { commands.throwerror("cannot recognize input.") } else { commands.execute(command, input[1:], nil) } fmt.printf("$ %s:", commands.pwd) } } i using init.go file in main package follows package main import ( "flag" "player/sourc
Read more

c - double free or corruption (fasttop) -

-
this code reads null terminated string file (file containing null terminated strings & should provide index wanted string). error: *** error in `./main': double free or corruption (fasttop): 0x090a4a80 *** aborted code: char *tmp_realloc=null; char *sstring=null; int i=1; /*set file pointer point string*/ fseek(fh,index,seek_set); sstring=(char*)malloc(sizeof(char)); if(sstring == null) return null; tmp_realloc=sstring; while( (*(sstring+i-1)=(char)fgetc((file*)fh)) != 0x0 ) { /*reallocate more memory*/ tmp_realloc=(char*)realloc((char*)sstring,++i); /*if not same address copy old new & set old point new*/ if(tmp_realloc != sstring){ strcpy((char*)tmp_realloc,(char*)sstring); free(sstring); sstring=tmp_realloc; } } could tell me pointer causing issue (& reason) ? remove part /*if not same address copy old new & set old point new*/ if(tmp_realloc != sstring){ strc
Read more