Comments (2)
Amazing answer, thank you. I like the coord_trans option a lot!I found two more interesting tweaks that I'll leave here as a note to others and probably my future self.
x <- data.frame(x=rlnorm(1000, 1,1))
# stat(density) achieves the same as stat(count*.1) for obvious reasons
ggplot(x, aes(x,stat(density)))+
geom_histogram(binwidth = 0.1 )+
stat_theodensity(distri = "lnorm") +
coord_trans(x = "log10")
# manually supplying breaks solves the xmin issue.
# If desired, one can make them logarithmic so they appear to have
# equal width in final plot:
logseq <- function(x, length = 100, abs_min=1e-5){
# exponentially spaces sequence covering the range of x (+/- ten percent)
pmax(abs_min, exp(seq(log(.9*min(x)), log(1.1*max(x)), length.out = length)))
}
ggplot(x, aes(x,stat(density)))+
geom_histogram(breaks = logseq(x) )+
stat_theodensity(distri = "lnorm") +
coord_trans(x = "log10")
This solves the issue, as far as I'm concerned. Thanks again for the nice work and helpful package!
from ggh4x.
Hi, thank you for bringing this to my attention.
I can indeed confirm that the scale_x_log10()
produces an error. This happens for understandable reasons; scale transformations occur prior to the calculation of any statistic, so the statistic is calculated on the transformed data. The log10 transformation makes the lognormal distribution a normal distribution with negative values, which the fitdistrplus::fitdist()
function (that estimates the parameterisation under the hood) can't handle for lognormal distributions.
For this I can see 2 possible solutions:
- Simply use the normal distribution if you use the log10 scale; but this work only for normal-lognormal scale/stat pairs.
- I should make an option to untransform the data prior to fitting the data and then re-transform the result.
With regards to the coordinate transformation, I don't think it is stat_theodensity()
that is causing the problems, but the histogram. If we compare the following 3 plots, the histogram fails when the binwidth is wide.
df <- data.frame(x = rlnorm(1000, 1, 1))
ggplot(df, aes(x)) +
geom_histogram(bins = 100) +
coord_trans(x = "log10") # Error
ggplot(df, aes(x)) +
geom_histogram(binwidth = 0.1) +
coord_trans(x = "log10") # Works
ggplot(df, aes(x)) +
stat_theodensity(distri = "lnorm") +
coord_trans(x = "log10") # Works
Because histograms are parameterised under the hood as rectangles, if the xmin
is below zero, that would produce an error. Maybe this is an issue that might be raised at ggplot2 itself? The solution I'd put forward is to choose a smaller binwidth for now.
As a small aside, the overlay of a histogram and a (scaled) density curve works best if the stat(count)
is multiplied by the width of a histogram's bin (which is variable when you set bins = 100
). I recommend the following for your example:
ggplot(df, aes(x)) +
geom_histogram(binwidth = 0.1) +
stat_theodensity(aes(y = stat(count * 0.1)),
distri = "lnorm") +
coord_trans(x = "log10")
Do you think the untransform option would solve the issue?
Best wishes
from ggh4x.
Related Issues (20)
- Cannot use `trans` after using `facetted_pos_scales` HOT 3
- Release ggh4x 0.2.7
- facet_nested have a different behaviour as compare to facet_grid when combining character and numeric HOT 2
- facet_wrap2 adding wrong x axis labels HOT 3
- Facet size = "variable" not working? HOT 2
- Release ggh4x 0.2.8
- match facet data with ggbarplot HOT 4
- Missing data prevents setting up secondary y-axis scale in help_sec_range and help_sec_max HOT 1
- Ggplotly will move x-axis back to default position HOT 1
- Add tag_facet function? HOT 5
- Error when facetting and using `ggtext` HOT 3
- geom_rectmargin won't work bug with fill_alpha HOT 3
- ggtext::element_markdown not working with ggh4x_0.2.8 HOT 3
- Remove x-axis labels with independent = "x" if they are identical HOT 4
- Add compatibility with ggtext::element_markdown() HOT 1
- Error in check_object(cols, is_quosures, "a {.fn vars} specification" HOT 3
- facet_grid with remove_labels removes labels but leaves the space between the panels HOT 2
- `guide_stringlegend()` throws `Error in get_layer_key(...) : argument "data" is missing, no default` HOT 4
- S3 guide system deprecated in ggplot2 3.5.0 HOT 2
- Compatibility with ggtext::element_markdown() issue HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from ggh4x.