Skip to contents
library(nhsbsa)
library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union

The NHSBSA Open Data Portal

The NHS Business Services Authority (NHSBSA) Open Data Portal publishes open datasets about NHS activity in England — prescribing, dental, pharmaceutical and contractor data among them. All of it is freely available under the Open Government Licence.

The portal runs on CKAN, a widely used open-source data catalogue. CKAN organises data into datasets (called “packages”) which each contain one or more resources — the individual files (usually CSV) that you can download, and, for tabular resources, query row by row.

nhsbsa is a thin, low-level client for this portal. It wraps the CKAN API and returns plain data for you to work with: tibbles for tabular results and lists for metadata. It contains no knowledge of any particular dataset, so you supply the dataset identifiers and interpret the results yourself. If you are familiar with the CKAN API, the package will feel familiar too: function names and arguments mirror the API.

The API and the response envelope

Every request goes to a CKAN action under

https://opendata.nhsbsa.net/api/3/action/<action>

and comes back as a JSON envelope of the form

{ "success": true, "result": ... }

or, on failure,

{ "success": false, "error": { "message": "..." } }

nhsbsa handles this envelope for you: it checks the success flag, returns the result, and otherwise raises an informative error. If the portal cannot be reached (for example with no internet connection) it fails gracefully with a clear message rather than an obscure low-level error.

The full set of actions is documented in the CKAN Action API reference, and the portal will return the documentation for any individual action it supports, e.g. https://opendata.nhsbsa.net/api/3/action/help_show?name=datastore_search_sql.

What the package wraps

The package wraps the useful read subset of the portal’s actions, in four small groups (the package reference index is organised the same way). If you need an action that is not yet wrapped, please open an issue.

Almost every function maps one-to-one onto a CKAN action of the same name. The two exceptions are the resource helpers, which combine an action with a little extra work and so are not pure wrappers:

(nhsbsa_group_list() is included for completeness; this portal currently defines no groups and so returns an empty vector — it organises data by organisation and tags instead.) Most read-only functions also accept .return_raw = TRUE, which returns the full parsed response envelope instead of the processed result — useful when you need fields the helper does not surface, such as a datastore query’s total.

Working with the returned objects

nhsbsa_package_show(), nhsbsa_resource_show() and nhsbsa_package_search() return potentially large nested lists. To make them easier to scan, they print a tidy summary, and tibble::as_tibble() turns them into a table — a dataset into its resources, and a search into one row per matching dataset:

pkg <- nhsbsa_package_show("english-prescribing-data-epd")
pkg
#> <nhsbsa dataset> "english-prescribing-data-epd"
#> Title: RETIRED - English Prescribing Dataset (EPD)
#> Organisation: Community Prescribing & Dispensing
#> Modified: 2026-03-06
#> Resources: 138
#> Tags: Prescribing, Prescriptions
#> Use `tibble::as_tibble()` for its resources.

tibble::as_tibble(pkg)
#> # A tibble: 138 × 7
#>    name       id                        format created last_modified url   size 
#>    <chr>      <chr>                     <chr>  <chr>   <chr>         <chr> <chr>
#>  1 EPD_201401 8ae6b792-2a0c-4f4b-826c-… CSV    2020-1… NA            http… 6618…
#>  2 EPD_201402 78b8c360-1048-4d30-82fd-… CSV    2020-1… NA            http… 6338…
#>  3 EPD_201403 54584067-3109-4f27-8c48-… CSV    2020-1… NA            http… 6566…
#>  4 EPD_201404 5e25a419-8334-4fb2-bcef-… CSV    2020-1… NA            http… 6501…
#>  5 EPD_201405 5763be28-0dde-430c-bddc-… CSV    2020-1… NA            http… 6623…
#>  6 EPD_201406 32c2c600-2e9d-47e2-9d45-… CSV    2020-1… NA            http… 6566…
#>  7 EPD_201407 9f185ff8-aad4-4555-b1fe-… CSV    2020-1… NA            http… 6716…
#>  8 EPD_201408 3afbd596-a511-4c4f-9a4e-… CSV    2020-1… NA            http… 6429…
#>  9 EPD_201409 d6fa3292-0ba8-468e-ae2d-… CSV    2020-1… NA            http… 6609…
#> 10 EPD_201410 e131e18d-9561-4ba2-98aa-… CSV    2020-1… NA            http… 6696…
#> # ℹ 128 more rows

