[R] cut variable within a loop

[Janka Vanschoenwinkel]

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?

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.
>>>>
>>>
>>>
>>
>>
>> --
>>
>> [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
>>
>>
>


-- 

[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]]