cut variable within a loop

Thank you all very much. A combination of the solutions suggested solved my
problem!

2015-08-16 22:31 GMT+02:00 David Winsemius <dwinsemius at comcast.net>:
On Aug 16, 2015, at 8:57 AM, Janka VANSCHOENWINKEL wrote:

Hi David,

Thanks for your comment. I'll explain what I want to do. I explained it
already earlier but the explanation might have gone lost in some of the
emails.

I now see that you did explain that you wanted the positional matching in
cut2 as a "break". The code runs without error on my machine, but delivers
a lot of warnings about masking. You are repeatedly using attach on the
same named objects. Using `attach` in programming is generally not a good
idea. In interactive use it is safer to use `with`, although that is not
generally considered safe in programming, either.

You need to do a better job of nailing down the source of the difficulty
what ever it might be. While you say the cut2 function "doesn't work", you
don't actually give evidence of "failure".

It's fairly simple to show that your theory about why your code fails in
some way as being due to cut2 failing to accept an "i" value inside an
lapply call is just wrong:

o <- lapply(1:3, function(i) { cut2( 0:10, i) } )
o
[[1]]
 [1]  0      [ 1,10] [ 1,10] [ 1,10] [ 1,10] [ 1,10] [ 1,10]
 [8] [ 1,10] [ 1,10] [ 1,10] [ 1,10]
Levels:  0 [ 1,10]

[[2]]
 [1] [ 0, 2) [ 0, 2) [ 2,10] [ 2,10] [ 2,10] [ 2,10] [ 2,10]
 [8] [ 2,10] [ 2,10] [ 2,10] [ 2,10]
Levels: [ 0, 2) [ 2,10]

[[3]]
 [1] [ 0, 3) [ 0, 3) [ 0, 3) [ 3,10] [ 3,10] [ 3,10] [ 3,10]
 [8] [ 3,10] [ 3,10] [ 3,10] [ 3,10]
Levels: [ 0, 3) [ 3,10]

You also have two different definitions of weight2 for your irrigation
model:

Alldata_Irrigation$weight2<-Alldata_Irrigation$sys02*Alldata_Irrigation$se025
Alldata_Irrigation$weight2<-Alldata_Irrigation$b48+Alldata_Irrigation$b50

--
David
The variable irrigation ranges from 0 to 100. (maybe not in de small
sample I gave, but in reality I have over 60000 observations and there the
variable ranges from 0 to 100). I want to make (and use) 100 different
samples. The sample is based each time on the "i" that I put at the
beginning of the loop.
So:

i = 1: this means there are 2 subsets. One from 0-1, another from 1-100
i = 2: this means there are 2 subsets. One from 0-2, another from 2-100
i = 3: this means there are 2 subsets. One from 0-3, another from 3-100
i = 4: this means there are 2 subsets. One from 0-4, another from 4-100
...
i = 96: this means there are 2 subsets. One from 0-96, another from
96-100
i = 97: this means there are 2 subsets. One from 0-97, another from
97-100
i = 98: this means there are 2 subsets. One from 0-98, another from
98-100
i = 99: this means there are 2 subsets. One from 0-99, another from
99-100
It might be possible that i = 1 and i = 2 give the same results in the
small dataset. But in the full dataset all numbers are represented.
The cut2 function is capable of "cutting" a sample based on a number
supplied. Yet, when I tell him this number is "i", then it doesn't work. If
instead I write that the number is 10, then it does work and it gives me 2
subsets from 0-10 and from 10-100.
Hope this is more clear!

Janka

2015-08-14 20:10 GMT+02:00 David Winsemius <dwinsemius at comcast.net>:

When using a function in R you may need to supply an argument name. Are
you expecting this to be the number of groups. I cannot decipher the intent
here with such sparse commentary, but this call to `cut2` does not make
sense to me. Perhaps you meant the number of groups? .... in which case you
need  cut2( Alldata$irrigation, g=i ), since the arguments to cut2 are not
that same as the arguments to cut.
At the moment you are implicitly sending on the first pass a 1 and then
on the second pass a 2 to the second argument of cut2 which is the `breaks`
argument. So you wold be getting two different factors each with different
cut-point levels. I looked at your data and in point of fact there would be
no difference since you have 29 zero values and no values between 0 and 1.

table(cut2(dat$irrigation, 1))
        0 [  1,100]
       29        21
table(cut2(dat$irrigation, 2))
        0 [  2,100]
       29        21

 levels(Alldata$irri)<-c("0","1")

 Alldata_Rainfed<-subset(Alldata, irri == 0)
 Alldata_Irrigation<-subset(Alldata, irri == 1)

 Alldata_Rainfed$w<-Alldata_Rainfed$b48+Alldata_Rainfed$b50
 Alldata_Irrigation$w<-Alldata_Irrigation$b48+Alldata_Irrigation$b50

 OLS_Rainfed <- lm(LnALVperHA~ps1+ps2+ps3+ps4+ts1+ts2+ts3+ts4+
                     ps1sq+ps2sq+ps3sq+ps4sq+ts1sq+ts2sq+ts3sq+ts4sq+
                     pdnsty+portsML+cities500k+rentedland+subsidies1+
                     elevmean+elevrange+
                     t_gravel+t_gravel+t_ph_h2o+t_silt+t_sand+
                     AT+BE+DK+ES+FI+FR+GR+IE+IT+LU+NL+PT+SE+WDE+EDE+UK,
                   weights=w,Alldata_Rainfed)

 attach(Alldata_Rainfed)

 CoefRainfed_ps1 <- OLS_Rainfed$coeff[2]
 CoefRainfed_ps2 <- OLS_Rainfed$coeff[3]
 CoefRainfed_ps3 <- OLS_Rainfed$coeff[4]
 CoefRainfed_ps4 <- OLS_Rainfed$coeff[5]
 CoefRainfed_ts1 <- OLS_Rainfed$coeff[6]
 CoefRainfed_ts2 <- OLS_Rainfed$coeff[7]
 CoefRainfed_ts3 <- OLS_Rainfed$coeff[8]
 CoefRainfed_ts4 <- OLS_Rainfed$coeff[9]
 CoefRainfed_ps1sq <- OLS_Rainfed$coeff[10]
 CoefRainfed_ps2sq <- OLS_Rainfed$coeff[11]
 CoefRainfed_ps3sq <- OLS_Rainfed$coeff[12]
 CoefRainfed_ps4sq <- OLS_Rainfed$coeff[13]
 CoefRainfed_ts1sq <- OLS_Rainfed$coeff[14]
 CoefRainfed_ts2sq <- OLS_Rainfed$coeff[15]
 CoefRainfed_ts3sq <- OLS_Rainfed$coeff[16]
 CoefRainfed_ts4sq <- OLS_Rainfed$coeff[17]

 attach(Alldata_Rainfed)

 ###### MARGINAL EFFECTS SEASONAL and YEARLY and REGIONAL (EU or
COUNTRY level)
 # Maar dit is dus de marginale impact per LnALVperHA?

 Alldata_Rainfed$MEts1 =
CoefRainfed_ts1+2*CoefRainfed_ts1sq*Alldata_Rainfed$ts1
 Alldata_Rainfed$MEts2 =
CoefRainfed_ts2+2*CoefRainfed_ts2sq*Alldata_Rainfed$ts2
 Alldata_Rainfed$MEts3 =
CoefRainfed_ts3+2*CoefRainfed_ts3sq*Alldata_Rainfed$ts3
 Alldata_Rainfed$MEts4 =
CoefRainfed_ts4+2*CoefRainfed_ts4sq*Alldata_Rainfed$ts4
 Alldata_Rainfed$MEt   = Alldata_Rainfed$MEts1 +
Alldata_Rainfed$MEts2 + Alldata_Rainfed$MEts3 + Alldata_Rainfed$MEts4

 Alldata_Rainfed$MEps1 =
CoefRainfed_ps1+2*CoefRainfed_ps1sq*Alldata_Rainfed$ps1
 Alldata_Rainfed$MEps2 =
CoefRainfed_ps2+2*CoefRainfed_ps2sq*Alldata_Rainfed$ps2
 Alldata_Rainfed$MEps3 =
CoefRainfed_ps3+2*CoefRainfed_ps3sq*Alldata_Rainfed$ps3
 Alldata_Rainfed$MEps4 =
CoefRainfed_ps4+2*CoefRainfed_ps4sq*Alldata_Rainfed$ps4
 Alldata_Rainfed$MEp   = Alldata_Rainfed$MEps1 +
Alldata_Rainfed$MEps2 + Alldata_Rainfed$MEps3 + Alldata_Rainfed$MEps4

 Alldata_Rainfed$weight2<-Alldata_Rainfed$b48+Alldata_Rainfed$b50
 attach(Alldata_Rainfed)
 library(stats)
 MEt_Rainfed<-weighted.mean(MEt,weight2)
 MEp_Rainfed<-weighted.mean(MEp,weight2)

 attach(Alldata_Irrigation)

 OLS_Irrigation <- lm(LnALVperHA~ps1+ps2+ps3+ps4+ts1+ts2+ts3+ts4+

 ps1sq+ps2sq+ps3sq+ps4sq+ts1sq+ts2sq+ts3sq+ts4sq+

 pdnsty+portsML+cities500k+rentedland+subsidies1+
                        elevmean+elevrange+
                        t_gravel+t_gravel+t_ph_h2o+t_silt+t_sand+

 AT+BE+DK+ES+FI+FR+GR+IE+IT+LU+NL+PT+SE+WDE+EDE+UK,
                      weights=w,Alldata_Irrigation)

 CoefIrrigation_ps1 <- OLS_Irrigation$coeff[2]
 CoefIrrigation_ps2 <- OLS_Irrigation$coeff[3]
 CoefIrrigation_ps3 <- OLS_Irrigation$coeff[4]
 CoefIrrigation_ps4 <- OLS_Irrigation$coeff[5]
 CoefIrrigation_ts1 <- OLS_Irrigation$coeff[6]
 CoefIrrigation_ts2 <- OLS_Irrigation$coeff[7]
 CoefIrrigation_ts3 <- OLS_Irrigation$coeff[8]
 CoefIrrigation_ts4 <- OLS_Irrigation$coeff[9]
 CoefIrrigation_ps1sq <- OLS_Irrigation$coeff[10]
 CoefIrrigation_ps2sq <- OLS_Irrigation$coeff[11]
 CoefIrrigation_ps3sq <- OLS_Irrigation$coeff[12]
 CoefIrrigation_ps4sq <- OLS_Irrigation$coeff[13]
 CoefIrrigation_ts1sq <- OLS_Irrigation$coeff[14]
 CoefIrrigation_ts2sq <- OLS_Irrigation$coeff[15]
 CoefIrrigation_ts3sq <- OLS_Irrigation$coeff[16]
 CoefIrrigation_ts4sq <- OLS_Irrigation$coeff[17]

 attach(Alldata_Irrigation)
 # gives the residual errors in Y
 Alldata_Irrigation$residuals <-resid(OLS_Irrigation)

 # gives the predicted values for Ln_Y
 Alldata_Irrigation$Ln_y_hat <-fitted(OLS_Irrigation)

 # Zelf functie rmse maken
 rmse <- function(error)
 {
   sqrt(mean(error^2))
 }
 Alldata_Irrigation$y_hat <-

exp(Alldata_Irrigation$Ln_y_hat)*exp(0.5*(rmse(OLS_Irrigation$residuals))^2)
 # absolute impact (landwaarde current)

Alldata_Irrigation$absolute.current<-Alldata_Irrigation$y_hat*Alldata_Irrigation$se025*Alldata_Irrigation$sys02

 ###### MARGINAL EFFECTS SEASONAL and YEARLY and REGIONAL (EU or
COUNTRY level)
 # Maar dit is dus de marginale impact per LnALVperHA?

 Alldata_Irrigation$MEts1 =
CoefIrrigation_ts1+2*CoefIrrigation_ts1sq*Alldata_Irrigation$ts1
 Alldata_Irrigation$MEts2 =
CoefIrrigation_ts2+2*CoefIrrigation_ts2sq*Alldata_Irrigation$ts2
 Alldata_Irrigation$MEts3 =
CoefIrrigation_ts3+2*CoefIrrigation_ts3sq*Alldata_Irrigation$ts3
 Alldata_Irrigation$MEts4 =
CoefIrrigation_ts4+2*CoefIrrigation_ts4sq*Alldata_Irrigation$ts4
 Alldata_Irrigation$MEt   = Alldata_Irrigation$MEts1 +
Alldata_Irrigation$MEts2 + Alldata_Irrigation$MEts3 +
Alldata_Irrigation$MEts4

 Alldata_Irrigation$MEps1 =
CoefIrrigation_ps1+2*CoefIrrigation_ps1sq*Alldata_Irrigation$ps1
 Alldata_Irrigation$MEps2 =
CoefIrrigation_ps2+2*CoefIrrigation_ps2sq*Alldata_Irrigation$ps2
 Alldata_Irrigation$MEps3 =
CoefIrrigation_ps3+2*CoefIrrigation_ps3sq*Alldata_Irrigation$ps3
 Alldata_Irrigation$MEps4 =
CoefIrrigation_ps4+2*CoefIrrigation_ps4sq*Alldata_Irrigation$ps4
 Alldata_Irrigation$MEp   = Alldata_Irrigation$MEps1 +
Alldata_Irrigation$MEps2 + Alldata_Irrigation$MEps3 +
Alldata_Irrigation$MEps4

Alldata_Irrigation$weight2<-Alldata_Irrigation$sys02*Alldata_Irrigation$se025

Alldata_Irrigation$weight2<-Alldata_Irrigation$b48+Alldata_Irrigation$b50
 attach(Alldata_Irrigation)
 library(stats)
 MEt_Irrigation<-weighted.mean(MEt,weight2)
 MEp_Irrigation<-weighted.mean(MEp,weight2)

 c(MEt_Rainfed,MEp_Rainfed,MEt_Irrigation,MEp_Irrigation)

 attach(Alldata)

 # And in the loop (index i):

 d[i, ] = c(MEt_Rainfed,MEp_Rainfed,MEt_Irrigation,MEp_Irrigation)

})
out<-as.data.frame(do.call(rbind, o))

And the data are:

structure(list(LnALVperHA = c(8.09964942932129, 9.53274631500244,
7.42697763442993, 8.25370121002197, 8.42619132995605, 8.0093936920166,
8.09785747528076, 8.49044704437256, 9.08215141296387, 8.38935947418213,
8.67814350128174, 8.38935947418213, 10.4056901931763, 8.48210144042969,
8.30281829833984, 8.92265796661377, 8.33178997039795, 4.54404163360596,
10.662184715271, 9.62167072296143, 7.98790407180786, 7.58244323730469,
7.23262739181519, 9.47037124633789, 8.93403625488281, 7.54256629943848,
9.40302467346191, 10.6290521621704, 8.59830188751221, 8.59585666656494,
9.10000514984131, 9.99381542205811, 9.54681301116943, 9.53055191040039,
8.67971229553223, 7.19780731201172, 8.90067958831787, 6.0509786605835,
6.55788946151733, 8.22567272186279, 9.05618953704834, 6.81858921051025,
8.46410751342773, 7.81292057037354, 8.38989448547363, 10.4709157943726,
8.06132888793945, 8.43629264831543, 10.3087100982666, 10.3218297958374
), ps1 = c(5.14855766296387, 4.71904611587524, 7.9462103843689,
10.6017990112305, 11.233078956604, 9.12952136993408, 12.6536712646484,
11.233078956604, 11.233078956604, 11.233078956604, 11.233078956604,
11.233078956604, 5.93759632110596, 10.6017990112305, 11.233078956604,
10.6017990112305, 7.95780467987061, 9.07744884490967, 4.29865598678589,
8.27481746673584, 3.25137901306152, 4.51061344146729, 6.34518480300903,
6.66202449798584, 6.66202449798584, 4.75249433517456, 6.28858852386475,
6.33270215988159, 10.3600759506226, 10.3600759506226, 18.7164611816406,
5.73318386077881, 7.92949104309082, 9.09823608398438, 11.233078956604,
10.4455404281616, 11.233078956604, 10.4455404281616, 10.4455404281616,
10.6017990112305, 9.19112777709961, 10.4455404281616, 11.233078956604,
11.064302444458, 11.233078956604, 5.93759632110596, 11.233078956604,
10.6017990112305, 6.05948448181152, 9.5645227432251), ps2 =
c(5.23111915588379,
4.86784505844116, 7.7175760269165, 4.34898376464844, 4.48626232147217,
9.57159423828125, 8.38174915313721, 4.48626232147217, 4.48626232147217,
4.48626232147217, 4.48626232147217, 4.48626232147217, 6.87198734283447,
4.34898376464844, 4.48626232147217, 4.34898376464844, 6.2098217010498,
7.5497522354126, 5.62545442581177, 5.57168531417847, 3.08954334259033,
6.6683931350708, 4.41767883300781, 6.11901044845581, 6.11901044845581,
4.06884765625, 6.35917854309082, 5.7121729850769, 8.55229663848877,
8.55229663848877, 11.8981914520264, 5.49351119995117, 5.34777498245239,
6.12420177459717, 4.48626232147217, 5.2967677116394, 4.48626232147217,
5.2967677116394, 5.2967677116394, 4.34898376464844, 4.51386308670044,
5.2967677116394, 4.48626232147217, 5.98725175857544, 4.48626232147217,
6.87198734283447, 4.48626232147217, 4.34898376464844, 5.58411026000977,
4.42436075210571), ps3 = c(4.95634937286377, 3.50353670120239,
6.01129817962646, 0.851324141025543, 0.816295921802521,
8.03804397583008,
5.56230783462524, 0.816295921802521, 0.816295921802521,
0.816295921802521,
0.816295921802521, 0.816295921802521, 6.01666784286499,
0.851324141025543,
0.816295921802521, 0.851324141025543, 3.45424580574036,
5.31899690628052,
7.45753812789917, 3.34133338928223, 6.61472988128662, 11.244439125061,
2.19617891311646, 5.29748106002808, 5.29748106002808, 1.63307499885559,
5.51272773742676, 6.78562116622925, 4.5334997177124, 4.5334997177124,
4.31791353225708, 7.10963106155396, 2.32198905944824, 2.74845194816589,
0.816295921802521, 1.47570741176605, 0.816295921802521,
1.47570741176605,
1.47570741176605, 0.851324141025543, 1.39068424701691,
1.47570741176605,
0.816295921802521, 1.85064959526062, 0.816295921802521,
6.01666784286499,
0.816295921802521, 0.851324141025543, 6.78009986877441,
1.21070051193237
), ps4 = c(5.66667366027832, 4.82342433929443, 7.40090322494507,
6.59299898147583, 7.33758926391602, 9.98004341125488, 10.3958940505981,
7.33758926391602, 7.33758926391602, 7.33758926391602, 7.33758926391602,
7.33758926391602, 8.31999015808105, 6.59299898147583, 7.33758926391602,
6.59299898147583, 7.05771064758301, 8.38344383239746, 4.75349426269531,
9.00399303436279, 5.48189449310303, 5.9071044921875, 5.30881881713867,
8.68398857116699, 8.68398857116699, 4.32339859008789, 8.57950687408447,
6.78787326812744, 8.68624305725098, 8.68624305725098, 12.9021902084351,
6.14854049682617, 6.71301507949829, 7.50605535507202, 7.33758926391602,
8.11069011688232, 7.33758926391602, 8.11069011688232, 8.11069011688232,
6.59299898147583, 5.92181205749512, 8.11069011688232, 7.33758926391602,
9.29954528808594, 7.33758926391602, 8.31999015808105, 7.33758926391602,
6.59299898147583, 6.16447877883911, 5.83903217315674), ts1 =
c(4.19949150085449,
2.46556353569031, 3.96805644035339, 9.05560302734375, 9.5199556350708,
1.18671488761902, 6.60286664962769, 9.5199556350708, 9.5199556350708,
9.5199556350708, 9.5199556350708, 9.5199556350708, 2.12847352027893,
9.05560302734375, 9.5199556350708, 9.05560302734375, 2.11432313919067,
6.49393510818481, -0.165110915899277, 7.78503036499023,
-7.71160411834717,
-0.979450941085815, 4.96369075775146, 4.28496122360229,
4.28496122360229,
6.35976600646973, 3.02656149864197, 2.80754446983337, 5.94739389419556,
5.94739389419556, 8.70161914825439, 1.57025468349457, 5.08782005310059,
4.27688789367676, 9.5199556350708, 8.49832916259766, 9.5199556350708,
8.49832916259766, 8.49832916259766, 9.05560302734375, 6.33359289169312,
8.49832916259766, 9.5199556350708, 7.99740839004517, 9.5199556350708,
2.12847352027893, 9.5199556350708, 9.05560302734375, 2.67069268226624,
7.33829879760742), ts2 = c(9.89923763275146, 10.9084701538086,
9.61682415008545, 13.6253662109375, 13.8121919631958, 6.19518041610718,
9.40560817718506, 13.8121919631958, 13.8121919631958, 13.8121919631958,
13.8121919631958, 13.8121919631958, 10.3912172317505, 13.6253662109375,
13.8121919631958, 13.6253662109375, 9.77112770080566, 11.5460777282715,
8.18180465698242, 12.9412984848022, 2.54625177383423, 8.29829216003418,
10.6650953292847, 10.1770324707031, 10.1770324707031, 12.4333782196045,
8.98324680328369, 8.45312309265137, 9.23384857177734, 9.23384857177734,
11.371600151062, 8.09108352661133, 12.0714511871338, 11.385799407959,
13.8121919631958, 13.912787437439, 13.8121919631958, 13.912787437439,
13.912787437439, 13.6253662109375, 12.0018119812012, 13.912787437439,
13.8121919631958, 14.0190010070801, 13.8121919631958, 10.3912172317505,
13.8121919631958, 13.6253662109375, 8.53981018066406, 12.7294788360596
), ts3 = c(17.718994140625, 21.1172523498535, 17.8669090270996,
23.1215572357178, 22.9536685943604, 15.3891229629517, 15.7000684738159,
22.9536685943604, 22.9536685943604, 22.9536685943604, 22.9536685943604,
22.9536685943604, 20.1229286193848, 23.1215572357178, 22.9536685943604,
23.1215572357178, 19.8251171112061, 19.3250198364258, 16.8351039886475,
22.2966594696045, 14.6743259429932, 17.1554985046387, 20.1656894683838,
20.0012702941895, 20.0012702941895, 23.2738876342773, 18.6255321502686,
16.2553405761719, 16.551155090332, 16.551155090332, 17.6266174316406,
16.1711521148682, 22.280725479126, 21.450382232666, 22.9536685943604,
23.5616970062256, 22.9536685943604, 23.5616970062256, 23.5616970062256,
23.1215572357178, 22.1113948822021, 23.5616970062256, 22.9536685943604,
23.5085678100586, 22.9536685943604, 20.1229286193848, 22.9536685943604,
23.1215572357178, 16.3595314025879, 22.7737102508545), ts4 =
c(11.661883354187,
12.7669324874878, 11.6320190429688, 17.2357921600342, 17.4911460876465,
9.09537506103516, 12.179615020752, 17.4911460876465, 17.4911460876465,
17.4911460876465, 17.4911460876465, 17.4911460876465, 12.0781927108765,
17.2357921600342, 17.4911460876465, 17.2357921600342, 11.9486837387085,
13.7441387176514, 8.9575023651123, 15.9984045028687, 4.02816677093506,
9.12790489196777, 13.0505475997925, 12.842321395874, 12.842321395874,
14.8937959671021, 11.5566177368164, 10.0515727996826, 12.2921047210693,
12.2921047210693, 14.2251281738281, 9.64802074432373, 14.6072359085083,
13.7993869781494, 17.4911460876465, 17.0232067108154, 17.4911460876465,
17.0232067108154, 17.0232067108154, 17.2357921600342, 15.045259475708,
17.0232067108154, 17.4911460876465, 16.7633666992188, 17.4911460876465,
12.0781927108765, 17.4911460876465, 17.2357921600342, 10.0954942703247,
15.9187803268433), ps1sq = c(26.5076465606689, 22.2693958282471,
63.1422576904297, 112.398139953613, 126.182060241699, 83.3481597900391,
160.11540222168, 126.182060241699, 126.182060241699, 126.182060241699,
126.182060241699, 126.182060241699, 35.2550506591797, 112.398139953613,
126.182060241699, 112.398139953613, 63.3266563415527, 82.4000778198242,
18.478443145752, 68.4726028442383, 10.5714654922485, 20.3456344604492,
40.2613716125488, 44.3825721740723, 44.3825721740723, 22.58620262146,
39.5463447570801, 40.1031150817871, 107.331176757812, 107.331176757812,
350.305908203125, 32.8693962097168, 62.8768272399902, 82.7779006958008,
126.182060241699, 109.109313964844, 126.182060241699, 109.109313964844,
109.109313964844, 112.398139953613, 84.4768295288086, 109.109313964844,
126.182060241699, 122.418785095215, 126.182060241699, 35.2550506591797,
126.182060241699, 112.398139953613, 36.7173538208008, 91.480094909668
), ps2sq = c(27.3646068572998, 23.695915222168, 59.560977935791,
18.9136600494385, 20.1265487670898, 91.6154174804688, 70.2537155151367,
20.1265487670898, 20.1265487670898, 20.1265487670898, 20.1265487670898,
20.1265487670898, 47.2242088317871, 18.9136600494385, 20.1265487670898,
18.9136600494385, 38.5618858337402, 56.9987602233887, 31.6457366943359,
31.0436763763428, 9.54527759552002, 44.4674682617188, 19.5158863067627,
37.4422874450684, 37.4422874450684, 16.5555210113525, 40.439151763916,
32.6289215087891, 73.1417770385742, 73.1417770385742, 141.566955566406,
30.1786651611328, 28.5986976623535, 37.5058479309082, 20.1265487670898,
28.0557479858398, 20.1265487670898, 28.0557479858398, 28.0557479858398,
18.9136600494385, 20.3749599456787, 28.0557479858398, 20.1265487670898,
35.8471832275391, 20.1265487670898, 47.2242088317871, 20.1265487670898,
18.9136600494385, 31.1822872161865, 19.5749683380127), ps3sq =
c(24.5653991699219,
12.27476978302, 36.1357040405273, 0.72475278377533, 0.666339039802551,
64.6101531982422, 30.9392681121826, 0.666339039802551,
0.666339039802551,
0.666339039802551, 0.666339039802551, 0.666339039802551,
36.2002906799316,
0.72475278377533, 0.666339039802551, 0.72475278377533,
11.9318141937256,
28.2917289733887, 55.614875793457, 11.1645088195801, 43.7546501159668,
126.437408447266, 4.82320165634155, 28.063304901123, 28.063304901123,
2.6669340133667, 30.3901672363281, 46.0446548461914, 20.552619934082,
20.552619934082, 18.6443767547607, 50.5468521118164, 5.39163303375244,
7.55398797988892, 0.666339039802551, 2.17771244049072,
0.666339039802551,
2.17771244049072, 2.17771244049072, 0.72475278377533, 1.93400263786316,
2.17771244049072, 0.666339039802551, 3.42490386962891,
0.666339039802551,
36.2002906799316, 0.666339039802551, 0.72475278377533,
45.9697532653809,
1.46579575538635), ps4sq = c(32.1111907958984, 23.2654228210449,
54.7733688354492, 43.4676361083984, 53.840217590332, 99.6012649536133,
108.074615478516, 53.840217590332, 53.840217590332, 53.840217590332,
53.840217590332, 53.840217590332, 69.2222366333008, 43.4676361083984,
53.840217590332, 43.4676361083984, 49.811279296875, 70.2821273803711,
22.5957069396973, 81.071891784668, 30.0511665344238, 34.8938827514648,
28.183557510376, 75.4116592407227, 75.4116592407227, 18.6917762756348,
73.6079406738281, 46.0752220153809, 75.4508209228516, 75.4508209228516,
166.466506958008, 37.8045501708984, 45.0645713806152, 56.3408660888672,
53.840217590332, 65.7832946777344, 53.840217590332, 65.7832946777344,
65.7832946777344, 43.4676361083984, 35.0678596496582, 65.7832946777344,
53.840217590332, 86.4815444946289, 53.840217590332, 69.2222366333008,
53.840217590332, 43.4676361083984, 38.0007972717285, 34.094295501709
), ts1sq = c(17.6357288360596, 6.07900333404541, 15.7454719543457,
82.0039443969727, 90.6295547485352, 1.40829217433929, 43.5978469848633,
90.6295547485352, 90.6295547485352, 90.6295547485352, 90.6295547485352,
90.6295547485352, 4.53039932250977, 82.0039443969727, 90.6295547485352,
82.0039443969727, 4.47036218643188, 42.1711921691895,
0.0272616147994995,
60.6066970825195, 59.4688377380371, 0.95932412147522, 24.6382255554199,
18.3608932495117, 18.3608932495117, 40.4466247558594, 9.16007423400879,
7.88230609893799, 35.3714942932129, 35.3714942932129, 75.7181777954102,
2.46569967269897, 25.8859119415283, 18.2917709350586, 90.6295547485352,
72.2215957641602, 90.6295547485352, 72.2215957641602, 72.2215957641602,
82.0039443969727, 40.1143989562988, 72.2215957641602, 90.6295547485352,
63.9585418701172, 90.6295547485352, 4.53039932250977, 90.6295547485352,
82.0039443969727, 7.13259935379028, 53.8506278991699), ts2sq =
c(97.9949035644531,
118.994720458984, 92.4833068847656, 185.650604248047, 190.776641845703,
38.3802604675293, 88.465461730957, 190.776641845703, 190.776641845703,
190.776641845703, 190.776641845703, 190.776641845703, 107.977394104004,
185.650604248047, 190.776641845703, 185.650604248047, 95.4749374389648,
133.311904907227, 66.9419250488281, 167.477203369141, 6.48339796066284,
68.8616561889648, 113.744255065918, 103.571990966797, 103.571990966797,
154.588897705078, 80.6987228393555, 71.4552917480469, 85.2639617919922,
85.2639617919922, 129.313293457031, 65.4656295776367, 145.719940185547,
129.636428833008, 190.776641845703, 193.565658569336, 190.776641845703,
193.565658569336, 193.565658569336, 185.650604248047, 144.043487548828,
193.565658569336, 190.776641845703, 196.53239440918, 190.776641845703,
107.977394104004, 190.776641845703, 185.650604248047, 72.9283599853516,
162.039627075195), ts3sq = c(313.962768554688, 445.938354492188,
319.226440429688, 534.606384277344, 526.870910644531, 236.825103759766,
246.492156982422, 526.870910644531, 526.870910644531, 526.870910644531,
526.870910644531, 526.870910644531, 404.932250976562, 534.606384277344,
526.870910644531, 534.606384277344, 393.035278320312, 373.456390380859,
283.420715332031, 497.141021728516, 215.335845947266, 294.311126708984,
406.655029296875, 400.050811767578, 400.050811767578, 541.673828125,
346.910461425781, 264.236083984375, 273.940734863281, 273.940734863281,
310.697631835938, 261.506164550781, 496.430725097656, 460.118896484375,
526.870910644531, 555.153564453125, 526.870910644531, 555.153564453125,
555.153564453125, 534.606384277344, 488.913787841797, 555.153564453125,
526.870910644531, 552.652770996094, 526.870910644531, 404.932250976562,
526.870910644531, 534.606384277344, 267.63427734375, 518.641906738281
), ts4sq = c(135.999526977539, 162.994567871094, 135.303863525391,
297.072540283203, 305.940185546875, 82.7258453369141, 148.343017578125,
305.940185546875, 305.940185546875, 305.940185546875, 305.940185546875,
305.940185546875, 145.882736206055, 297.072540283203, 305.940185546875,
297.072540283203, 142.771041870117, 188.901351928711, 80.2368469238281,
255.948944091797, 16.2261276245117, 83.3186492919922, 170.316787719727,
164.925216674805, 164.925216674805, 221.825164794922, 133.555419921875,
101.034118652344, 151.095840454102, 151.095840454102, 202.354278564453,
93.0843048095703, 213.371337890625, 190.423080444336, 305.940185546875,
289.789581298828, 305.940185546875, 289.789581298828, 289.789581298828,
297.072540283203, 226.359832763672, 289.789581298828, 305.940185546875,
281.010467529297, 305.940185546875, 145.882736206055, 305.940185546875,
297.072540283203, 101.919006347656, 253.407562255859), pdnsty =
c(0.616999983787537,
0.0850000008940697, 0.068000003695488, 0.025000000372529,
0.0549999997019768,
0.0230000000447035, 0.133000001311302, 0.0549999997019768,
0.0549999997019768,
0.0549999997019768, 0.0549999997019768, 0.0549999997019768,
0.25900000333786,
0.025000000372529, 0.0549999997019768, 0.025000000372529,
0.0140000004321337,
0.14300000667572, 0.140000000596046, 0.777999997138977,
0.0329999998211861,
0.316000014543533, 0.0179999992251396, 0.105999998748302,
0.105999998748302,
0.046000000089407, 0.108000002801418, 0.310999989509583,
0.101000003516674,
0.101000003516674, 0.14300000667572, 0.168999999761581,
0.0439999997615814,
0.0379999987781048, 0.0549999997019768, 0.063000001013279,
0.0549999997019768,
0.063000001013279, 0.063000001013279, 0.025000000372529,
0.0640000030398369,
0.063000001013279, 0.0549999997019768, 0.209000006318092,
0.0549999997019768,
0.25900000333786, 0.0549999997019768, 0.025000000372529,
0.257999986410141,
0.0469999983906746), portsML = c(0.0900330692529678,
0.0604440234601498,
0.168490216135979, 0.275995850563049, 0.269018620252609,
0.175392478704453,
0.0350189469754696, 0.269018620252609, 0.269018620252609,
0.269018620252609,
0.269018620252609, 0.269018620252609, 0.11026918143034,
0.275995850563049,
0.269018620252609, 0.275995850563049, 0.145082741975784,
0.00440915673971176,
0.426146239042282, 0.0686663240194321, 0.103511147201061,
0.289726078510284,
0.234196603298187, 0.123688526451588, 0.123688526451588,
0.315173029899597,
0.112561739981174, 0.0461684986948967, 0.179993003606796,
0.179993003606796,
0.0438785217702389, 0.096462681889534, 0.0934395045042038,
0.121217466890812,
0.269018620252609, 0.212490051984787, 0.269018620252609,
0.212490051984787,
0.212490051984787, 0.275995850563049, 0.162760972976685,
0.212490051984787,
0.269018620252609, 0.270619571208954, 0.269018620252609,
0.11026918143034,
0.269018620252609, 0.275995850563049, 0.108705826103687,
0.196496397256851
), cities500k = c(0.0360943526029587, 0.0577861145138741,
0.183606043457985,
0.150749072432518, 0.185974538326263, 0.0923599153757095,
0.353672504425049,
0.185974538326263, 0.185974538326263, 0.185974538326263,
0.185974538326263,
0.185974538326263, 0.0887016654014587, 0.150749072432518,
0.185974538326263,
0.150749072432518, 0.144800990819931, 0.00326321297325194,
0.0622526630759239,
0.00816718116402626, 0.181859150528908, 0.163181975483894,
0.204970955848694,
0.129742562770844, 0.129742562770844, 0.0783679932355881,
0.0559677332639694,
0.0293320622295141, 0.248573184013367, 0.248573184013367,
0.174525216221809,
0.092569001019001, 0.176346719264984, 0.16088992357254,
0.185974538326263,
0.280431807041168, 0.185974538326263, 0.280431807041168,
0.280431807041168,
0.150749072432518, 0.088722825050354, 0.280431807041168,
0.185974538326263,
0.189705356955528, 0.185974538326263, 0.0887016654014587,
0.185974538326263,
0.150749072432518, 0.0712414756417274, 0.0842432081699371), rentedland
= c(0.571943998336792,
0, 0.5929936170578, 0, 0, 0.755691230297089, 0.440930217504501,
0, 0, 0, 0.229885056614876, 0, 0, 0, 0, 0, 0.890581607818604,
0.212423488497734, 0.386227518320084, 0, 0.11130790412426,
0.483032256364822,
0.444395005702972, 0, 0, 0.253378361463547, 0, 0.10909091681242,
0.181818187236786, 0.666666686534882, 0, 0.94951194524765,
0.846153855323792,
0.403846144676208, 0, 0, 0.155963316559792, 0, 0, 0.408163279294968,
0.699570834636688, 0, 0, 0, 0, 0, 0.0476190522313118, 0, 0, 0
), subsidies1 = c(361.835754394531, 0, 368.242034912109,
345.636352539062,
701.746032714844, 488.922821044922, 344.918609619141, 790.392150878906,
795.3125, 631.666687011719, 193.563217163086, 565.75, 0,
577.586181640625,
395.681823730469, 192, 371.963653564453, 9.9977331161499,
310.838317871094,
905.764709472656, 1745.76293945312, 359.003814697266, 163.204330444336,
427.94970703125, 204.842727661133, 52.2592887878418, 0, 0,
3022.24243164062,
80.2666702270508, 445.366577148438, 925.681640625, 824.769226074219,
625.192321777344, 850.441162109375, 280.891723632812, 619.266052246094,
333.962249755859, 376.304351806641, 317.551025390625, 166.652359008789,
171.224487304688, 526.119445800781, 253.191497802734, 334.470581054688,
107.277839660645, 431.428588867188, 0, 107.245544433594,
339.701507568359
), elevmean = c(0.121736958622932, 0.46412268280983, 0.344255149364471,
0.466430068016052, 0.43000802397728, 1.15364873409271,
0.0955904126167297,
0.43000802397728, 0.43000802397728, 0.43000802397728, 0.43000802397728,
0.43000802397728, 0.370405077934265, 0.466430068016052,
0.43000802397728,
0.466430068016052, 0.849120080471039, 0.0433186627924442,
0.335433751344681,
0.271958351135254, 0.125564843416214, 0.376024007797241,
0.815701544284821,
0.525435268878937, 0.525435268878937, 0.62959760427475,
0.518330037593842,
0.00362438289448619, 0.628515422344208, 0.628515422344208,
0.274942100048065,
0.0728112533688545, 0.496583759784698, 0.739268243312836,
0.43000802397728,
0.321640431880951, 0.43000802397728, 0.321640431880951,
0.321640431880951,
0.466430068016052, 0.585907399654388, 0.321640431880951,
0.43000802397728,
0.147326037287712, 0.43000802397728, 0.370405077934265,
0.43000802397728,
0.466430068016052, 0.0183117985725403, 0.414920538663864), elevrange =
c(0.180000007152557,
1.99300003051758, 0.611000001430511, 2.35199999809265,
2.29999995231628,
2.94199991226196, 0.354999989271164, 2.29999995231628,
2.29999995231628,
2.29999995231628, 2.29999995231628, 2.29999995231628, 2.01799988746643,
2.35199999809265, 2.29999995231628, 2.35199999809265, 1.7389999628067,
0.160999998450279, 0.314000010490417, 1.76300001144409,
0.17399999499321,
0.653999984264374, 1.63399994373322, 2.19099998474121,
2.19099998474121,
1.14100003242493, 1.34800004959106, 0.00899999961256981,
2.41300010681152,
2.41300010681152, 0.787999987602234, 0.26800000667572,
1.92200005054474,
2.02600002288818, 2.29999995231628, 1.05099999904633, 2.29999995231628,
1.05099999904633, 1.05099999904633, 2.35199999809265, 2.35999989509583,
1.05099999904633, 2.29999995231628, 0.772000014781952,
2.29999995231628,
2.01799988746643, 2.29999995231628, 2.35199999809265,
0.0649999976158142,
1.75399994850159), t_gravel = c(4.58953237533569, 13.3146963119507,
10.0136280059814, 13.8894920349121, 13.9366893768311, 13.5653190612793,
7.71220588684082, 13.9366893768311, 13.9366893768311, 13.9366893768311,
13.9366893768311, 13.9366893768311, 11.4818019866943, 13.8894920349121,
13.9366893768311, 13.8894920349121, 13.4321727752686, 5.71388387680054,
8.03888702392578, 9.01077747344971, 4.58924961090088, 8.14134693145752,
11.8983144760132, 9.96716785430908, 9.96716785430908, 11.1739711761475,
10.4019403457642, 5.16821479797363, 10.7357034683228, 10.7357034683228,
9.23897457122803, 4.3336238861084, 10.9520101547241, 12.9722995758057,
13.9366893768311, 13.1780118942261, 13.9366893768311, 13.1780118942261,
13.1780118942261, 13.8894920349121, 12.7335777282715, 13.1780118942261,
13.9366893768311, 12.315260887146, 13.9366893768311, 11.4818019866943,
13.9366893768311, 13.8894920349121, 6.68424606323242, 14.101095199585
), t_ph_h2o = c(6.07352828979492, 6.72695684432983, 5.60523176193237,
6.13967752456665, 6.86059141159058, 7.40929126739502, 5.68151950836182,
6.86059141159058, 6.86059141159058, 6.86059141159058, 6.86059141159058,
6.86059141159058, 6.51894521713257, 6.13967752456665, 6.86059141159058,
6.13967752456665, 6.98909568786621, 5.5628228187561, 6.68793487548828,
6.57724285125732, 4.67033195495605, 6.32772016525269, 6.4612717628479,
6.73934555053711, 6.73934555053711, 6.80293703079224, 6.17414236068726,
7.03696584701538, 5.93052577972412, 5.93052577972412, 5.43228578567505,
5.5989408493042, 6.86088180541992, 6.68706750869751, 6.86059141159058,
6.00043678283691, 6.86059141159058, 6.00043678283691, 6.00043678283691,
6.13967752456665, 6.89467239379883, 6.00043678283691, 6.86059141159058,
6.81896543502808, 6.86059141159058, 6.51894521713257, 6.86059141159058,
6.13967752456665, 5.63159275054932, 6.13170003890991), t_silt =
c(34.2329025268555,
33.4969100952148, 34.4774589538574, 27.8914813995361, 31.9258117675781,
39.6254501342773, 34.7939414978027, 31.9258117675781, 31.9258117675781,
31.9258117675781, 31.9258117675781, 31.9258117675781, 26.6626663208008,
27.8914813995361, 31.9258117675781, 27.8914813995361, 29.7444763183594,
21.3432540893555, 37.4038734436035, 28.1513748168945, 19.4936828613281,
33.5968360900879, 32.8024406433105, 33.313850402832, 33.313850402832,
28.3197917938232, 33.3154563903809, 38.103458404541, 36.0389099121094,
36.0389099121094, 34.9229164123535, 26.5577545166016, 30.9245643615723,
31.1334323883057, 31.9258117675781, 27.1493148803711, 31.9258117675781,
27.1493148803711, 27.1493148803711, 27.8914813995361, 31.3038387298584,
27.1493148803711, 31.9258117675781, 31.6541061401367, 31.9258117675781,
26.6626663208008, 31.9258117675781, 27.8914813995361, 15.6523361206055,
27.803352355957), t_sand = c(47.0063323974609, 37.0355186462402,
45.8286781311035, 36.0810203552246, 39.9931793212891, 39.3664970397949,
46.2948226928711, 39.9931793212891, 39.9931793212891, 39.9931793212891,
39.9931793212891, 39.9931793212891, 49.3508529663086, 36.0810203552246,
39.9931793212891, 36.0810203552246, 39.2436943054199, 65.7813262939453,
35.8039131164551, 51.2884674072266, 66.2952728271484, 46.6789817810059,
41.4505424499512, 44.4590721130371, 44.4590721130371, 48.7276763916016,
43.3654098510742, 33.999683380127, 43.040699005127, 43.040699005127,
43.2519073486328, 59.4827156066895, 43.8675765991211, 41.7124671936035,
39.9931793212891, 34.94921875, 39.9931793212891, 34.94921875,
34.94921875, 36.0810203552246, 39.1853942871094, 34.94921875,
39.9931793212891, 39.8589019775391, 39.9931793212891, 49.3508529663086,
39.9931793212891, 36.0810203552246, 75.7048721313477, 33.5687866210938
), AT = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), BE = c(0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0), DE = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), DK = c(0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0), ES = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), FI = c(0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0), FR = c(1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), GR = c(0, 1, 0,
1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
0, 1, 1, 0, 1), IE = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), IT = c(0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0), LU = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), NL = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0), PT = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), SE = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), WDE = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), EDE = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), UK = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), CY = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), BG = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), CZ = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), EE = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), HU = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), LT = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), LV = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), PL = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), RO = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), SI = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), SK = c(0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), b48 = c(70, 2.70000004768372, 63.9000015258789,
5.5, 6.30000019073486, 8.80000019073486, 48.0800018310547,
5.09999990463257,
6.40000009536743, 6, 6.69999980926514, 4, 6.30000019073486,
5.80000019073486,
8.80000019073486, 2, 13, 0.5, 10.25, 34, 65.2300033569336,
37.7799987792969,
74.9400024414062, 31.0200004577637, 20.0300006866455, 70.7200012207031,
40, 4.90000009536743, 13.5, 5, 26.8700008392334, 3, 2,
3.09999990463257,
6.80000019073486, 15.6999998092651, 9.19999980926514, 5.30000019073486,
4.59999990463257, 17.3999996185303, 7, 4.90000009536743,
13.3999996185303,
2.34999990463257, 8.5, 24.8700008392334, 4, 1.39999997615814,
34.7799987792969, 6.69999980926514), b50 = c(0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 34.2400016784668, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0), irrigation = c(0, 100, 0, 5.45454584062099,
7.9365074634552, 89.3392562866211, 0, 17.6470592617989, 0, 0,
65.5172407627106, 0, 61.904764175415, 34.4827562570572,
7.95454531908035,
75, 0, 0, 0, 0, 0, 0, 5.26393800973892, 0, 0, 0, 0, 0, 0, 0,
0, 0, 74.6153831481934, 84.6153914928436, 0, 5.09554147720337,
0, 0, 0, 21.0884347558022, 18.4549376368523, 6.1224490404129,
25.3731369972229, 2.12765969336033, 0, 84.3988716602325, 0, 0,
0, 100), awc_class = c(106.228088378906, 78.2306137084961,
80.9311141967773,
32.4921531677246, 54.8475151062012, 80.6665878295898, 116.331588745117,
54.8475151062012, 54.8475151062012, 54.8475151062012, 54.8475151062012,
54.8475151062012, 56.3101806640625, 32.4921531677246, 54.8475151062012,
32.4921531677246, 59.3034172058105, 101.193893432617, 96.5840377807617,
54.2786560058594, 87.1388244628906, 66.1907730102539, 57.205738067627,
55.4114303588867, 55.4114303588867, 80.9288787841797, 63.6008758544922,
150, 30.3404140472412, 30.3404140472412, 19.8318557739258,
104.236854553223,
79.2445755004883, 57.0045547485352, 54.8475151062012, 34.320426940918,
54.8475151062012, 34.320426940918, 34.320426940918, 32.4921531677246,
65.1337509155273, 34.320426940918, 54.8475151062012, 73.6748657226562,
54.8475151062012, 56.3101806640625, 54.8475151062012, 32.4921531677246,
127.726959228516, 27.9528160095215), sys02 = c(18.8571434020996,
303.529418945312, 30.2469139099121, 104.305557250977, 86.4935073852539,
51.25, 83.0927810668945, 453.118286132812, 42.5, 104.305557250977,
48.461540222168, 86.4935073852539, 55.1851844787598, 104.305557250977,
104.305557250977, 185.277770996094, 17.9775276184082, 25.2777786254883,
64, 21.6666660308838, 30, 24.2372875213623, 47.0285720825195,
16.1904754638672, 33.75, 22.5423736572266, 10.2857141494751,
39.230770111084, 6.06741571426392, 1, 28.3255805969238,
21.6000003814697,
69.2592620849609, 86.6666641235352, 48.5185203552246, 44.4186058044434,
48.6538467407227, 437.105255126953, 437.105255126953, 19.1666660308838,
48.461540222168, 437.105255126953, 48.6538467407227, 453.118286132812,
48.6538467407227, 14.2857141494751, 453.118286132812, 453.118286132812,
95.2380981445312, 63), se025 = c(163.529998779297, 2.70000004768372,
157, 5.5, 6.30000019073486, 36.0200004577637, 86, 5.09999990463257,
6.40000009536743, 6, 8.69999980926514, 4, 6.30000019073486,
5.80000019073486,
8.80000019073486, 2, 118.809997558594, 44.1100006103516,
16.7000007629395,
34, 73.4000015258789, 73.0800018310547, 134.880004882812,
31.0200004577637,
20.0300006866455, 94.7200012207031, 40, 5.5, 16.5, 15,
26.8700008392334,
59.4199981689453, 13, 5.19999980926514, 6.80000019073486,
15.6999998092651,
10.8999996185303, 5.30000019073486, 4.59999990463257, 29.3999996185303,
23.2999992370605, 4.90000009536743, 13.3999996185303, 2.34999990463257,
8.5, 24.8700008392334, 4.19999980926514, 1.39999997615814,
34.7799987792969,
6.69999980926514)), .Names = c("LnALVperHA", "ps1", "ps2", "ps3",
"ps4", "ts1", "ts2", "ts3", "ts4", "ps1sq", "ps2sq", "ps3sq",
"ps4sq", "ts1sq", "ts2sq", "ts3sq", "ts4sq", "pdnsty", "portsML",
"cities500k", "rentedland", "subsidies1", "elevmean", "elevrange",
"t_gravel", "t_ph_h2o", "t_silt", "t_sand", "AT", "BE", "DE",
"DK", "ES", "FI", "FR", "GR", "IE", "IT", "LU", "NL", "PT", "SE",
"WDE", "EDE", "UK", "CY", "BG", "CZ", "EE", "HU", "LT", "LV",
"PL", "RO", "SI", "SK", "b48", "b50", "irrigation", "awc_class",
"sys02", "se025"), row.names = c("2", "3", "4", "5", "6", "7",
"8", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20",
"21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31",
"32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42",
"43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53"
), class = "data.frame")

