How can I plot a histogram of a long-tailed data using R? -

-

i have data centered in small range (1-10) there significant number of points (say, 10%) in (10-1000). plot histogram data focus on (1-10) show (10-1000) data. log-scale th histogram.

yes, know means not bins of equal size

a simple hist(x) gives alt text while hist(x,breaks=c(0,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,3,4,5,7.5,10,15,20,50,100,200,500,1000,10000))) gives alt text

none of want.

update following answers here produce want (i went continuous plot instead of bar-histogram):

breaks <- c(0,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,4,8) ggplot(t,aes(x)) + geom_histogram(colour="darkblue", size=1, fill="blue") + scale_x_log10('true size/predicted size', breaks = breaks, labels = breaks)![alt text][3]

alt text problem i'd match between scale , actual bars plotted. there 2 options doing : 1 use actual margins of plotted bars (how?) "ugly" x-axis labels 1.1754,1.2985 etc. other, prefer, control actual bins margins used match breaks.

using ggplot2 seems easy option. if want more control on axes , breaks, can following :

edit : new code provided

x <- c(rexp(1000,0.5)+0.5,rexp(100,0.5)*100)  breaks<- c(0,0.1,0.2,0.5,1,2,5,10,20,50,100,200,500,1000,10000) major <- c(0.1,1,10,100,1000,10000)   h <- hist(log10(x),plot=f)   plot(h$mids,h$counts,type="n",       xaxt="n",       xlab="x",ylab="counts",       main="histogram of x",       bg="lightgrey" ) abline(v=log10(breaks),col="lightgrey",lty=2) abline(v=log10(major),col="lightgrey") abline(h=pretty(h$counts),col="lightgrey") plot(h,add=t,freq=t,col="blue") #position of ticks @ <- log10(breaks)  #creation x axis axis(1,at=at,labels=10^at)

this close can ggplot2. putting background grey not straightforward, doable if define rectangle size of plot screen , put background grey.

check functions used, , ?par. allow build own graphs. hope helps.

alt text


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 - UML - How would you draw a try catch in a sequence diagram? -

-
Image
is there standard way of drawing try catch sequence diagrams? i've made attempt based on this feel end result doesn't feel right. code diagram based around: public static void save () { try { filehandle filehandle = gdx.files.external(file); filehandle.writestring(boolean.tostring(constantshandler.soundenabled)+"\n", false); (int = 0; < 5; i++) { filehandle.writestring(integer.tostring(constantshandler.highscores[i])+"\n", true); } } catch (throwable e) { } } my attempt note: know still need add loop. that's ok. fragments meant show conditional control flow in sequence diagram. superstructures puts limitations on use of fragments (see pp. 467 of ss2.1.1). should use critical region try part above. , option catch . however, keep telling can take quite freedom in using uml. it's language , such changes in life time. long reader ge...
Read more

c++ - Gamma correction doesn't look properly corrected, is this linear? -

-
Image
i want implement gamma correction opengl lighting, gamma correction applied, results not seem linear @ all. i found opengl: gamma corrected image doesn't appear linear similar issue, hasn't yet received answer nor discussed actual diffuse lights. as illustration, have following 4 light colors defined in linear space: glm::vec3 lightcolors[] = { glm::vec3(0.25), glm::vec3(0.50), glm::vec3(0.75), glm::vec3(1.00) }; with each light source seperated , basic linear attenuation applied diffuse lighting equation following results: this fragment shader: void main() { vec3 lighting = vec3(0.0); for(int = 0; < 4; ++i) lighting += diffuse(normalize(fs_in.normal), fs_in.fragpos, lightpositions[i], lightcolors[i]); // lighting = pow(lighting, vec3(1.0/2.2)); fragcolor = vec4(lighting, 1.0f); } not barely see difference in brightness gamma corrected lights, attenuation distorted gamma correction. far understanding goe...
Read more