Next Article in Journal
Predictors of Adolescent Depressive Symptoms
Next Article in Special Issue
Preliminary Examination of an Appropriate Price Calculation Method and Medical Treatment Costs for Foreign Visitors in Japan
Previous Article in Journal
Current Knowledge, Attitudes, and Practice among Health Care Providers in OSCC Awareness: Systematic Review and Meta-Analysis
Previous Article in Special Issue
Developing a Stigma Scale for the Workplace: Focus on an Airline Cabin Crew
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validation of a Satisfaction Questionnaire on Spa Tourism

Department of Economics and Business Administration, University of Malaga, Campus El Ejido, 29071 Malaga, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(9), 4507; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094507
Submission received: 12 April 2021 / Accepted: 21 April 2021 / Published: 23 April 2021
(This article belongs to the Special Issue Tourism/Hospitality and Well-being)

Abstract

:
The authors’ line of research is within the existing methodological debate around the concepts of quality of services, destinations, and quality measurements methods. The authors consider that the most appropriate way to measure quality is to develop instruments according to the destination and context in question, defining the quality of the tourist destination for practical purposes based on the satisfaction experienced by the tourist or the SERVPERF model, weighted and used to measure the quality of sun and beach tourist destinations. The authors of this work propose the knowledge of spa tourism, its quality and its level of satisfaction as a research gap and consider it as a starting point to validate a questionnaire that would allow the measurement and comparison of parameters with other segments already studied and that can also serve as a measuring instrument for tourist segments with similar characteristics, not as well known in the international literature as inland, ecological or nature tourism. Good internal reliability results were obtained in all items and in all dimensions. The factor analysis distributed the weights of the variables in the theoretical model, and construct validity was obtained with an association between the global evaluation by dimension and the general significance. The score of the main questionnaire was statistically significant.

1. Introduction

The evolution of research in medical tourism includes different groups of topics, such as issues related to health, well-being, thermal tourism, and quality of service, as well as topics related to medical treatments and tourism [1]. The basis for the implementation of a quality system in any company is the strategic quality plan. This involves objectively evaluating the current situation and contemplating the client’s vision [2,3]. There is a flood of definitions about the concept of quality, suggesting that it is a broad and multidimensional term with different interpretations. Quality has been defined from perspectives that are, in some cases, complementary (objective quality versus subjective quality; internal quality versus external quality) and, in other cases, antagonistic (static quality versus dynamic quality; absolute quality versus relative quality). Summarizing the literature, six concepts of quality have been identified [4]:
  • Quality as excellence;
  • Quality as conformity with specifications;
  • Quality as uniformity;
  • Quality as fitness for use;
  • Quality as satisfaction of customer expectations;
  • Quality as value creation, understood as the degree of satisfaction of all the key stakeholders of the organization or total quality.
Intuitively, a relationship can be determined between customer opinion and future results. Measuring perceived quality is “listening” to the customer and using it as a tool to organize opinions and determine areas for improvement on which to act. Quality is a multidimensional concept that encompasses many independent attributes (any item that is part of the definition), and the generic grouping of attributes defines the dimension. The analysis of perceived quality is a highly topical subject and an open field of research [4].
In relation to tourism experiences with water, some authors find a positive correlation with the perception of the quality of life, satisfaction and loyalty of people, both towards the experience and towards the destination [5].
Andalucía (Spain) is a region in Southern Spain with 87,268 km2 [6] and with a population of 8 million people [7], and it is well known as one of the main tourist destinations for sun and beach (e.g., Costa del Sol, Marbella or Malaga) and for cultural tourism (e.g., Córdoba, Granada and Alhambra) [8], in which young tourists and foreigners are predominate in summer and retirees constitute the majority in winter [9,10,11], but spas have not been extensively studied [12]. Therefore, the knowledge of spa tourism, its quality and level of satisfaction is considered in the present work as a research gap, and it is proposed as a starting point to validate a questionnaire that would enable measuring and comparing parameters with other segments and similar studies [13] and methodological approaches [13,14,15,16]. It is also interesting that it can serve as a measuring instrument for tourist segments with similar characteristics not so well known in the international literature, such as inland, ecological or nature tourism [17].

2. Literature Review

Many methodologies based on questionnaires establish theoretical satisfaction as a response to quality services, and there is a broad debate on the advantages and disadvantages of the use of each. Likert-type scales are easy to construct [18] and have high reliability and validity [15,19].
The debate stems from the fact that to fully measure the image of the destination, several components must be captured. These include attribute-based imagery; holistic impressions; and functional, psychological, unique, and common characteristics that require a combination of structured and unstructured methodologies to measure the image of the destination as envisioned in the conceptual framework [18].
Sasser et al. [20] consider that expectations result in attributes. Grönroos [21], representative of the Nordic model, states that the client compares the expected service with the service received. Parasuraman et al. [22,23], representatives of the American school, consider that to determine quality, it would be necessary to compare the difference between perceptions and expectations, and they elaborated a scale, SERVQUAL, with five dimensions. Cronin and Taylor [14] establish quality based solely on perceptions; that is, it eliminates the measurement of expectations but maintains the SERVPERF scale created by Parasuraman et al. [22,23]. Rust et al. [24] consider only three dimensions: result, delivery, and service environment. Dabholkar et al. [25] propose a hierarchy of quality in three levels or primary dimensions, and these, in turn, are divided into sub-dimensions. Brady and Cronin [26] consider dimensions, sub-dimensions, and level: reliability, responsiveness, and empathy. Therefore, the lack of consensus means that there are different models to assess quality.
The first aspect that arises when evaluating the quality of services is the measurement of expectation, which has been modeled on two widely used scales, SERVQUAL, which compares expectations and perceptions [27], and SERVPERF based solely on the measurement of perceptions [28]. Some authors doubt its validity due to theoretical difficulties and practical application measuring expectations [28]; others consider that the comparison between expectations and perceptions of quality is the most appropriate way to measure it [27]. Alén [29] makes a comparison of scales for the measurement of perceived quality in thermal establishments in which he concludes that the SERVPERF scale based on perceptions has better psychometric properties than the one based on the subtractive paradigm (perceptions–expectations).
Another point of debate is whether, when evaluating quality, the importance that users attach to each attribute or dimension should be measured. Cronin and Taylor [27] deduced from the results of their research that weighted models are different for studying the quality of services. Teas [30] makes a comparison between weighted and unweighted models and concludes that the latter are more appropriate in statistical terms. Quester et al. [31] make a comparison between weighted and unweighted models and conclude that the predictive power of SERVQUAL and SERVPERF improves with the inclusion of importance scores. On the other hand, among those in favor of including importance scores, there are also disagreements about how they should be measured. Parasuraman et al. [28] believe that scores should only be obtained for each factor [28], while others [31] argue that the importance should be asked individually for each attribute. In addition, there is a debate about whether importance scores should be obtained directly by asking the respondents [22] or if, on the contrary, they must be obtained indirectly through statistical procedures [13] from another type of information collected in the questionnaires.
The diversity of authors who have tried to define the tourist destination construct [4,32,33,34,35,36] indicates that the concept of tourist destination encompasses a great diversity of services offered by both private companies and public administrations, and even infrastructure and natural resources. Therefore, to evaluate the quality of a tourist destination, the methodologies developed in the field of quality can be applied, although taking into account that it is necessary to adapt them to a specific context [13].
Kozak and Rimmington [11] in Mallorca followed the line of measuring quality by developing and validating a destination measurement instrument composed of a multiattribute questionnaire, and according to Robinson [37], the best way to measure quality is by adapting the different instruments to the context. Otero and Otero [13], on Costa del Sol, used a multiattribute questionnaire that they validated using the CFA technique and developed a model to measure the quality of the tourist destination. They found high levels of satisfaction with the different attributes of the destination and other investigations obtained results along the same line [11,38,39,40,41].

