May I ask to solve this question in R?

An experiment to determine whether different colored candles (red, white, blue, yellow) burn at different speeds was conducted by Hsing-Chuan Tsai, Mei-Chiao Yang, Derek Wheeler, and Tom Schultz in 1989. Each experimenter collected four observations on each color in a random order, and "experimenter" was used as a blocking factor. Thus, the design was a general complete block design with v = 4, k = 16, b = 4, and s = 4. The resulting burning times (in seconds) are shown in Table 17.21. A pilot experiment indicated that treatments and blocks do interact. The candles used in the experiment were cake candles made by a single manufacturer. Analyze the experiment as though the experimenters represent a random sample from a large population of people who might use these candles in practice. Use a two-way mixed model with interaction.

Table 17.21 Data for the candle experiment (seconds) Person Color Red White Blue Yellow 1 989 1032 1044 979 1011 951 974 998 1077 1019 987 1031 928 1022 1033 1041 2 899 912 847 880 899 800 886 859 911 943 879 830 820 812 901 907 3 898 840 840 952 909 790 950 992 955 1005 961 915 871 905 920 890 4 993 957 987 960 864 925 949 973 1005 982 920 1001 824 790 978 938