Re: [R] Input and output of time series data - any function or packages that helps?

I think this is a very vague question. If you are exporting to some other software, csv is probably a very good choice, unless there is something more specific about that software. As for "manual", I think we left the scribes in the middle ages. R is eminently scriptable for automating tasks, b

Re: [R] Input and output of time series data - any function or packages that helps?

On Oct 1, 2012, at 4:38 AM, jpm miao wrote: > Hello, > > I work with time series data. From time to time I run programs to > produce results that are in time series form (e.g., quarterly or monthly > data). After a few days I might need to access part of the results and to > run another prog

[R] Input and output of time series data - any function or packages that helps?

Hello, I work with time series data. From time to time I run programs to produce results that are in time series form (e.g., quarterly or monthly data). After a few days I might need to access part of the results and to run another program. Is there any function or package (like dataframe or zo

[R] nlme: spatial autocorrelation on a sphere

I have spatial data on a sphere (the Earth) for which I would like to run an gls model assuming that the errors are autcorrelated, i.e. including a corSpatial correlation in the model specification. In this case the distance metric should be calculated on the sphere, therefore metric = "euclid

Re: [R] Unexpected behavior with weights in binomial glm()

On Sep 30, 2012, at 4:47 AM, Josh Browning wrote: > Hi David, > > Yes, I agree that the results are "very similar" but I don't > understand why they are not exactly equal given that the data sets are > identical. > > And yes, this 1% numerical difference is hugely important to me. I believe yo

Re: [R] Unexpected behavior with weights in binomial glm()

Hi Josh, I would not trust either result in this case. Additional numerical precision is not likely to help (well, perhaps agreement between the two sets of results, but agreement doesn't help if neither results are meaningful). If this were my data, I would try an approach like this: # load re

Re: [R] Please help on Sweave question

On 12-09-30 2:18 PM, Andras Farkas wrote: Dear All, I have reviewed some examples over the net on this issue, but still not getting the results. I have the following text and code i would like to place into a pdf with sweave, but I would like to change the margins on the produced document. T

Re: [R] Unexpected behavior with weights in binomial glm()

Bert Gunter gene.com> writes: > > I haven't followed this thread closely, but if perfect separation in a > binomial glm is the problem, google it. e.g. > > http://www.ats.ucla.edu/stat/mult_pkg/faq/general >/complete_separation_logit_models.htm > > This presumably explains your concerns ab

Re: [R] Strategy for building a robust rare event logistic model - Questions.

As this does not appear to be an R question. I suggest you post on a statistics list, like stats.stackexchange.com instead. -- Bert On Sun, Sep 30, 2012 at 8:22 AM, Mehul Makwana wrote: > Hello, > > I am reaching out for a problem I have seem to hit wall on. I am trying to > build a rare event

[R] Please help on Sweave question

Dear All,     I have reviewed some examples over the net on this issue, but still not getting the results. I have the following text and code i would like to place into a pdf with sweave, but I would like to change the margins on the produced document. The options(width=60) I guess should do thi

Re: [R] Anova and tukey-grouping

Hello again ! It worked exactly the same way as yours! I'm kind of astonished, because in "my" opinion our skripts are the same (obviously not...). Whatever! It works :-D Tomorrow (as where I live it's already evening...sunday) will be a great day, indeed, and I can go on with analysis ! I really

Re: [R] Generating by inverting function

On Sep 30, 2012, at 5:34 AM, R. Michael Weylandt wrote: > On Sun, Sep 30, 2012 at 10:51 AM, Salma Wafi wrote: >> Dear All, >> Thank you very much for your kind and response. >> >> Actually, S(t) = exp(b* F(t)) represent one of survival models, where F(t) >> is known cumulative distribtion func

Re: [R] Anova and tukey-grouping

typohne.csv Hello arun kirshna, I checked str() but I do not see any mistake. Find attached a subset of my data (did you mean this with "dput()" ?) Typ= treatment abun= abundance div= diversity, number of species best wishes. -- View

[R] Strategy for building a robust rare event logistic model - Questions.

Hello, I am reaching out for a problem I have seem to hit wall on. I am trying to build a rare event model:- Default event occurred in NY/NJ area. Some key statistics: 138 default aka 4bp (.04%). Obviously my logistic and robust logistic models are failing Hosmer GOF tests and so I am thinking

Re: [R] Anova and tukey-grouping

Hi, I got your data from Nabble.  I meant to do it by dput() because the .csv files may not get through the R-help mailing list. This is what I got: dat1<-read.delim("typohne.csv")  str(dat1) #'data.frame':    60 obs. of  3 variables: # $ typ : Factor w/ 5 levels "A","B","C","D",..: 1 2 3 4 5 1

Re: [R] Unexpected behavior with weights in binomial glm()

I haven't followed this thread closely, but if perfect separation in a binomial glm is the problem, google it. e.g. http://www.ats.ucla.edu/stat/mult_pkg/faq/general/complete_separation_logit_models.htm This presumably explains your concerns about coefficient agreement. -- Bert On Sun, Sep 30,

Re: [R] Merging multiple columns into one column

Hi Meredith, Why not just use paste()? cuyahoga_nf$newcolumn <- paste(cuyahoga_nf[,1], cuyahoga_nf[,2], cuyahoga_nf[,3], . , cuyahoga_nf[,4]) Best, Steve Politzer_ahles - Original Message - > From: Meredith Ballard LaBeau > To: r-help@r-project.org > Cc: > Sent: Friday, September

Re: [R] norm tables

I don't quite understand what you are asking. It sounds more like a psychometrics or test construction question than an R question. Do you have the data and just need to analyse it? You might want to have a look at CRAN Task View: Psychometric Models and Methods for some of the available

Re: [R] modified Ward method

On 30.09.2012 04:28, eliza botto wrote: Dear useRs,Is there a command or package in R to do heirarchical clustering by "modified ward method"??regardseliza No idea - maybe someone can help given you explain us which modification of the ward method is m

Re: [R] Generating by inverting function

On Sun, Sep 30, 2012 at 10:51 AM, Salma Wafi wrote: > Dear All, > Thank you very much for your kind and response. > > Actually, S(t) = exp(b* F(t)) represent one of survival models, where F(t) is > known cumulative distribtion function. So, that is why I need to generate the > survival times T b

Re: [R] Unexpected behavior with weights in binomial glm()

Hi David, Yes, I agree that the results are "very similar" but I don't understand why they are not exactly equal given that the data sets are identical. And yes, this 1% numerical difference is hugely important to me. I have another data set (much larger than this toy example) that works on the

Re: [R] Generating by inverting function

Dear All, Thank you very much for your kind and response.   Actually, S(t) = exp(b* F(t)) represent one of survival models, where F(t) is known cumulative distribtion function. So, that is why I need to generate the survival times T based on this model. Then for generating time I have to invert

Re: [R] plot selected scales on +/- axes

On 09/29/2012 11:33 PM, Sabine Dietmann wrote: Hi All, I wonder whether it is possible to use the plot function and reduce the scales on both +/- quadrants, i.e. plot the data from -15. -8 and 8..15 only. Hi Sabine, As David noted, you might want gap.plot: library(plotrix) testdat<-c(sample(-

Re: [R] Generating by inverting function

On Sep 29, 2012, at 10:23 AM, Salma Wafi wrote: > Hello, > I am trying to generate random survival times by inverting the function, > S(t)= exp(b*F(t)), This looks a bit confused. S(t) = 1- exp( -integral( h(t)dt ) where h(t) is the instantaneous hazard. That integral is the cumulative hazar

Re: [R] Generating by inverting function

On Sat, Sep 29, 2012 at 6:23 PM, Salma Wafi wrote: > Hello, > I am trying to generate random survival times by inverting the function, > S(t)= exp(b*F(t)), where b is constant and F(t) is some cumulative > distribution function, let say that F(t) is cdf of normal distribution or any > others d

Re: [R] Query

On Sun, Sep 30, 2012 at 7:15 AM, ghoshm ghosh wrote: > Dear All, > > > If possible please let me know whether I can find solution of following > problem using R or not? > > > I have to select 6-10 parameters whose range are given but I need to select > the parameter values which satisfies 5-6 know

[R] Query

Dear All, If possible please let me know whether I can find solution of following problem using R or not? I have to select 6-10 parameters whose range are given but I need to select the parameter values which satisfies 5-6 known conditions i.e. to find a suitable parameters satisfying given co

[R] Generating by inverting function

Hello, I am trying to generate random survival times by inverting the function,  S(t)= exp(b*F(t)), where b is constant and F(t) is some cumulative distribution function, let say that F(t) is cdf of normal distribution or any others distributions.   as we know that S(t) has uniform distribution

[R] norm tables

Dear all, Is anybody knows simple method for creation norm tables for psychological tests for several groups distinguished by sex and age? best regards, Krzysztof __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS