Databases
StackOverflow2013
Postsdbo
POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
 

Note that there are some explanatory texts on larger screens.

plurals
  1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
    primarykey
    data
    text
    <p>I have an R function which produces 95% confidence ellipses for scatterplots. The output looks like this, having a default of 50 points for each ellipse (50 rows):</p> <pre><code> [,1] [,2] [1,] 0.097733810 0.044957994 [2,] 0.084433494 0.050337990 [3,] 0.069746783 0.054891438 </code></pre> <p>I would like to superimpose a number of such ellipses for each level of a factor called 'site' on a <code>ggplot2</code> scatterplot, produced from this command:</p> <pre><code>&gt; plat1 &lt;- ggplot(mapping=aes(shape=site, size=geom), shape=factor(site)); plat1 + geom_point(aes(x=PC1.1,y=PC2.1)) </code></pre> <p>This is run on a dataset, called <code>dflat</code> which looks like this:</p> <pre><code>site geom PC1.1 PC2.1 PC3.1 PC1.2 PC2.2 1 Buhlen 1259.5649 -0.0387975838 -0.022889782 0.01355317 0.008705276 0.02441577 2 Buhlen 653.6607 -0.0009398704 -0.013076251 0.02898955 -0.001345149 0.03133990 </code></pre> <p>The result is fine, but when I try to add the ellipse (let's say for this one site, called "Buhlen"):</p> <pre><code>&gt; plat1 + geom_point(aes(x=PC1.1,y=PC2.1)) + geom_path(data=subset(dflat, site="Buhlen"),mapping=aes(x=ELLI(PC1.1,PC2.1)[,1],y=ELLI(PC1.1,PC2.1)[,2])) </code></pre> <p>I get an error message: <code>"Error in data.frame(x = c(0.0977338099339815, 0.0844334944904515, 0.0697467834016782, : arguments imply differing number of rows: 50, 211</code></p> <p>I've managed to fix this in the past, but I cannot remember how. It seems that geom_path is relying on the same points rather than plotting new ones. Any help would be appreciated.</p>
    singulars
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PTQuestion
    1. This table or related slice is empty.
    1. USradu
    plurals
    1. This table or related slice is empty.
    1. PL
      singulars
      1. LTLinked
    2. PL
      singulars
      1. LTLinked
    3. PL
      singulars
      1. LTLinked
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    2. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. USunutbu
      1. VTFavorite
    2. VO
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. This table or related slice is empty.
      1. VTUpMod
    3. VO
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. This table or related slice is empty.
      1. VTUpMod
    1. CODid you try to change the default of 50 points to 211? Does it work? You might have to add another argument to your function (the number of points)
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. USGeorge Dontas
    2. COHi, thanks for the quick reply. The function can change number of points, and I did try it with 211 points. It produces a strange very thick circle. I think it is not subsetting the data, first of all, and it should be able to plot it with 50 points - at least from the documentation, you can use different datasets on the same plot, so naturally, different numbers of points should be ok too.
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. USradu
    3. COit will be much easier for us if you provide a minimal reproducible example.
      singulars
      1. POHow can a data ellipse be superimposed on a ggplot2 scatterplot?
      1. USxiechao
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload