r - Graphing a histogram overlaid with a fitted 2 parameter Weibull function -

-

i plot both histogram fitted weibull function on same graph. code plot histogram is:

    hist(data$grddia2, prob=true,breaks=5)

the code fitted weibull function is:(need mass package)

    fitdistr(data$grddia2,densfun=dweibull,start=list(scale=1,shape=2))

how plot both on same graph. i've attached data set.

also, bonus can provide code can achieve same thing, create graph each column of data. many columns within data set. nice have graphs on same page.

https://www.dropbox.com/s/ra9c2kkk49vyyyc/diameter%20distribution.csv?dl=0

here code

library("ggplot2") library("dplyr") library("tidyr") library("mass")  # import dataset , filter column "treeno" # use namespace dplyr:: explicitly because of conflict mass:: function "select" data <- read.csv("diameter distribution.csv") %>%    dplyr::select(-treeno)  # function provide weibull distribution each column # distribution calculated based on estimated scale , shape parameters of input fitweibull <- function(column) {   x <- seq(0,7,by=0.01)   fitparam <- column %>%     unlist %>%      fitdistr(densfun=dweibull,start=list(scale=1,shape=2))   return(dweibull(x, scale=fitparam$estimate[1], shape=fitparam$estimate[2])) }  # apply function each column consolidate in data.frame fitdata <-data %>%   apply(2, as.list) %>%    lapply(fun = fitweibull) %>%    data.frame()  # display graphs multiplyingfactor<-10 ggplot() +   geom_histogram(data=gather(data), aes(x=value, group=key, fill=key), alpha=0.2) +   geom_line(data=gather(fitdata), aes(x=rep(seq(0,7,by=0.01),ncol(fitdata)), y=multiplyingfactor*value, group=key, color=key))

and output figure enter image description here

variant: wonderful ggplot2 package can have graphs apart adding final line of code

+ facet_wrap(~ key) + theme(legend.position = "none")

which gives other figure: enter image description here


Comments

Post a Comment

Popular posts from this blog

c++ - No viable overloaded operator for references a map -

-
i trying use map can make tag names reference number. when try , use it, in code error (6 in total every time reference map): src/main.cpp:25:45: error: no viable overloaded operator[] type 'std::map<std::string, std::string>' const char* idcs = node.child_value(tagmap[3]); this code: #include "pugi/pugixml.hpp" #include <iostream> #include <string> #include <map> int main() { pugi::xml_document doca, docb; std::map<std::string, pugi::xml_node> mapa, mapb; std::map<std::string, std::string> tagmap {{"1", "data"}, {"2", "entry"}, {"3", "id"}, {"4", "content"}}; if (!doca.load_file("a.xml") || !docb.load_file("b.xml")) { std::cout << "can't find input files"; return 1; } (auto& node: doca.child(tagmap[1]).children(tagmap[2])) { co
Read more

java - Custom OutputStreamAppender not run: LOGBACK: No context given for <MYAPPENDER> -

-
i wish write custom appender on basis of outputstreamappender . wrote following class package tests; import java.io.ioexception; import java.io.outputstream; import ch.qos.logback.core.outputstreamappender; public class myappender<e> extends outputstreamappender<e> { public myappender() { system.out.println("myappender created"); setoutputstream(new outputstream() { @override public void write(int b) throws ioexception { throw new unsupportedoperationexception(); } }); } } then wrote following runner: package tests; import org.slf4j.loggerfactory; import ch.qos.logback.classic.logger; public class runner { private static final logger log = (logger) loggerfactory.getlogger(runner.class); public static void main(string[] args) { //log.addappender(new myappender<iloggingevent>()); log.info("hello world"); } } finall
Read more

java - Cannot secure connection using TLS -

-
i'm trying establish secure socket connection between java client applet (built jdk 1.7.0_75-b13) , vc++ server application. as test vehicle, used vc++ client/server sample found in msdn forums, modified use schannel , able establish socket using cipher suite tls_rsa_with_aes_128_cbc_sha. works of tls 1.0/1.1/1.2. when try opening socket java applet same server application, connection rejected server reporting following: tls 1.0 acceptsecuritycontext failed: 0x80090327 tls 1.1 acceptsecuritycontext failed: 0x80090331 tls 1.2 acceptsecuritycontext failed: 0x80090331 this java code used create socket: debugprint("setting secure connection"); sslsocketfactory sslsocketfactory = (sslsocketfactory) sslsocketfactory.getdefault(); sslsocket sslsocket = (sslsocket) sslsocketfactory.createsocket("127.0.0.1", socketnumber); debugprint("starting handshake"); sslsocket.settcpnodelay(true); sslsocket.setsolinger(false, 0); sslsocket.s
Read more