[R] Error running caret's gbm train function with new version of caret

2013-05-04 Thread Katrina Bennett
I am running caret for model exploration. I developed my code a number of months ago and I've been running it with no issues. Recently, I updated my version of caret however, and now I am getting a new error. I'm wondering if this is due to the new release. The error I am getting is when I am

[R] Find position of asymptote

2012-04-19 Thread Katrina Bennett
Hi all, I would like to find the x position of an two asymptotes. Here is a sample of what I would like to do: x - seq(1, 153,, 153) a - 85 m - 65 s =-1.5 fn - function (x, a, m, s) { a * (exp((m - x)/s) * (1/s))/((1 + exp((m - x)/s)))^2 } plot.deriv1 - fn(1:153, a, m, s) I can find the

Re: [R] Find position of asymptote

2012-04-19 Thread Katrina Bennett
be appreciated. Sincerely, Katrina On Thu, Apr 19, 2012 at 4:34 AM, David Winsemius dwinsem...@comcast.net wrote: On Apr 19, 2012, at 4:41 AM, Katrina Bennett wrote: Hi all, I would like to find the x position of an two asymptotes. Here is a sample of what I would like to do: x - seq(1, 153,, 153

Re: [R] Find position of asymptote

2012-04-19 Thread Katrina Bennett
Dear R Help, Sorry I wasn't more clear before. Here is another crack at this. What I am still trying to do is estimate the point on a line when the slope changes or asymptotes. I have found some similar postings talking about this but no answers.

Re: [R] Sum matrix by rows, conditional on value

2011-11-29 Thread Katrina Bennett
Yes, that worked for me. Thank you. On Mon, Nov 28, 2011 at 10:16 PM, David Winsemius dwinsem...@comcast.netwrote: On Nov 29, 2011, at 1:08 AM, Katrina Bennett wrote: I'd like to sum a matrix only where the matrix meets a specific condition. The matrix has 365 rows and about 50,000

Re: [R] Retain parts of a matrix

2011-11-28 Thread Katrina Bennett
) get(t.mean.1.c) + (x * -0.0065) ) Finally, I want to output everything into one single matrix array. I hope that is more clear. Thanks for your assistance. Katrina On Sun, Nov 27, 2011 at 7:35 PM, David Winsemius dwinsem...@comcast.netwrote: On Nov 27, 2011, at 10:15 PM, Katrina

Re: [R] Retain parts of a matrix

2011-11-28 Thread Katrina Bennett
switched it in my example to my actual value (since I am going to use 300m for this initial attempt at this). I will need to make it a lot more complicated eventually! On Mon, Nov 28, 2011 at 6:57 PM, Katrina Bennett kebenn...@alaska.eduwrote: Hi David and Jim! Thanks very much for your help

[R] Sum matrix by rows, conditional on value

2011-11-28 Thread Katrina Bennett
I'd like to sum a matrix only where the matrix meets a specific condition. The matrix has 365 rows and about 50,000 columns. str(cdem.mat.yr) num [1:365, 1:41772] -43.5 -48.4 -45.9 -38.4 -32 ... I'm having trouble replicating this because it is not working out for me, so I'm unsure I can

[R] Retain parts of a matrix

2011-11-27 Thread Katrina Bennett
Hi all, I'm working to apply a function that will generate a matrix of results only when a specific criteria is met. I want my final results to be a matrix with both the values that meet the criteria (the results of the function), and those that to do in the same positions in the matrix (the

Re: [R] Obtaining a derivative of nls() SSlogis function

2011-11-19 Thread Katrina Bennett
, 2011 at 6:02 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Thu, Nov 17, 2011 at 4:40 PM, Katrina Bennett kebenn...@alaska.edu wrote: Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its

[R] Obtaining a derivative of nls() SSlogis function

2011-11-17 Thread Katrina Bennett
Hello, I am wondering if someone can help me. I have the following function that I derived using nls() SSlogis. I would like to find its derivative. I thought I had done this using deriv(), but for some reason this isn't working out for me. Here is the function: asym - 84.951 xmid - 66.90742 scal