2015-08-14 14:58 GMT+02:00 PIKAL Petr <petr.pikal at precheza.cz>:
Hi Janka

Sorry, but we are limited in connecting to web services so I am not
able to restore your data and see your code. Result of dput(somedata)
coppied to email is preferable for sharing data and code can be copied to
email too. But do not use HTML as it usually scrambles  text.

Answer in line

From: Janka Vanschoenwinkel [mailto:janka.vanschoenwinkel at uhasselt.be
]
Sent: Friday, August 14, 2015 2:17 PM
To: Thierry Onkelinx; PIKAL Petr
Cc: r-help at r-project.org
Subject: Re: [R] cut variable within a loop

Hi Thierry and Petr,

I really appreciate the comments you already gave. Thank you very
much for that.

Below you can find a link to the data and the code. Hopefully this
helps in spotting the error.

I still think the issue is that the cut2 function only accepts
numbers, and not an "i" that refers to the number at the start of the loop.
To answer Petr his question, yes, column 3 and 4 are NA (these are the
columns of the second interval). But I don't really understand your point
so could you clarify this please?

If you use NA as a number of intervals you will get such errors

k<-c(2,4,NA,5)

ii<-vector(4, mode="list")

for (i in 1:4) {

ii[[i]] <- cut2(iris[,i], k[i])

}