3. Materials and Methods

3.1. Questionnaire Reliability: Test–Retest and Internal Reliability

A cross-sectional survey design was carried out, and the validated questionnaire by Otero and Otero [13], which in turn was inspired by the SERVPERF methodological school [28] study, was the basic instrument to obtain the sample information. The cultural adaptation to the specific spa tourism was carried out based on the contribution of a group of experts made up of university professors, professionals from the sector, and regular users of the spas. The result was a questionnaire with slight modifications to the original that serves to measure a complex construct, such as satisfaction in spa tourism or in segments with similar conditions, for which we have provided an operational definition by meeting with experts. Content validity could be obtained by following Streiner and Norman [42].
Under the previous premises, it was considered a priori that the questionnaire was made up of 7 factors or dimensions, which was later confirmed using the statistical technique of Confirmatory Factor Analysis (n = 725).
Once the adaptation was made, the diagnostic test–retest reproducibility test was carried out on 30 users (12 men and 18 women, aged 25–58 years), who did not belong to the main sample, and three study spas among the cases with a valid response with a separation between one and two days. The intraclass correlation coefficient (ICC) or Kappa index was used depending on the type of variable (quantitative or qualitative), valued with the Landis and Koch scale [43]: <0.00: poor (poor); 0.00–0.20: slight (slightly); 0.21–0.40: fair (fair); 0.41–0.60: moderate (moderate); 0.61–0.80; substantial (substantial); 0.81–1.00: almost perfect (almost perfect).
The internal consistency of the 52-item questionnaire was measured with Cronbach’s alpha coefficient excluding the item, which indicates that the questions that measure the same phenomenon are correlated with one another. The global alpha was obtained for each of the dimensions (accommodation, restoration, spa, sports facilities, leisure/culture/shopping, public roads/urban and natural environment, and transport infrastructure and other services) and by item in relation to its dimension.
The literature suggests that the Likert-type and semantic differential scales are easy to construct and manage [18]. The results of the empirical research show that the Likert-type and semantic differential scales have high reliability and validity [19]. The use of the “delighted–terrible scale” has been reported as the possibility of reducing satisfaction responses [15,16]. “Do not know” was also included in the scale of these possibilities for those who might not have an opinion due to no direct experience with the destination attribute. To obtain a weighting factor, the user was asked to grant a scale score of 1 to 100 on how important each factor was to him or her. The third section of the questionnaire, made up of four questions, was designed to determine overall satisfaction with the destination with the same seven-point scale, the intention to recommend it with three categories, and the intention to return with five, in which the user was only given asked by perceptions.

3.2. Data Collection Procedures

The sample was designed based on information gathered from external sources, in particular, a telephone and/or email survey of individuals responsible for managing spas. The relevant data were used to calculate the number of users per year per spa, as well as their distribution by age group and length of stay. The overall population size was estimated at 53,231 users per year with an age equal to or greater than 15 years old. These individuals were thus defined as this study’s primary subjects.
Stratified sampling was used, with each stratum coinciding with one of the participating spas. The sample size was estimated using Equation (1):
n c = D e f f ×   n
in which n c is the size adjusted to ensure complete sampling (i.e., strata) and D e f f is the design’s effect or the relationship between the variance under stratified sampling and under simple random sampling. This variance was arbitrarily estimated a priori as 1.5, which was shown to be an overestimate. Any proportion calculated in this study had to meet the requirements of an accuracy of 5% and a 95% confidence level, with the worst possibility considered to be p = 0.5. This implied that n = 384 and nc = 576, which was set as the minimum.
The sample size and/or stratum was the result of a proportional fixation (i.e., the largest spas received more than one visit during the survey). The final result was 725 valid questionnaires, which coincides with some authors’ suggestion that sample size [44] be based on 10 subjects per analyzed variable. The sample size’s adequacy was thus confirmed, since 500 subjects would be an acceptable sample size and 1000 or more subjects would be an excellent sample size.

3.3. Factor Analysis, Calculation of Dimensions, and Validity of the Construct

In the factorial model of Otero and Otero [13] applied to sun and beach tourism, seven dimensions were obtained (accommodation, restaurants, beaches, sports facilities, leisure/culture/shopping, public roads/urban and natural environment, and transport infrastructures). In the present study, CFA was carried out to determine the degree to which the group of factors identified for sun and beach tourism is capable of representing the data from the spa tourism matrix.
To obtain a model with parsimony—that is, to maximize the amount of variance of the variables that can be explained by the minimum underlying factors or components—the Cattell slope test was performed, which consists of graphically representing the extracted factors (placed on the abscissa axis) against their eigenvalues (placed on the ordinate axis) to establish an inflection point on the graph. To obtain a more easily interpretable grouping of variables, the orthogonal rotation of factors was performed using the Varimax method, as it is the most frequently used [3,45]. To calculate the weight of each dimension, weights derived from the importance that each user gave to each of them were constructed, so that in each user, it would add up to 100 with all the weights in such a way that the weighting given by the i-th user was defined as follows:
w i j = I i j J = 1 7 I i j × 100  
where I i j is the importance given by use i to factor j.
With this definition, the sum of the weights for each individual should be equal to 100.
J = 1 7 w i j = 100
The weight variable multiplies the value of each dimension to estimate the total score. The percentage distribution was not calculated for the weights, since it is the means that are relevant, taken from Otero and Otero [13]. The score of the items of global assessment by dimension and general satisfaction in users (n = 725) corresponded to the single question by dimension and global satisfaction questions and are different from the score by dimension estimated from the main questionnaire of 52 items.
The internal consistency of the questionnaire was complemented with questions of its future intention. A priori, it was to be expected that the response of intention to recommend it and intention to return would be related to the level of satisfaction, as observed in the satisfaction study of sun and beach tourism [10]. To determine the association between global assessment by dimension and general satisfaction and the score of the main questionnaire, the Pearson correlation coefficient was calculated.
To calculate the association between the variables of future intention and the scores estimated in the main questionnaire (52 items and 725 users), the comparison of means was used. The aim was to check whether the difference between the two means was due to satisfaction influencing the intention to recommend it and the intention to return, or if the observed differences could simply be due to random variability. For this, the Student’s t test was applied for independent samples, with weighted samples. The null hypothesis was that both means are equal in the population. The alternative hypothesis maintained that the means were different and that both had different effects.
In all statistical analyses, the weighting of the observations was taken into account. For descriptive statistics, the SPSS Windows 15.0 program (SPSS Inc., Chicago, IL, USA) was used, and for analytical statistics (calculation of p values, as well as confidence intervals), the SUDAAN 7.0 program (RTI, RTP, Durham, NC, USA), indicating the STRWR sample design (stratified with replacement, the conglomerates being the strata). It was not corrected for a finite population, since the sample fraction (ratio between sample size and population) was much less than 10%. All the statistical procedures used, both those discussed in terms of reliability and validity and those carried out to evaluate the association of sociodemographic variables and characteristics of the type of visit by tourists to spas, are expressed in detail at the bottom of each results table.

4. Results

4.1. Sociodemographic Description of the Users and the Visit to the Spa