[R] Seasonal and 11-day subset for zoo object

2011-09-07 Thread Katrina Bennett
I have a zooreg object and I want to be able to generate a value for seasons and 11-day composites paste it onto my zoo data frame, along with year, month and days. Right now I have the following to work from: eg. dat.zoo.mdy - with(month.day.year(time(dat.zoo)), cbind(dat.zoo, year, month, day,

Re: [R] sapply to bind columns, with repeat?

2011-08-13 Thread Katrina Bennett
On Fri, Aug 12, 2011 at 2:35 AM, Katrina Bennett kebenn...@alaska.edu wrote: Hi R-help, I am working with US COOP network station data and the files are concatenated in single rows for all years, but I need to pull these apart into rows for each day. To do this, I need to extract part of each

[R] sapply to bind columns, with repeat?

2011-08-12 Thread Katrina Bennett
Hi R-help, I am working with US COOP network station data and the files are concatenated in single rows for all years, but I need to pull these apart into rows for each day. To do this, I need to extract part of each row such as station id, year, mo, and repeat this against other variables in the

Re: [R] nls, how to determine function?

2011-08-09 Thread Katrina Bennett
Hi Petr, thanks for your help on this. I will most definitely get this book as it appears to be a good one. I am happy with how nls() and the self starting function appears to be fitting the data set. What I am wondering about is how to write out this function outside of R. For example, if I sat

[R] nls, how to determine function?

2011-08-08 Thread Katrina Bennett
Hi R help, I am trying to determine how nls() generates a function based on the self-starting SSlogis and what the formula for the function would be. I've scoured the help site, and other literature to try and figure this out but I still am unsure if I am correct in what I am coming up with.

Re: [R] Translate Sine Function in R?

2011-08-05 Thread Katrina Bennett
Hi Petr, Thank you for the response. I am not familiar with the nls package. I am interested in obtaining a function to represent the relationship observed in dat over x.seq (days) so that I can calculate a minimum and maximum of the function. In the example below, it seems like you are using a

Re: [R] Translate Sine Function in R?

2011-08-05 Thread Katrina Bennett
On Fri, Aug 5, 2011 at 4:41 PM, Katrina Bennett kebenn...@alaska.eduwrote: Hi Petr, Thank you for the response. I am not familiar with the nls package. I am interested in obtaining a function to represent the relationship observed in dat over x.seq (days) so that I can calculate a minimum

[R] Translate Sine Function in R?

2011-08-04 Thread Katrina Bennett
Hello, I'm trying to generate a sine wave in R to fit my observations using the general formula: y=a*sin(b[x+h*pi)]+k where a = amplitude, b=period, h=phase shift, and k=vertical shift I want to use following translation to bring the sine function up onto the y-axis to range from 0-1, and this

[R] Installing rgdal in R: correct -configure flags for GDAL install on Linux Redhat

2011-05-05 Thread Katrina Bennett
Hi, I'm installing rgdal but I keep having failures because I have not been able to find a good source of information for the correct configuration settings when installing GDAL. My error from the R install.packages(rgdal) is below. Can someone point me to a good source to tell me how to set

[R] zooreg and window

2011-04-14 Thread Katrina Bennett
all.indexes = end) : Incompatible methods (Ops.POSIXt, Ops.Date) for = 2: In which(in.index all.indexes = start all.indexes = end) : Incompatible methods (Ops.POSIXt, Ops.Date) for = Why is this happening? Thanks for your help. Katrina Bennett __ R

Re: [R] zooreg and window

2011-04-14 Thread Katrina Bennett
to be applied appropriately! Thank you, Katrina On Thu, Apr 14, 2011 at 12:26 AM, Achim Zeileis achim.zeil...@uibk.ac.at wrote: On Wed, 13 Apr 2011, Katrina Bennett wrote: Hello, I have a following time series data head(mend.dat) ID PARAM Year Month Day Value SYM 1 15052500     1 1965     5  15   128

[R] Timeseries Data Plotted as Monthly Boxplots

2011-02-16 Thread Katrina Bennett
Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list