Error in if (r[1] < cuts[1]) cuts <- c(r[1], cuts) :

 missing value where TRUE/FALSE needed

for (i in 1:4) {

ii[[i]] <- cut(iris[,i], k[i])

}

Error in cut.default(iris[, i], k[i]) : invalid number of intervals

If you remove NA from k definition error is gone.

k<-c(2,4,3,5)

ii<-vector(4, mode="list")

for (i in 1:4) {

ii[[i]] <- cut(iris[,i], k[i])

}

You can try it yourself. The error is not related to cycle; whenever
number of intervals in cut call is NA you always get an error.

Cheers

Petr

https://drive.google.com/folderview?id=0By9u5m3kxn9yfkxxeVNMdnRQQXhoT05CRlJlZVBCWWF2NURMMTNmVFVFeXJXXzhlMWE4SUk&usp=sharing

Thank you very much once again!

Janka

2015-08-11 15:10 GMT+02:00 Thierry Onkelinx <thierry.onkelinx at inbo.be
:
You'll need to send a reproducible example of the code. We can't run
the code that you send. Hence it is hard to help you. See e.g.
http://adv-r.had.co.nz/Reproducibility.html

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for
Nature and Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no
more than asking him to perform a post-mortem examination: he may be able
to say what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does
not ensure that a reasonable answer can be extracted from a given body of
data. ~ John Tukey