The average age of spa tourists in Andalusia is 56 years old, and the tourists are predominantly female. Fifty-three percent of users are retired, and 34% are employed. Moreover, 29.8% have an average monthly income of EUR 500, 33.1% have an average monthly income between EUR 501 and 1000, and 26.6% are in the interval between EUR 1001 and 1500 (Table 1). The origin of these users is mainly local, with 82.9% from the Andalusian community, with no foreigners in the sample (Table 2). In addition, 51.3% of spa tourists in Andalusia go as a couple and 19.7% as a family. In the accommodation modality, 49.6% choose a three-star hotel. The highest percentage of overnight stays is among the IMSERSO (48.4%), followed by weekend and long-weekend tourists (28.1%), with an average stay of one week. Moreover, 46.8% of users who go to a spa repeat their visit, contracting the trip without intermediaries (99.1%) (Table 1).
At the population level, 92.4% express the intention of recommending the spa compared to 2.0% who have no intention of recommending it. Additionally, 63.6% of the population express the intention to return next year compared to 1.6% who have no intention of ever returning (Table 2).

4.2. Validation of the Satisfaction Questionnaire in Spa Tourism

4.2.1. Test–Retest

For the diagnostic agreement of the 52 items of the main questionnaire, a test–retest was applied (with a separation between 1 and 2 days), analyzed with the Intraclass Correlation Coefficient (ICC), valued with the Landis and Koch scale (39): <0.0: poor; 0.00–0.20: slight; 0.21–0.40: fair; 0.41–0.60: moderate; 0.61–0.80: substantial; 0.81–1.00: almost perfect.
To calculate the diagnostic concordance of the 52 items of the main questionnaire, the test–retest (with a separation between one and two days), the Intraclass Correlation Coefficient (ICC) was used, assessed with the Landis and Koch scale [43]: <0.0: poor; 0.00–0.20: slight; 0.21–0.40: fair; 0.41–0.60: moderate; 0.61–0.80: substantial; 0.81–1.00: almost perfect.
Diagnostic concordance was found with scores higher than 0.65 for almost all the items of the accommodation dimension. In the restoration dimension, values higher than those of accommodation were found, mostly above 0.67; when analyzing the spa dimension, we found values higher than 0.51; in the sports facilities dimension, these values were above 0.67; regarding leisure/culture/shopping, values higher than 0.61 were obtained; for the dimension of public roads, and urban and natural environments, the values were in almost all respects above 0.36; and finally, for the dimension of transport, infrastructures and other services, the CCI value was higher in all variables than 0.61 (Table 3).
In the general assessment questions (global assessment by dimension, importance, general satisfaction, satisfaction in quality/price, intention to recommend it, and intention to return), an ICC of greater than 0.61 was obtained in all the items of the quantitative variables and a kappa index higher than 0.70 for the qualitative ones (Table 4).
Table 4. Reliability of the questionnaire. General valuation questions. Test–retest with a separation between one and two days n = 30.
Table 4. Reliability of the questionnaire. General valuation questions. Test–retest with a separation between one and two days n = 30.
QuestionInitialRepetition (24–48 h.)Concordance a
Global assessment by dimension
 Accommodation5.87 ± 1.255.77 ± 1.10ICC = 0.75
 Catering (meals)5.83 ± 0.995.53 ± 1.01ICC = 0.46
 Spa5.50 ± 0.905.37 ± 1.03ICC = 0.50
 Sports facilities4.53 ± 1.964.43 ± 1.77ICC = 0.83
 Leisure, culture, and shopping4.37 ± 1.404.43 ± 1.65ICC = 0.66
 Urban and natural environment5.27 ± 1.235.07 ± 1.34ICC = 0.70
 Infraest. transport and other serv.4.43 ± 1.434.93 ± 1.36ICC = 0.67
Importance of the factor (mean ± of b)
 Accommodation85 ± 1887 ± 13ICC = 0.62
 Catering (meals)80 ± 2280 ± 14ICC = 0.61
 Spa80 ± 2476 ± 18ICC = 0.35
 Sports facilities45 ± 2450 ± 22ICC = 0.68
 Leisure, culture, and shopping59 ± 2161 ± 21ICC = 0.73
 Urban and natural environment76 ± 2277 ± 20ICC = 0.85
 Infraest. transport and other serv.58 ± 2360 ± 23ICC = 0.77
General stay satisfaction (mean ± de)5.97 ± 1.035.90 ± 0.99ICC = 0.81
Satisfaction in quality/price (mean ± de)5.83 ± 1.095.77 ± 0.94ICC = 0.71
Intention to recommend it (%) kappa c = 1
 Yes93.393.3
 No3.36.7
 Does not know3.30.0
Intention to return (%) kappa = 0.70
 Next year13.313.3
 This year13.310.0
 Ever43.360.0
 Does not know26.713.3
 Never3.33.3
a: Measured by ICC or kappa, depending on the type of variable. Both ICC and kappa are interpreted with the Landis and Koch scale (see footnote d in Table 5). b: Arithmetic mean ± SD unweighted. c: The kappa was calculated after eliminating the case (3.3%) that did not answer in the initial assessment.
Table 5. Internal reliability of the questionnaire (Cronbach’s alpha) (n = 725) a.
Table 5. Internal reliability of the questionnaire (Cronbach’s alpha) (n = 725) a.
By Dimension By Item (in Relation to Its Size)
DimensionNo. ItemsnbGlobal AlphaItemCorrected Item- Total CorrelationAlpha without Item
Accommodation86360.87101-Cleaning, hygiene0.7560.842
02-Security0.6880.850
03-Environmental quality0.5560.862
04-Room0.6890.848
05-Aesthetics0.7300.843
06-Ease for children and disabled0.3890.892
07-Personal professionalism0.6880.850
08-Value for money0.6670.851
Catering67100.91909-Variety0.8230.897
10-Cleaning0.7760.904
11-Environment0.6640.918
12-Food quality0.8400.894
13-Personal professionalism0.6960.914
14-Quality/price ratio0.8290.896
Spa96820.83315-Medicinal mineral waters0.5770.817
16-Cleaning0.6690.808
17-Security0.5410.820
18-Equipment0.6010.809
19-Environment0.5790.812
20-Medical service0.5130.819
21-Handicapped access0.4250.851
22-Personal professionalism0.6450.805
23-Value for money0.6450.804
Sports facilities65200.90724-For children0.5330.917
25-Various0.6940.898
26-Soft (paddle, etc.)0.7950.883
27-Tennis0.8350.877
28-Pool0.7300.893
29-Others0.8790.869
Leisure/culture/shopping75620.83030-Serv. complem. (parking)0.4400.828
31-Pub/discos/cafes0.7800.775
32-Excursions0.6000.803
33-Culture/show0.6900.789
34-Shops0.7860.773
35-Animation of streets0.4910.820
36- Tourist office0.3570.853
Public roads/urban and natural environment97130.86937-Cleaning0.5850.857
38-Absence of beggars0.6230.857
39-Feeling of security0.7050.847
40-Little noise0.6270.853
41-Traffic0.6830.847
42-Urban furniture0.5600.863
43-Attitude of the people0.5270.861
44-Urban landscape0.6680.848
45-Natural environment0.5750.858
Transport infrastructures and other services75420.80346-Taxis0.3090.816
47-Buses0.4840.787
48-Roads0.4490.793
49-Line buses0.5710.771
50-Post office0.7100.742
51-Public telephones0.5790.770
52-Banking services0.6600.753
a: Note that it is a weighted sample. b: Number (unweighted) of users out of the total of 725 who responded to all items in the dimension.

4.2.2. Internal Reliability

The internal consistency of the questionnaire measured with Cronbach’s alpha coefficient indicates a high correlation between the items. The global values for accommodation were at 0.871; restoration: 0.919; spa: 0.833; sports facilities: 0.907; leisure/culture/shopping: 0.830; public roads/urban and natural environment: 0.869; and transport infrastructure and other services: 0.803 (Table 6).
By item (in relation to its dimension), we found a corrected item–total correlation within the accommodation dimension for “cleanliness and hygiene” of 0.756 and for “aesthetics” of 0.730 as the highest values. The lowest value was obtained by “facility for children and the disabled” at 0.389. For the restoration dimension, the value was 0.840 for “food quality” and 0.829 for “value for money” as high values, and 0.664 for “environment” as the lowest value. For the spa dimension, these values were between 0.669 for “cleanliness” and 0.425 for “disabled access”. In the sports facilities dimension, we obtained 0.879 for “other sports facilities” and the lowest value of 0.533 for “children”. In the case of the leisure/culture/shopping dimension, the highest value was 0.786 for “shops”, and the lowest value was 0.357 for “tourist office”. In public roads/urban and natural environment, we obtained 0.705 for “sense of security” and 0.527 for “people’s attitude”. Finally, in the case of the dimension of transport infrastructures and other services, the highest value was 0.710 for “post office”, and the lowest values were 0.309 for “taxis” and 0.449 for “roads” (Table 5).

4.2.3. Factor Analysis and Calculation of Dimensions

The KMO test, which compares the Pearson correlation coefficients between each pair of variables with their respective partial correlation coefficients and whose values oscillate between 0 and 1, obtained an excellent value of adequacy of the weighted sample, KMO = 0.923 (Table 6). The Bartlett sphericity test, whose null hypothesis is that there is no correlation between variables, obtained a p value <0.001, which rejects the null hypothesis of the absence of correlation between variables, which indicates the adequacy of the CFA (Table 7).
The Cattel slope test (Figure 1) indicates a number of factors to be extracted from eight factors, coinciding with the theoretical model (38).
After the Varimax orthogonal rotation, a grouping of variables was obtained that distributed the weights over the dimensions indicated in the theoretical model, with slight imbalances in the variables that distributed a similar weight in two different factors, such as the variables ease for children and the disabled, medical service, complementary services, people’s attitudes, and natural environment. The infrastructure and other services dimension was divided into two, leaving, on the one hand, roads and other services and, on the other, the variables referring to transport, such as taxis, urban buses, and line buses (Table 7). The variance explained with rotated factors (Varimax rotation) is shown in Table 8.

4.2.4. Validity of the Construct

The association (Pearson correlation) between global assessment by dimension, general satisfaction and score of the main questionnaire (n = 725) of the weighted sample in all cases was statistically significant (p < 0.001). Pearson’s linear correlation coefficient r showed positive values close to 1, 0.91 for the restoration dimension and 0.69 for transport infrastructure and other services, which showed a close and direct association between them (Table 9).
Regarding the association (comparison of means) between the variables of future intention and the scores estimated with the main questionnaire of 52 items (n = 725) of the weighted sample, they are significant (p < 0.001) with respect to greater satisfaction (Table 10).

5. Discussion and Conclusions

The questionnaire had three structured sections for the purpose of measuring the satisfaction of spa tourists with general identification questions, specific questions of satisfaction, their probability of returning, and their probability of recommending the destination. The second section, made up of 52 items, was based on a seven-point scale from delighted to terrible, measuring seven factors, namely, accommodation, restaurants, spa, sports facilities, leisure/culture/shopping, public roads/urban and natural environment, and infrastructures and other services.
Our research coincides with other studies that have identified specific attributes of the destination at a lower cost [11,13].
The reference model was taken into account and contrasted with the suggestions of the experts for adaptation to the context, finding a good general fit with the theoretical model and some particularities, such as the variables of ease for children and the disabled, medical service, complementary services, the attitude of the people, and the natural environment, which distributed their weight in two different factors and which the authors define as complex variables [46]. The dimension “infrastructures and other services” was divided into two, leaving, on the one hand, roads and other services and, on the other, the variables related to transport, such as taxis, urban buses, and scheduled buses. Otero and Otero [13] justify the categorization of the items based on the advantage of considering that the dimensions within several factors obtain more useful information to carry out quality policies, since it allows the identification of those responsible. On the other hand, it would be logical to disaggregate transport infrastructures and other services based on the practical knowledge acquired in this work. The CFA for five dimensions explains 48% of the total variance compared to the study by Otero and Otero [13], where it only explains 36%. However, in the analysis of the results with eight components, it explains 65.1% of the total variance.
All this questions the validity of the construct. However, the results of the global assessment obtained by dimension in the weighted sample corresponding to the single question per dimension of the questionnaire and different from the one estimated from the main questionnaire of 52 items, as well as the general questions regarding general satisfaction, stood out with higher average values. The study found statistical significance between the highest satisfaction values and the variables of future intention. Other investigations obtained similar results [11,13,32,39,40,41]. The validated questionnaire of satisfaction on sun and beach tourism in Costa del Sol is adapted to the spa sector, adjusting to the seven initial dimensions, except for transport infrastructures, which is divided into two, on the one hand, roads and public telephones, and on the other hand, taxis and buses.
The findings of this research show a solid instrument that can be used by different institutions to replicate and use as a barometer of spa tourism and culturally adapted to be used for similar segments. With this instrument, one can obtain both the level of satisfaction and the factors associated with it, and even segment the market.
The analysis is limited only to spas located in Southern Spain, so the results can be extrapolated to any spa in Spain or any country with similar climatic conditions.

Author Contributions

Conceptualization, R.A.-A.; methodology, R.A.-A., G.G. and C.A.-A.; software, R.A.-A., G.G. and C.A.-A.; validation, R.A.-A., G.G. and C.A.-A.; formal analysis, R.A.-A., G.G. and C.A.-A.; investigation, R.A.-A., G.G. and C.A.-A.; resources, R.A.-A., G.G. and C.A.-A.; data curation, R.A.-A., G.G. and C.A.-A.; writing—original draft preparation, R.A.-A., G.G. and C.A.-A.; writing—review and editing, R.A.-A., G.G. and C.A.-A.; visualization, R.A.-A., G.G. and C.A.-A.; supervision, R.A.-A., G.G. and C.A.-A.; project administration, R.A.-A., G.G. and C.A.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank the professionals and experts of the Andalusian spa sector for their participation in this research project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De la Hoz-Correa, A.; Muñoz-Leiva, F.; Bakucz, M. Past themes and future trends in medical tourism research: A co-word analysis. Tour. Manag. 2018, 65, 200–211. [Google Scholar] [CrossRef]
  2. Anaya-Aguilar, R.; Anaya-Aguilar, C.; Bravo-Pérez, M. Balneario: La definición de una tradición milenaria en Andalucía. Rev. Derecho Empres. Soc. 2017, 11, 237–246. [Google Scholar]
  3. Anaya Aguilar, R.M. Diagnóstico y Tendencias del Turismo de Balnearios en Andalucía.Universidad de Málaga, Servicio de Publicaciones. 2011. Available online: https://riuma.uma.es/xmlui/handle/10630/5077 (accessed on 30 December 2020).
  4. Camisón, C.; Cruz, S.; González, T. Gestión de la Calidad: Conceptos, Enfoques, Modelos y Sistemas; Pearson Prentice Hall: Madrid, Spain, 2006; ISBN 978-84-205-4262-1. [Google Scholar]
  5. Campón-Cerro, A.M.; Di-Clemente, E.; Hernández-Mogollón, J.M.; Folgado-Fernández, J.A. Healthy Water-Based Tourism Experiences: Their Contribution to Quality of Life, Satisfaction and Loyalty. Int. J. Environ. Res. Public Health 2020, 17, 1961. Available online: https://0-www-mdpi-com.brum.beds.ac.uk/1660-4601/17/6/1961 (accessed on 2 April 2021). [CrossRef] [PubMed] [Green Version]
  6. Consejería de Turismo y Deporte de la Junta de Andalucía. Inicio—Web Oficial de Turismo de Andalucía. 2020. Available online: https://www.andalucia.org/es/conoce-andalucia/situacion-geografica (accessed on 8 February 2021).
  7. Instituto Estadistica y Cartografia de Andalucia (IECA). Junta de Andalucía. Padrón Municipal de Habitantes. Cifras Oficiales de Población Municipal | Instituto de Estadística y Cartografía de Andalucía. 2021. Available online: https://www.juntadeandalucia.es/institutodeestadisticaycartografia/padron/index.htm (accessed on 3 March 2021).
  8. Instituto Estadistica y Cartografia de Andalucia (IECA). Junta de Andalucia. Encuesta de Coyuntura Turística de Andalucía. Cuarto Trimestre 2019 y año 2020 | Instituto de Estadística y Cartografía de Andalucía. 2021. Available online: https://www.juntadeandalucia.es/institutodeestadisticaycartografia/turismo/notaprensa.htm (accessed on 3 March 2021).
  9. Consejería de Turismo y Comercio de Andalucía. Turismo de Salud y Bienestar. 2015, p. 1. Available online: https://www.juntadeandalucia.es/turismoydeporte/opencms/areas/servicios/centro-documentacion/publicaciones/turismo/Turismo-de-salud-y-bienestar-en-Andalucia-00001/ (accessed on 31 January 2021).
  10. Otero Cordones, C. Evaluación de la Calidad de Destinos Turísticos de Soly Playa: Una Aplicación a la Costa del Sol. Universidad de Málaga. 2003. Available online: http://www.biblioteca.uma.es/bbldoc/tesisuma/16699336.pdf (accessed on 2 April 2021).
  11. Kozak, M.; Rimmington, M. Tourist Satisfaction with Mallorca, Spain, as an Off-Season Holiday Destination. J. Travel Res. 2000, 38, 260–269. [Google Scholar] [CrossRef]
  12. Cong, L.C. A formative model of the relationship between destination quality, tourist satisfaction and intentional loyalty: An empirical test in Vietnam. J. Hosp. Tour. Manag. 2016, 26, 50–62. Available online: https://linkinghub.elsevier.com/retrieve/pii/S1447677015300528 (accessed on 2 April 2021). [CrossRef]
  13. Otero Moreno, J.; Otero Cordones, C. Evaluación de la calidad de destinos turísticos: El caso de la Costa del Sol. Pap. Econ. Esp. 2004, 236–249. Available online: https://dialnet.unirioja.es/servlet/articulo?codigo=1068213 (accessed on 2 April 2021).
  14. Cronin, J.J.; Taylor, S.A. Servperf versus Servqual: Reconciling Performance-Based and Perceptions-Minus-Expectations Measurement of Service Quality. J. Mark. 1994, 58, 125–131. [Google Scholar] [CrossRef]
  15. Maddox, R.N. Measuring Satisfiaction with Tourism. J. Travel Res. 1985, 23, 2–5. [Google Scholar] [CrossRef]
  16. Westbrook, R.A. A Rating Scale for Measuring Product/ Service Satisfaction. J. Mark. 1980, 44, 68–72. [Google Scholar] [CrossRef]
  17. Consejería de Turismo y Deporte de la Junta de Andalucía. www.andalucia.org—Web Oficial de Turismo de Andalucía. 2021. Available online: https://www.andalucia.org/es/segmentos-turisticos (accessed on 30 March 2021).
  18. Echtner, C.M.; Ritchie, J.R.B. The Measurement of Destination Image: An Empirical Assessment. J. Travel Res. 1993, 31, 3–13. [Google Scholar] [CrossRef]
  19. Westbrook, R.A.; Oliver, R.L. The Dimensionality of Consumption Emotion Patterns and Consumer Satisfaction. J. Consum. Res. 1991, 18, 84. [Google Scholar] [CrossRef]
  20. Sasser, W.E.; Olsen, R.P.; Wyckoff, D.D.; Harvard University. Graduate School of Business Administration. Management of Service Operations: Text, Cases, and Readings; Allyn and Bacon: Boston, MA, USA, 1978; Available online: https://www.mendeley.com/catalogue/c8d0fa2c-05aa-395d-8001-261f6f9575cc/?utm_source=desktop (accessed on 24 December 2020).
  21. Gronroos, C.; Grönroos, C.; Gronroos, C. A Service Quality Model and its Marketing Implications. Eur. J. Mark. 1984, 18, 36–44. [Google Scholar] [CrossRef]
  22. Parasuraman, A.; Zeithaml, V.A.; Berry, L.L. A Conceptual Model of Service Quality and Its Implications for Future Research. J. Mark. 1985, 49, 41. Available online: https://0-www-jstor-org.brum.beds.ac.uk/stable/1251430?origin=crossref (accessed on 2 April 2021). [CrossRef]
  23. Parasuraman, A.; Zeithaml, V.A.; Berry, L.L. Servqual: A Multiple-Item Scale for Measuring Consumer Perceptions of service quality. J. Retail. 1988, 64, 12–40. [Google Scholar]
  24. Rust, R.T.; Zahorik, A.J.; Keiningham, T.L. Return on Quality (ROQ): Making Service Quality Financially Accountable. J. Mark. 1995, 59, 58–70. [Google Scholar] [CrossRef]
  25. Dabholkar, P.A.; Thorpe, D.I.; Rentz, J.O. A measure of service quality for retail stores: Scale development and validation. J. Acad. Mark. Sci. 1996, 24, 3–16. Available online: http://0-link-springer-com.brum.beds.ac.uk/10.1007/BF02893933 (accessed on 2 April 2021). [CrossRef]
  26. Brady, M.K.; Cronin, J.J. Some new thoughts on conceptualizing perceived service quality: A hierarchical approach. J. Mark. 2001, 65, 34–49. [Google Scholar] [CrossRef] [Green Version]
  27. Cronin, J.J.; Taylor, S.A. Measuring Service Quality: A Reexamination and Extension. J. Mark. 1992, 56, 55. Available online: https://0-www-jstor-org.brum.beds.ac.uk/stable/1252296?origin=crossref (accessed on 2 April 2021). [CrossRef]
  28. Parasuraman, A.; Berry, L.L.; Zeithaml, V.A. Perceived service quality as a customer-based performance measure: An empirical examination of organizational barriers using an extended service quality model. Hum. Resour. Manag. 1991, 30, 335–364. [Google Scholar] [CrossRef]
  29. Alén González, M.E. Comparison of scales for the measurement of perceived quality in thermal spas. Rev. Galega Econ. Publ. Interdiscip. Fac. Cienc. Econ. Empres. 2006, 15, 25–44. [Google Scholar]
  30. Teas, R.K. Expectations as a Comparison Standard in Measuring Service Quality: An Assessment of a Reassessment. J. Mark. 1994, 58, 132–139. [Google Scholar] [CrossRef]
  31. Quester, P.G.; Romaniuk, S.; Wilkinson, J.W. A Test of Four Service Quality Measurement Scales: The Case of the Australian Advertising Industry. In Proceedings of the 1995 World Marketing Congress; Springer: Cham, Switzerland, 2015; pp. 384–391. [Google Scholar]
  32. Bigné, J.E.; Andreu, L.; Gnoth, J. The theme park experience: An analysis of pleasure, arousal and satisfaction. Tour. Manag. 2005, 26, 833–844. [Google Scholar] [CrossRef]
  33. Hu, Y.; Ritchie, J.R.B. Measuring Destination Attractiveness: A Contextual Approach. J. Travel Res. 1993. [Google Scholar] [CrossRef]
  34. Kim, S.; Yoon, Y. The Hierarchical Effects of Affective and Cognitive Components on Tourism Destination Image. J. Travel Tour. Mark. 2003, 14, 1–22. [Google Scholar] [CrossRef]
  35. Rodríguez-Díaz, M.; Espino-Rodríguez, T.F. A Model of Strategic Evaluation of a Tourism Destination Based on Internal and Relational Capabilities. J. Travel Res. 2008, 46, 368–380. [Google Scholar] [CrossRef]
  36. Giddens, A. Time, space and regionalisation. In Social Relations and Spatial Structures; Palgrave: London, UK, 1985; pp. 265–295. [Google Scholar]
  37. Robinson, S. Measuring service quality: Current thinking and future requirements. Mark. Intell. Plan. 1999, 17, 21–32. [Google Scholar] [CrossRef]
  38. Bigné, J.E.; Sánchez, M.I.; Sánchez, J. Tourism image, evaluation variables and after purchase behaviour: Inter-relationship. Tour. Manag. 2001, 22, 607–616. [Google Scholar] [CrossRef]
  39. Castro, C.B.; Armario, E.M.; Ruiz, D.M. The influence of market heterogeneity on the relationship between a destination’s image and tourist’s future behaviour. Tour. Manag. 2007, 28, 175–187. Available online: https://0-www-sciencedirect-com.brum.beds.ac.uk/science/article/pii/S0261517705002001?via%3Dihub (accessed on 2 April 2021). [CrossRef]
  40. Chen, C.-F.; Tsai, D. How destination image and evaluative factors affect behavioral intentions? Tour. Manag. 2007, 28, 1115–1122. Available online: https://0-linkinghub-elsevier-com.brum.beds.ac.uk/retrieve/pii/S0261517706001397 (accessed on 2 April 2021). [CrossRef]
  41. Lee, T.-H.H. A structural model for examining how destination image and interpretation services affect future visitation behavior: A case study of Taiwan’s Taomi eco-village. J. Sustain. Tour. 2009, 17, 727–745. [Google Scholar] [CrossRef]
  42. Streiner, D.L.; Norman, G.R. Health Measurement Scales: A practical Guide to Their Development and Use; Oxford University Press: New York, NY, USA, 2008; ISBN 9780191724015. [Google Scholar]
  43. Landis, J.R.; Koch, G.G. The Measurement of Observer Agreement for Categorical Data. Biometrics 1977, 33, 159. [Google Scholar] [CrossRef] [Green Version]
  44. Martínez-González, M.A.; Sánchez-Villegas, A.; López del Burgo, C. Introducción a los modelos multivariables. In Bioestadística Amigable; Martínez-González, M.A., Sánchez-Villegas, A., Faulín-Fajardo, F.J., Eds.; Ediciones Díaz de Santos, S.A.: Madrid, Spain, 2006; pp. 685–750. ISBN 84-7978-791-0. [Google Scholar]
  45. Hair, J.F.; Anderson, R.E.; Tatham, R.L.; Black, W.C. Análisis Multivariante; Prentice Hall: Madrid, Spain, 1999; ISBN 9788578110796. [Google Scholar]
  46. Norman, G.R.; Streiner, D.L. Bioestadística. 1996. Available online: https://0-scholar-google-com.brum.beds.ac.uk/scholar?hl=es&as_sdt=0%2C5&q=Norman%2C+G.R.%2C+Streiner%2C+DL.%281996%29.+Bioestadística.+Madrid%3A+Mosby%2FDoyma+Libros%2C+S.A.&btnG= (accessed on 27 December 2020).
