# VPC plot with observations

**URL:** <https://discourse.pumas.ai/t/vpc-plot-with-observations/640>\
**Category:** Uncategorized\
**Created:** [March 8, 2022, 5:22am UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640 "2022-03-08T05:22:14Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![jayashreed1412](https://yyz2.discourse-cdn.com/flex030/user_avatar/discourse.pumas.ai/jayashreed1412/32/178_2.png) [@jayashreed1412](https://discourse.pumas.ai/u/jayashreed1412)\
**Post date:** [March 8, 2022, 5:22am UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/1 "2022-03-08T05:22:14Z")

</div>

Hi,  
I was trying to plot VPC for my estimation, but in the VPC Plot, there are no observations; only the bandwidth is visible.

This is the code I used to run the VPC.

```auto
vpc1=vpc(est53; 
        samples=500,
        observations=[keys(est53.data[1].observations)[1]],
        quantiles=(0.1,0.5,0.9),
        level=0.95,
        ensemblealg=EnsembleSerial(),
        bandwidth=60,
        smooth=true)

final_vpc=vpc_plot(vpc1)

```

 ![final model vpc](https://canada1.discourse-cdn.com/flex030/uploads/pumas/original/1X/4ffdf89c95117520014bb3929e397152acd70a67.png)

---

<div class="post-metadata">

**Author:** ![bobbrown](https://avatars.discourse-cdn.com/v4/letter/b/c2a13f/32.png) [@bobbrown](https://discourse.pumas.ai/u/bobbrown)\
**Post date:** [March 8, 2022, 11:00am UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/2 "2022-03-08T11:00:28Z")

</div>

Lately there have been several research reports with Visual Predictive Check (VPC). Can anybody provide a layman intro to what these plots are allowing the analyst diagnose?  
Bob

---

<div class="post-metadata">

**Author:** ![patrick](https://yyz2.discourse-cdn.com/flex030/user_avatar/discourse.pumas.ai/patrick/32/206_2.png) [@patrick](https://discourse.pumas.ai/u/patrick)\
**Post date:** [March 8, 2022, 12:06pm UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/3 "2022-03-08T12:06:02Z")

</div>

> [@jayashreed1412](#):
>
> ```auto
> final_vpc=vpc_plot(vpc1)
> 
> ```

You would have to pass`observations=true` to the plotting function.

```auto
final_vpc=vpc_plot(vpc1; observations=true)

```

You can see the options for the plot by writing `?vpc_plot` in the REPL (terminal).

---

<div class="post-metadata">

**Author:** ![patrick](https://yyz2.discourse-cdn.com/flex030/user_avatar/discourse.pumas.ai/patrick/32/206_2.png) [@patrick](https://discourse.pumas.ai/u/patrick)\
**Post date:** [March 8, 2022, 12:18pm UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/4 "2022-03-08T12:18:02Z")

</div>

It’s an interesting question, but also a bit off-topic in terms of the actual question being asked. Could you open a separate question at [https://discourse.pumas.ai/c/diagnostics/vpc/](https://discourse.pumas.ai/c/diagnostics/vpc/). I’m sure someone will be happy to discuss the merits/pros/cons of VPCs or at least provide some references.

Best,  
Patrick

---

<div class="post-metadata">

**Author:** ![bobbrown](https://avatars.discourse-cdn.com/v4/letter/b/c2a13f/32.png) [@bobbrown](https://discourse.pumas.ai/u/bobbrown)\
**Post date:** [March 8, 2022, 2:09pm UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/5 "2022-03-08T14:09:33Z")

</div>

I apologize @patrick.

---

<div class="post-metadata">

**Author:** ![patrick](https://yyz2.discourse-cdn.com/flex030/user_avatar/discourse.pumas.ai/patrick/32/206_2.png) [@patrick](https://discourse.pumas.ai/u/patrick)\
**Post date:** [March 8, 2022, 6:28pm UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/6 "2022-03-08T18:28:18Z")

</div>

No worries, I think it’s an interesting discussion. If you start a new tthread we can discuss there!

---

<div class="post-metadata">

**Author:** ![jayashreed1412](https://yyz2.discourse-cdn.com/flex030/user_avatar/discourse.pumas.ai/jayashreed1412/32/178_2.png) [@jayashreed1412](https://discourse.pumas.ai/u/jayashreed1412)\
**Post date:** [March 9, 2022, 12:33am UTC](https://discourse.pumas.ai/t/vpc-plot-with-observations/640/7 "2022-03-09T00:33:51Z")

</div>

Thank you, Patrick. It works.