Re: [R] Bootstrapping in R

Hi, I have read the help page and it was helpful but, I am having concerns because each time I run this code I get the same value. I expected that each time I run the code, I will get different values due to random sampling. How do I get this randomization? The values shouldn’t be the same eac

[R] transform

Hi all, I want to standardize a variable based on certain condition Here is my sample of data obs, Year, bb, kk, y 1, 2001, 25 ,100, 12.6 2, 2001, 15 ,111, 24.7 3, 2001, 53, 110, 13.8 4, 2001, 50, 75, 9.6 5, 2001, 125, 101, 31.5 6, 2001, 205, 407, 65.7 7, 2001, 250, 75, 69.1 If the

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

Hi Deb: > > gribfile <- > > 'http://thredds.ucar.edu/thredds/ncss/grib/NCEP/GFS/Global_0p5deg/best?north=47.0126&west=-114.841&east=-112.641&south=44.8534&time_start=present&time_duration=PT3H&accept=netcdf&var=v-component_of_wind_height_above_ground,u-component_of_wind_height_above_ground' > > d

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

Hi Now I am using netcdfSubset and I am able to download the file but not sure how to read the files. here my scripts library("ncdf4") gribfile<-" http://thredds.ucar.edu/thredds/ncss/grib/NCEP/GFS/Pacific_40km/best/dataset.html " download.file(gribfile,basename(gribfile),mode = "wb") x<-nc_open(g

[R] Optimization issue - Solution converges for any initial value

Narendra You should check the definition of your objective function, something seems to be wrong. When you evaluate the objective function with your initial guess, and then with some different value, you should expect the returned value will be different, but it is not: > print(Koval.Error.

Re: [R] How to test a difference in ratios of count data in R

> On Sep 30, 2016, at 9:40 AM, Shuhua Zhan wrote: > > > Thank you, David and Greg for your help! > I drew conclusion that the treatment B significantly increases the ratio of x > group (X/n) from based on p values from the treatmentB line of the outputs > at logistic reg. and Poisson reg.(p=

Re: [R] Elegant way to get specific dates within prespecified period

I would probably try a different strategy. First, construct a sequence of January 1 dates that is a little bit too long. For example, start with the January in the same year as "open" and finish with the January in the same year as "close". You can construct the vector using seq(), like this: > s

Re: [R] How to test a difference in ratios of count data in R

Thank you, David and Greg for your help! I drew conclusion that the treatment B significantly increases the ratio of x group (X/n) from based on p values from the treatmentB line of the outputs at logistic reg. and Poisson reg.(p=6.11e-07, Logistic; p=0.000152, Poisson). I'm wondering whether

Re: [R] How to implement "zero-overhead" code re-use (a la Python, Perl, etc.) in R?

I believe that you should post this on the R-devel list. I think that the expertise to give you an authoritative answer is likely to reside there, and they may not monitor or wish to respond on this list. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming

[R] How to implement "zero-overhead" code re-use (a la Python, Perl, etc.) in R?

I'm collaborating in a long-running research project that, over the years, has accummulated source code (written in-house) in several languages: Python, Perl, Mathematica, MATLAB. Recently I have started writing source code in R for this project, and I am having trouble incorporating it into our e

[R] isssues with predict.coxph, offset, type = "expected", and newdata?

Hi, It seems there might be two issues with predict.coxph(), and I'd appreciate confirmation before submitting a bug report to the package author. (1) predict() seems to produce incorrect predictions when using type = "expected" from a Cox model with an offset and specifying a new data for predi