2015-08-11 14:57 GMT+02:00 Janka Vanschoenwinkel <
janka.vanschoenwinkel at uhasselt.be>:
Hi Thierry!

Thanks for your answer. I tried this, but I get this error:

"Error in cut.default(x, k2) : invalid number of intervals"

Which is strange because I am not specifying intervals, but the
number at where the sample has to be cut?

Greetings from Belgium! :-)

2015-08-11 14:52 GMT+02:00 Thierry Onkelinx <thierry.onkelinx at inbo.be
:
Dear Janka,

You loop goes for 0 to 100. It should probably go from 1:99

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for
Nature and Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no
more than asking him to perform a post-mortem examination: he may be able
to say what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does
not ensure that a reasonable answer can be extracted from a given body of
data. ~ John Tukey

2015-08-11 14:38 GMT+02:00 Janka Vanschoenwinkel <
janka.vanschoenwinkel at uhasselt.be>:
Dear list members,

I have a loop where I want to do several calculations for different
samples
and save the results for each sample. These samples are for each loop
different. I want to use the "i" in the loop to cut the samples.

So for instance:

  - In loop 1 (i=1), I have a sample from 0-1 and a sample from 1-100.
  - In loop 2 (i=2), I have a sample from 0-2 and a sample from 2-100.
  - In loop 99 (i=99), I have a sample from 0-99 and a sample from
99-100.
I built the following function, but there is *a problem with the cut2
function* since it doesn't recognize the "i". Outside the lapply loop
it
works, but not inside the loop.

Could somebody please help me with this problem? Thanks a lot!

d=data.frame(MEt_Rainfed=rep(0,100),MEp_Rainfed=rep(0,100),MEt_Irrigation=rep(0,100),MEp_Irrigation=rep(0,100))

   o<-lapply(0:100, function(i){

       Alldata$irri=cut2(Alldata$irrigation,i)

       levels(Alldata$irri)<-c("0","1")

      Alldata_Rainfed<-subset(Alldata, irri == 0)

      Alldata_Irrigation<-subset(Alldata, irri == 1)

   #calculations per sample, then store all the values per i and per
variable in a dataframe: (the calculations are not shown in this
example)

    d[i, ] = c(MEt_Rainfed,MEp_Rainfed,MEt_Irrigation,MEp_Irrigation)

  })

  out<-as.data.frame(do.call(rbind, o))

--
P Please consider the environment before printing this e-mail

       [[alternative HTML version deleted]]

______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

--

Mevrouw Janka Vanschoenwinkel
Doctoraatsbursaal - PhD
Milieueconomie - Environmental economics

T +32(0)11 26 87 42 | GSM +32(0)476 28 21 40