Figure 1. Cattel slope.
Figure 1. Cattel slope.
Ijerph 18 04507 g001
Table 1. Sociodemographic distribution of users (n = 725).
Table 1. Sociodemographic distribution of users (n = 725).
SamplePopulation
Variablen (%)% ± se a
Age (years)
 15–35137 (18.9)18.1 ± 1.6
 36–49116 (16.0)14.6 ± 1.5
 50–64149 (20.6)21.2 ± 1.8
 65–88323 (44.6)46.1 ± 2.1
 Mean ± SD/se b55.2 ± 18.356.0 ± 0.8
Sex
 Male286 (39.4)38.2 ± 2.1
 Woman439 (60.6)61.8 ± 2.1
Occupation
 Occupied268 (37.0)34.0 ± 2.0
 Housewife45 (6.2)5.8 ± 1.0
 Student (n = 48) or unemployed (n = 3)51 (7.0)7.0 ± 1.1
 Retired361 (49.8)53.1 ± 2.1
Monthly income (euros)
 <500207 (28.6)29.8 ± 2.0
 501–1000229 (31.6)33.1 ± 2.0
 1001–1500194 (26.8)26.6 ± 1.9
 501–200077 (10.6)8.6 ± 1.1
 2001–250011 (1.5)1.0 ± 0.3
 2501–30004 (0.6)0.7 ± 0.4
 >30003 (0.4)0.2 ± 0.1
Region of origin c
 Total Andalusia596 (82.2)82.9 ± 1.6
 Canary Islands2 (0.3)0.1 ± 0.1
 Cantabria1 (0.1)0.0 ± 0.0
 Castilla Leon2 (0.3)0.1 ± 0.1
 Castilla Mancha13 (1.8)2.6 ± 0.8
 Catalonia22 (3.0)2.3 ± 0.5
 Estremadura3 (0.4)0.2 ± 0.1
 The Rioja2 (0.3)0.2 ± 0.1
 Madrid55 (7.6)7.9 ± 1.1
 Murcia4 (0.6)1.1 ± 0.6
 Basque Country2 (0.3)0.2 ± 0.2
 Valencia16 (2.2)1.7 ± 0.5
 Ceuta/Melilla7 (1.0)0.5 ± 0.2
Andalusian Provinces (on national total)
 Almeria73 (10.1)8.3 ± 0.8
 Cadiz39 (5.4)4.5 ± 0.8
 Cordova25 (3.4)2.7 ± 0.6
 Granada127 (17.5)18.3 ± 1.6
 Huelva5 (0.7)0.9 ± 0.5
 Jaen104 (14.3)20.4 ± 1.6
 Malaga178 (24.6)21.9 ± 1.7
 Seville45 (6.2)5.7 ± 1.0
a: Population percentage ± standard error, taking into account the weight of the sample and stratified sampling (Spas) with SUDAAN. b The standard deviation (SD) was calculated in the case of the sampling distribution and the standard error (se) in the population distribution. c In the sample, there were no users of foreign origin or from the regions of Aragon, Asturias, Balearic Islands, Galicia, or Navarra.
Table 2. Distribution of the variables of future intention of users (n = 725) a.
Table 2. Distribution of the variables of future intention of users (n = 725) a.
SamplePopulation
Variablen (%)% ± se b
Intent to recommend it
 Yes665 (91.7)92.4 ± 1.1
 No17 (2.3)2.0 ± 0.6
 Does not know43 (5.9)5.6 ± 1.0
Intention to return
 Next year470 (64.8)63.6 ± 2.0
 This year96 (13.2)11.0 ± 1.1
 Ever93 (12.8)16.0 ± 1.7
 Does not know54 (7.4)7.7 ± 1.2
 Never12 (1.7)1.6 ± 0.5
