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

PHP while loop dynamic rowspan -

-
Image
consider table example using these table wanted print table these adding rowspan item id 1002. here php code $temp_val = ''; $counter = 1; $sql_sel = mysqli_query($con,"select * item_table"); while($res_sel = mysqli_fetch_array($sql_sel)){ if($temp_val == $res_sel['item_id']){ $counter++; echo "<tr></tr>"; } else{ echo "<tr><td rowspan='".$counter."'>".$res_sel['item_id']."</td></tr>"; } $temp_val = $res_sel['item_id']; } echo "</table>"; it's not correct, it's adding rowspan item id 1003 you should switch code inside else , if statement
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

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