www.uhasselt.be/eec

Universiteit Hasselt | Campus Diepenbeek
Agoralaan Gebouw D | B-3590 Diepenbeek
Kantoor F11

Postadres: Universiteit Hasselt | Martelarenlaan 42 | B-3500 Hasselt

P Please consider the environment before printing this e-mail

--

Mevrouw Janka Vanschoenwinkel
Doctoraatsbursaal - PhD
Milieueconomie - Environmental economics

T +32(0)11 26 87 42 | GSM +32(0)476 28 21 40

www.uhasselt.be/eec

Universiteit Hasselt | Campus Diepenbeek
Agoralaan Gebouw D | B-3590 Diepenbeek
Kantoor F11

Postadres: Universiteit Hasselt | Martelarenlaan 42 | B-3500 Hasselt

P Please consider the environment before printing this e-mail

________________________________
Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a
jsou ur?eny pouze jeho adres?t?m.
Jestli?e jste obdr?el(a) tento e-mail omylem, informujte laskav?
neprodlen? jeho odes?latele. Obsah tohoto emailu i s p??lohami a jeho kopie
vyma?te ze sv?ho syst?mu.
Nejste-li zam??len?m adres?tem tohoto emailu, nejste opr?vn?ni tento
email jakkoliv u??vat, roz?i?ovat, kop?rovat ?i zve?ej?ovat.
Odes?latel e-mailu neodpov?d? za eventu?ln? ?kodu zp?sobenou
modifikacemi ?i zpo?d?n?m p?enosu e-mailu.

V p??pad?, ?e je tento e-mail sou??st? obchodn?ho jedn?n?:
- vyhrazuje si odes?latel pr?vo ukon?it kdykoliv jedn?n? o uzav?en?
smlouvy, a to z jak?hokoliv d?vodu i bez uveden? d?vodu.
- a obsahuje-li nab?dku, je adres?t opr?vn?n nab?dku bezodkladn?
p?ijmout; Odes?latel tohoto e-mailu (nab?dky) vylu?uje p?ijet? nab?dky ze
strany p??jemce s dodatkem ?i odchylkou.
- trv? odes?latel na tom, ?e p??slu?n? smlouva je uzav?ena teprve
v?slovn?m dosa?en?m shody na v?ech jej?ch n?le?itostech.
- odes?latel tohoto emailu informuje, ?e nen? opr?vn?n uzav?rat za
spole?nost ??dn? smlouvy s v?jimkou p??pad?, kdy k tomu byl p?semn? zmocn?n
nebo p?semn? pov??en a takov? pov??en? nebo pln? moc byly adres?tovi tohoto
emailu p??padn? osob?, kterou adres?t zastupuje, p?edlo?eny nebo jejich
existence je adres?tovi ?i osob? j?m zastoupen? zn?m?.

This e-mail and any documents attached to it may be confidential and
are intended only for its intended recipients.
If you received this e-mail by mistake, please immediately inform its
sender. Delete the contents of this e-mail with all attachments and its
copies from your system.
If you are not the intended recipient of this e-mail, you are not
authorized to use, disseminate, copy or disclose this e-mail in any manner.
The sender of this e-mail shall not be liable for any possible damage
caused by modifications of the e-mail or by delay with transfer of the
email.

In case that this e-mail forms part of business dealings:
- the sender reserves the right to end negotiations about entering
into a contract in any time, for any reason, and without stating any
reasoning.
- if the e-mail contains an offer, the recipient is entitled to
immediately accept such offer; The sender of this e-mail (offer) excludes
any acceptance of the offer on the part of the recipient containing any
amendment or variation.
- the sender insists on that the respective contract is concluded
only upon an express mutual agreement on all its aspects.
- the sender of this e-mail informs that he/she is not authorized to
enter into any contracts on behalf of the company except for cases in which
he/she is expressly authorized to do so in writing, and such authorization
or power of attorney is submitted to the recipient or the person
represented by the recipient, or the existence of such authorization is
known to the recipient of the person represented by the recipient.

--

Mevrouw Janka Vanschoenwinkel
Doctoraatsbursaal - PhD
Milieueconomie - Environmental economics

T +32(0)11 26 87 42 | GSM +32(0)476 28 21 40

www.uhasselt.be/eec

Universiteit Hasselt | Campus Diepenbeek
Agoralaan Gebouw D | B-3590 Diepenbeek
Kantoor F11

Postadres: Universiteit Hasselt | Martelarenlaan 42 | B-3500 Hasselt

P Please consider the environment before printing this e-mail

______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA

--

      Mevrouw Janka Vanschoenwinkel
Doctoraatsbursaal - PhD
Milieueconomie - Environmental economics

T +32(0)11 26 87 42 | GSM +32(0)476 28 21 40

www.uhasselt.be/eec

Universiteit Hasselt | Campus Diepenbeek
Agoralaan Gebouw D | B-3590 Diepenbeek
Kantoor F11

Postadres: Universiteit Hasselt | Martelarenlaan 42 | B-3500 Hasselt

P Please consider the environment before printing this e-mail

David Winsemius
Alameda, CA, USA

[image: Logo UHasselt]Mevrouw Janka Vanschoenwinkel
*Doctoraatsbursaal - PhD *
Milieueconomie - Environmental economics

T +32(0)11 26 87 42 | GSM +32(0)476 28 21 40

www.uhasselt.be/eec

Universiteit Hasselt | Campus Diepenbeek
Agoralaan Gebouw D | B-3590 Diepenbeek
Kantoor F11

Postadres: Universiteit Hasselt | Martelarenlaan 42 | B-3500 Hasselt

P Please consider the environment before printing this e-mail

	[[alternative HTML version deleted]]

cut variable within a loop

Thread (5 messages)