a There is no casa in which they do not answer these questions. b Population percentage ± standard error, taking into account the weight of the sample and stratified sampling (Spas) with SUDAAN.
Table 3. Reliability of the 52 items of the main questionnaire. Test–retest (1–2 days apart) (n = 30 a).
Table 3. Reliability of the 52 items of the main questionnaire. Test–retest (1–2 days apart) (n = 30 a).
Dimension b: Item–Item No.InitialRepetition (24–48 h.)
Mean ± SD c
ICC d
ACC: 01-Cleaning, hygiene5.90 ± 1.325.87 ± 1.200.80
ACC: 02-Security5.67 ± 1.215.57 ± 1.300.65
ACC: 03-Environmental quality5.87 ± 1.075.43 ± 1.190.65
ACC: 04-Room5.63 ± 1.165.33 ± 1.180.46
ACC: 05-Aesthetics5.37 ± 1.565.33 ± 1.420.66
ACC: 06-Facility for children and disabled4.93 ± 1.744.87 ± 1.530.84
ACC: 07-Personal professionalism5.57 ± 1.435.57 ± 1.040.74
ACC: 08-Quality/price ratio5.77 ± 1.145.67 ± 1.030.67
CAT: 09-Variety5.10 ± 1.275.17 ± 1.150.68
CAT: 10-Cleaning5.90 ± 0.995.83 ± 0.910.67
CAT: 11-Environment5.23 ± 1.335.17 ± 1.260.68
CAT: 12-Food quality5.63 ± 1.075.43 ± 1.070.56
CAT: 13-Personal professionalism5.90 ± 1.165.57 ± 1.040.71
CAT: 14-Ratio quality/price5.77 ± 1.075.60 ± 1.160.53
SPA: 15- Medicinal mineral waters6.00 ± 1.085.70 ± 1.260.51
SPA: 16-Cleaning6.07 ± 1.055.93 ± 1.110.68
SPA: 17-Security5.57 ± 1.015.83 ± 1.150.43
SPA: 18-Equipment5.50 ± 1.045.20 ± 1.270.67
SPA: 19-Environment5.33 ± 1.305.43 ± 1.070.51
SPA: 20-Medical service4.87 ± 1.634.90 ± 1.400.64
SPA: 21-Handicapped access4.40 ± 1.814.50 ± 1.590.67
SPA: 22-Personal professionalism5.60 ± 1.135.27 ± 1.440.51
SPA: 23-Value for money5.50 ± 0.905.30 ± 1.340.23
SPO: 24-For children5.30 ± 1.495.27 ± 1.390.67
SPO: 25-Various4.87 ± 1.594.87 ± 1.430.72
SPO: 26-Soft (paddle, etc.)3.93 ± 1.824.23 ± 1.700.74
SPO: 27-Tennis3.20 ± 2.113.14 ± 2.030.82
SPO: 28-Pool5.30 ± 1.955.10 ± 1.770.81
SPO: 29-Others4.13 ± 2.014.03 ± 1.990.84
LEI: 30-Serv. complem. (parking)5.03 ± 1.835.17 ± 1.460.82
LEI: 31-Pub/discos/cafes3.97 ± 1.474.10 ± 1.540.76
LEI: 32-Excursions4.27 ± 1.724.37 ± 1.520.61
LEI: 33-Culture/shows4.00 ± 1.724.20 ± 1.830.71
LEI: 34-Shops4.43 ± 1.704.27 ± 1.760.65
LEI: 35-Street animation4.43 ± 1.774.77 ± 1.520.61
LEI: 36-Tourist office4.57 ± 1.774.90 ± 1.630.78
PUB: 37-Cleaning5.50 ± 1.415.33 ± 1.370.73
PUB: 38-Absence of beggars5.87 ± 1.315.37 ± 1.450.63
PUB: 39-Feeling of security5.30 ± 1.375.10 ± 1.400.58
PUB: 40-Little noise4.87 ± 1.414.57 ± 1.430.36
PUB: 41-Traffic4.23 ± 1.704.60 ± 1.330.38
PUB: 42-Urban furniture4.77 ± 1.504.93 ± 1.170.55
PUB: 43-Attitude of the people5.80 ± 1.305.47 ± 1.430.67
PUB: 44-Urban landscape5.10 ± 1.815.20 ± 1.450.64
PUB: 45-Natural environment5.67 ± 1.305.47 ± 1.380.67
INF: 46-Taxis3.93 ± 1.684.03 ± 1.500.73
INF: 47-Buses3.93 ± 1.704.07 ± 1.460.73
INF: 48-Roads5.07 ± 1.015.20 ± 1.100.64
INF: 49-Line buses3.93 ± 1.784.63 ± 1.710.69
INF: 50-Post office4.43 ± 1.794.70 ± 1.580.69
INF: 51-Public telephones4.53 ± 1.504.63 ± 1.400.61
INF: 52-Banking services4.60 ± 1.754.87 ± 1.430.71
a: Corresponds to 30 users, 12 men and 18 women, aged 15–88 years, not belonging to the main sample, and from three spas in the study. The concordance of the variables age, sex, occupation, income, transport used, and type of accommodation was virtually 100% (results not shown in the table). b: ACC: accommodation; CAT: catering; SPA: spa; SPO: sports facilities; LEI: leisure/culture/shopping; PUB: public roads/urban and natural environment; INF: transport infrastructures and other services. c: Arithmetic mean ± SD unweighted only with valid cases. d: ICC between cases with valid response. It can be assessed with the Landis and Koch scale [43]: <0.00: poor; 0.00–0.20: slight (slightly); 0.21–0.40: fair (fair); 0.41–0.60: moderate (moderate); 0.61–0.80: substantial (substantial); 0.81–1.00: almost perfect (almost perfect). Based on 29 individuals, since one did not respond in the repetition.
Table 6. Diagnostic criteria of the factor analysis (n = 725) a.
Table 6. Diagnostic criteria of the factor analysis (n = 725) a.
CriterionOutcome
Kaiser–Meyer–Olkin (adequacy measure)0.923
Bartlett’s sphericity testχ2 = 18082.2 (1326 gl.), p < 0.001
a Note that it is a weighted sample.
Table 7. Factor analysis (n = 725). Factors (8) after Varimax rotation, with all weights a.
Table 7. Factor analysis (n = 725). Factors (8) after Varimax rotation, with all weights a.
Factors (from Factor Analysis) and Approximate Interpretation
Questionnaire
Dimension b: Item-Item No.
1
≈SPO
2
≈ACC
3
≈LEI
4
≈SPA
5
≈PUB
6
≈CAT
7
≈PUB
8
≈INF
ACC: 01-Cleaning, hygiene0.080.770.040.180.160.190.060.10
ACC: 02-Security0.040.750.100.160.200.150.06−0.06
ACC: 03-Environmental quality−0.120.640.000.100.230.21−0.060.13
ACC: 04-Room0.190.680.060.170.100.150.030.20
ACC: 05-Aesthetics−0.050.750.170.260.160.13−0.010.09
ACC: 06-Facility for children and disabled.0.550.340.210.28−0.040.04−0.19−0.12
ACC: 07-Personal professionalism0.220.660.050.140.190.280.060.00
ACC: 08-Quality/price ratio0.240.580.070.310.020.280.110.04
CAT: 09-Variety0.230.230.190.280.150.720.030.08
CAT: 10-Cleaning0.080.290.130.280.200.700.010.12
CAT: 11-Environment−0.020.350.050.250.330.570.07−0.15
CAT: 12-Food quality0.230.240.080.270.160.740.050.08
CAT: 13-Personal professionalism0.150.340.070.090.170.690.12−0.07
CAT: 14-Ratio quality/price0.230.260.120.290.080.730.130.04
SPA: 15-Medicinal mineral waters−0.070.270.060.600.190.230.12−0.04
SPA: 16-Cleaning0.020.310.040.650.160.220.100.08
SPA: 17-Security0.020.280.070.430.310.240.15−0.22
SPA: 18-Equipment0.360.130.120.600.020.12−0.100.17
SPA: 19-Environment0.030.340.120.520.250.150.050.02
SPA: 20-Medical service0.040.360.140.340.180.250.220.00
SPA: 21-Handicapped access0.570.080.140.38−0.040.17−0.25−0.01
SPA: 22-Personal professionalism0.310.26−0.170.550.100.190.100.04
SPA: 23-Value for money0.310.19−0.030.610.010.210.110.09
SPO: 24-For children0.550.100.320.280.030.41−0.12−0.01
SPO: 25-Various0.760.05−0.100.330.010.02−0.060.07
SPO: 26-Soft (paddle, etc.)0.790.060.100.04−0.030.080.210.09
SPO: 27-Tennis0.80−0.030.08−0.100.060.150.260.06
SPO: 28-Pool0.760.07−0.17−0.060.210.11−0.100.13
SPO: 29-Others0.840.000.11−0.060.030.140.200.01
LEI: 30-Serv. complem. (parking)0.510.400.360.120.020.030.050.02
LEI: 31-Pub/discos/cafes0.010.120.870.050.050.090.100.03
LEI: 32-Excursions0.270.160.61−0.070.190.32−0.040.27
LEI: 33-Culture/show0.200.140.70−0.050.130.300.060.21
LEI: 34-Shops0.050.100.870.030.000.070.140.07
LEI: 35-Street animation0.000.080.550.320.360.15−0.030.17
LEI: 36-Tourist office−0.180.020.590.14−0.13−0.150.230.10
PUB: 37-Cleaning0.260.010.030.230.590.12−0.020.22
PUB: 38-Absence of beggars0.070.130.040.370.630.040.05−0.16
PUB: 39-Feeling of security0.030.220.140.190.740.110.08−0.17
PUB: 40-Little noise−0.140.29−0.06−0.050.750.12−0.050.19
PUB: 41-Traffic0.040.29−0.07−0.030.770.180.010.02
PUB: 42-Urban furniture0.340.120.010.020.580.130.280.29
PUB: 43-Attitude of the people0.030.160.240.460.430.16−0.01−0.08
PUB: 44-Urban landscape−0.020.090.270.430.570.140.100.12
PUB: 45-Natural environment−0.120.090.230.520.460.150.130.02
INF: 46-Taxis0.060.150.220.100.170.01−0.190.76
INF: 47-Buses0.140.110.160.07−0.07−0.030.150.81
INF: 48-Roads−0.090.110.200.18−0.020.170.700.05
INF: 49-Line buses0.060.100.23−0.040.120.160.490.57
INF: 50-Post office0.33−0.040.610.130.050.020.540.07
INF: 51-Public telephones0.300.120.270.100.27−0.010.65−0.01
INF: 52- Banking services0.38−0.060.550.130.070.130.47−0.09
a: The largest weights in each variable are in bold. Note that it is a weighted sample. b: ACC: accommodation; CAT: catering; SPA: spa; SPO: sports facilities; LEI: leisure/culture shopping; PUB: public roads/urban and natural environment; INF: transport infrastructures and other services.
Table 8. Variance explained with rotated factors (Varimax rotation) (n = 725) a.
Table 8. Variance explained with rotated factors (Varimax rotation) (n = 725) a.
FactorInterpretation ApproximateSum of Saturations Squared after Rotation
Total% Variance% Accumulated
1≈Sports facilities5.711.011.0
2≈Accommodation5.210.121.1
3≈Leisure, culture, and shopping4.79.130.2
4≈Spa4.68.939.1
5≈Public roads, urban and natural environment4.68.848.0
6≈Catering (meals)4.38.256.2
7≈Infraest. transp. and other serv. (roads, public phones) b2.44.660.8
8≈Infraest. transp. and other serv. (taxis, buses) b2.24.365.1
a: Note that it is a weighted sample. b: Note that regarding the theoretical model (questionnaire) in the factor analysis, the dimension “infraest. transp. and other serv.” is divided into two factors (7 and 8).
Table 9. Construct validity: association (Pearson correlation) a between global assessment by dimension, general satisfaction and score of the main questionnaire (n = 725).
Table 9. Construct validity: association (Pearson correlation) a between global assessment by dimension, general satisfaction and score of the main questionnaire (n = 725).
VariableScores by Dimension Calculated with All ItemsTotal Score
ACCCATSPASPOLEIPUBINF
Global valuation items b
ACC0.84
CAT 0.91
SPA 0.81
SPO 0.85
LEI 0.9
PUB 0.76
INF 0.69
Satisf. general0.570.620.670.290.40.530.380.7
Satisf. price quality0.450.480.620.310.260.350.330.54
a: In the table, Pearson’s linear correlation coefficient (r) is shown. Note that it is a weighted sample and that the cases with NS/NC values of the variables in the first column are excluded (see Table 10 **). In all cases, it was statistically significant (p < 0.001). b: ACC: accommodation; CAT: catering; SAP: spa; SPO: sports facilities; LEI: leisure/culture/shopping; PUB: public roads/urban and natural environment; INF: transport infrastructures and other services.
Table 10. Construct validity: association (comparison of means) between the variables of future intention and the scores estimated a.
Table 10. Construct validity: association (comparison of means) between the variables of future intention and the scores estimated a.
VariableScores by Dimension b Calculated with all ItemsPunctuation
Total
ACCCATSPASPOLEIPUBINF
Intent to recommend it
Yes5.95 ± 0.866.02 ± 1.056.01 ± 0.753.95 ± 1.714.52 ± 1.416.00 ± 0.944.31 ± 1.225.47 ± 0.75
No do not know5.15 ± 0.994.59 ± 1.394.86 ± 1.092.68 ± 1.253.16 ± 1.325.00 ± 1.203.11 ± 1.294.29 ± 0.90
p-value cp < 0.001p < 0.001p < 0.001p < 0.001p < 0.001p < 0.001p < 0.001p < 0.001
Intention to return
Next year/this year6.03 ± 0.856.13 ± 1.036.12 ± 0.723.93 ± 1.764.54 ± 1.456.12 ± 0.884.32 ± 1.245.56 ± 0.73
ever/do not know/never5.45 ± 0.875.25 ± 1.185.36 ± 0.923.61 ± 1.564.05 ± 1.365.35 ± 1.093.92 ± 1.304.85 ± 0.87
p-value c<0.001<0.001<0.0010.019<0.001<0.001<0.001<0.001
a: Based on the 725 users, since they all answered the questions in column 1. The sample sizes are in Table 7 **. The mean ± standard deviation is given in the table. Note that it is a weighted sample. b: ACC: accommodation; CAT: catering; SPA: spa; SPO: sports facilities; LEI: leisure/culture/shopping; PUB: public roads/urban and natural environment; INF: transport infrastructures and other services. c: T-student for independent samples, with weighted sample
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anaya-Aguilar, R.; Gemar, G.; Anaya-Aguilar, C. Validation of a Satisfaction Questionnaire on Spa Tourism. Int. J. Environ. Res. Public Health 2021, 18, 4507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094507

AMA Style

Anaya-Aguilar R, Gemar G, Anaya-Aguilar C. Validation of a Satisfaction Questionnaire on Spa Tourism. International Journal of Environmental Research and Public Health. 2021; 18(9):4507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094507

Chicago/Turabian Style

Anaya-Aguilar, Rosa, German Gemar, and Carmen Anaya-Aguilar. 2021. "Validation of a Satisfaction Questionnaire on Spa Tourism" International Journal of Environmental Research and Public Health 18, no. 9: 4507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094507

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop