library(sandwich)
library(lmtest)

d <- read.csv("sample-conjoint-profile-data.csv", check.names = FALSE)
d$x1 <- as.integer(d[["Attribute 1"]] == "Level 1")
d$x2 <- as.integer(d[["Attribute 2"]] == "Level 1")
d$x3 <- as.integer(d[["Attribute 3"]] == "Level 1")

model <- lm(chosen ~ x1 + x2 + x3, data = d)
print(coeftest(model, vcov. = vcovCL(model, cluster = d$response_id)))