They are still plain lists underneath, so $, [[ and str() work as usual, and .return_raw = TRUE (or unclass()) gives the unclassed list.

tibble::as_tibble() on a dataset gives the same table as nhsbsa_list_resources() — the difference is just the entry point: nhsbsa_list_resources(id, pattern) fetches (and optionally filters) in one call, while as_tibble() reuses metadata you have already fetched, avoiding a second request.

From the portal website to the API

It helps to think of the package as a programmatic version of the portal website. The things you click there map onto API calls:

  • Browsing all datasets is nhsbsa_package_list() (identifiers) or nhsbsa_package_search() (richer records, with searching and paging).

  • Clicking a tag, say #Prescribing, takes the website to /dataset/?tags=Prescribing. A filter query on the tags field returns the datasets carrying that tag:

    nhsbsa_package_search(fq = 'tags:"Prescribing"')$count
    #> [1] 17

    This returns more datasets than the website shows for the same tag, because the website’s dataset view hides the Freedom of Information disclosure log by default. To match what the website displays, exclude that organisation as well. tibble::as_tibble() turns the result into one row per dataset, so you can pull out the titles to compare against the website directly:

    nhsbsa_package_search(
      fq = 'tags:"Prescribing" -organization:freedom-of-information-disclosure-log'
    ) |>
      as_tibble() |>
      pull(title)
    #> [1] "English Prescribing Dataset (EPD) with SNOMED Code"          
    #> [2] "Prescription Cost Analysis (PCA) Monthly Administrative Data"
    #> [3] "Prescription Cost Analysis (PCA) Annual Statistics"          
    #> [4] "Missing Scottish Dispensing Data from June 2023 to June 2024"
    #> [5] "RETIRED - English Prescribing Dataset (EPD)"

    Tags are case-sensitive; list them with nhsbsa_tag_list(), and the organisations you can filter on with nhsbsa_organization_list().

  • Opening a dataset’s page (e.g. /dataset/english-prescribing-data-epd) corresponds to nhsbsa_package_show("english-prescribing-data-epd"), whose resources element lists the files shown on that page. nhsbsa_list_resources() tidies those resources into a tibble.

  • The “Download” button on a resource fetches the file that nhsbsa_download_resource() streams to disk; the data preview / “Data API” for a resource is backed by the datastore that nhsbsa_datastore_search() queries.

A worked example

Find a dataset

nhsbsa_package_list() returns the identifier of every dataset; use it when you want to scan or search the ids yourself:

datasets <- nhsbsa_package_list()
length(datasets)
#> [1] 2218
head(datasets)
#> [1] "03449"                                           
#> [2] "03500"                                           
#> [3] "baby-loss-certificate-key-performance-indicators"
#> [4] "bnf-code-information-current-year"               
#> [5] "bnf-code-information-historic"                   
#> [6] "bnf-code-information-monthly-changes"

When you do not already know the id, search for one with nhsbsa_package_search(). It prints a tidy summary — the match count and a table of the matching datasets; tibble::as_tibble() returns that table to work with:

nhsbsa_package_search(q = "prescribing", rows = 5)
#> <nhsbsa package search> 661 datasets found
#> Showing the first 5; increase `rows` for more.
#> # A tibble: 5 × 5
#>   name                        title organisation num_resources metadata_modified
#>   <chr>                       <chr> <chr>                <int> <chr>            
#> 1 prescriber-details          Pres… community_p…            48 2026-07-29T11:29…
#> 2 foi-03835                   FOI-… freedom-of-…            41 2026-06-16T12:59…
#> 3 english-prescribing-data-e… RETI… community_p…           138 2026-03-06T12:49…
#> 4 hospital-prescribing-dispe… Hosp… community_p…           114 2026-07-22T08:04…
#> 5 english-prescribing-datase… Engl… community_p…            67 2026-07-22T14:12…

Inspect a dataset’s resources

A dataset is a container of resources (files). nhsbsa_list_resources() lists them as a tibble, including the download url of each:

resources <- nhsbsa_list_resources("english-prescribing-data-epd")
nrow(resources)
#> [1] 138
resources |>
  select(name, format, url) |>
  slice_head(n = 6)
#> # A tibble: 6 × 3
#>   name       format url                                                         
#>   <chr>      <chr>  <chr>                                                       
#> 1 EPD_201401 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…
#> 2 EPD_201402 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…
#> 3 EPD_201403 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…
#> 4 EPD_201404 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…
#> 5 EPD_201405 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…
#> 6 EPD_201406 CSV    https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d…

Filter by a pattern matched against the resource name:

nhsbsa_list_resources("english-prescribing-data-epd", pattern = "202401") |>
  select(name, id, last_modified)
#> # A tibble: 1 × 3
#>   name       id                                   last_modified             
#>   <chr>      <chr>                                <chr>                     
#> 1 EPD_202401 fe7c75f9-7ac6-4d03-8941-74f596db4a5a 2024-03-19T09:39:16.950431

nhsbsa_resource_show() returns the full metadata for a single resource (by its id), and prints a summary if you need more detail than the table above:

nhsbsa_resource_show(resources$id[[1]])
#> <nhsbsa resource> "EPD_201401"
#> Format: CSV
#> Size: 6618466913
#> Datastore: FALSE
#> Modified: —
#> URL:
#> https://opendata.nhsbsa.net/dataset/65050ec0-5abd-48ce-989d-defc08ed837e/resource/8ae6b792-2a0c-4f4b-826c-dc6483dc32a7/download/epd_201401.csv

Download a resource file

Identify a single resource — by resource_id, or by a pattern that matches exactly one resource name — and stream its file to disk. You choose the destination directory (it must already exist), and the file is saved there under its own name; here we use a (smaller) resource from the BNF code dataset and save to a temporary directory:

bnf <- nhsbsa_list_resources("bnf-code-information-current-year")
path <- nhsbsa_download_resource(
  "bnf-code-information-current-year",
  resource_id = bnf$id[[1]],
  directory = tempdir()
)
#>  Downloading "BNF_CODE_CURRENT_202503_VERSION_88" to
#>   /tmp/RtmppEvjVl/bnf_code_current_202503_version_88.csv.
basename(path)
#> [1] "bnf_code_current_202503_version_88.csv"

Query rows without downloading the whole file

Not every resource can be queried row by row. The datastore is a separate, queryable copy of the tabular resources (CSVs); non-tabular files such as PDFs can only be downloaded. Where a resource is in the datastore, you can query it directly.

Two things to know about this portal specifically:

  • The datastore identifies a resource by its name (e.g. "EPD_202401", shown in the name column of nhsbsa_list_resources()), not by its id.
  • The datastore_active metadata flag is unreliable here (it is often FALSE even for resources that are queryable), so rather than trusting it, simply try the query — CSV resources are generally queryable by name, and a non-tabular resource returns an error.

Ways to query datastore data

There are two functions, and on this portal they have a clear division of labour: use nhsbsa_datastore_search() to read rows, and nhsbsa_datastore_search_sql() to filter or aggregate them.

nhsbsa_datastore_search() returns rows from a resource. You can choose and order columns with fields, sort with sort, and page with limit/offset. Field names are case-sensitive and must match the resource’s columns exactly (for the EPD they are upper case, e.g. PCO_CODE):

nhsbsa_datastore_search(
  resource_id = "EPD_202401",
  fields = c("PCO_CODE", "BNF_CHEMICAL_SUBSTANCE", "ITEMS"),
  sort = "ITEMS desc",
  limit = 5
)
#> Warning: ! Retrieved 5 of 18080573 matching rows; 18080568 not returned.
#>  Fetch the next page with `offset = 5` (reusing your other arguments),
#>   increasing `offset` until all rows are retrieved.
#>  Raising `limit` returns more rows per request, up to the server-side maximum.
#> # A tibble: 5 × 3
#>   PCO_CODE BNF_CHEMICAL_SUBSTANCE ITEMS
#>   <chr>    <chr>                  <int>
#> 1 11J00    1404000H0               3584
#> 2 06H00    0212000B0               3571
#> 3 02Y00    0212000B0               3469
#> 4 12F00    1404000H0               3160
#> 5 11M00    1404000H0               3038

The datastore returns at most one page of rows per request. When more rows exist than were returned, nhsbsa_datastore_search() warns you and explains how to page through the rest by increasing offset:

nhsbsa_datastore_search(
  resource_id = "EPD_202401",
  fields = c("PCO_CODE", "ITEMS"),
  limit = 5,
  offset = 5
)
#> Warning: ! Retrieved 5 of 18080573 matching rows; 18080563 not returned.
#>  Fetch the next page with `offset = 10` (reusing your other arguments),
#>   increasing `offset` until all rows are retrieved.
#>  Raising `limit` returns more rows per request, up to the server-side maximum.
#> # A tibble: 5 × 2
#>   PCO_CODE ITEMS
#>   <chr>    <int>
#> 1 -            1
#> 2 -            2
#> 3 -            1
#> 4 -            1
#> 5 -           21

CKAN’s datastore_search also defines filters (exact field matching) and q (full-text search) parameters, and nhsbsa_datastore_search() exposes them for completeness. Be aware that this portal’s datastore does not apply them — they return no matching rows — so to filter by value, use SQL instead.

Filtering and aggregating with SQL

nhsbsa_datastore_search_sql() runs a read-only SQL query, which is the reliable way to filter, compute expressions, aggregate and sort on this portal. The portal requires the resource_id alongside the query, and you reference the same resource name in the FROM clause:

# Filter to one organisation
nhsbsa_datastore_search_sql(
  resource_id = "EPD_202401",
  sql = "SELECT PCO_CODE, BNF_CHEMICAL_SUBSTANCE, ITEMS
         FROM `EPD_202401`
         WHERE PCO_CODE = 'W2U3Z'
         LIMIT 5"
)
#> # A tibble: 5 × 3
#>   PCO_CODE BNF_CHEMICAL_SUBSTANCE ITEMS
#>   <chr>    <chr>                  <int>
#> 1 W2U3Z    0407020AD                  1
#> 2 W2U3Z    2122                       1
#> 3 W2U3Z    020400080                  1
#> 4 W2U3Z    0101021B0                  1
#> 5 W2U3Z    0205052AE                  1
# Aggregate: total items prescribed per organisation
nhsbsa_datastore_search_sql(
  resource_id = "EPD_202401",
  sql = "SELECT PCO_CODE, SUM(ITEMS) AS items
         FROM `EPD_202401`
         GROUP BY PCO_CODE
         ORDER BY items DESC
         LIMIT 5"
)
#> # A tibble: 5 × 2
#>   PCO_CODE   items
#>   <chr>      <int>
#> 1 91Q00    3151968
#> 2 W2U3Z    3129964
#> 3 A3A8R    3124609
#> 4 D9Y0V    2697474
#> 5 15N00    2413359

The SQL string is sent to the API verbatim, so you are responsible for paging (via LIMIT/OFFSET) and for quoting identifiers correctly.