[R] Generating by inverting function

2012-09-30 Thread Salma Wafi
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

2012-09-30 Thread Krzysztof Fronczyk
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

[R] Query

2012-09-30 Thread ghoshm ghosh
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

Re: [R] Query

2012-09-30 Thread R. Michael Weylandt
On Sun, Sep 30, 2012 at 7:15 AM, ghoshm ghosh ghoshm2...@gmail.com 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

Re: [R] Generating by inverting function

2012-09-30 Thread R. Michael Weylandt
On Sat, Sep 29, 2012 at 6:23 PM, Salma Wafi salmawaf...@yahoo.com 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

Re: [R] Generating by inverting function

2012-09-30 Thread David Winsemius
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 hazard

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

2012-09-30 Thread Jim Lemon
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)

Re: [R] Generating by inverting function

2012-09-30 Thread Salma Wafi
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] Unexpected behavior with weights in binomial glm()

2012-09-30 Thread Josh Browning
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

2012-09-30 Thread R. Michael Weylandt
On Sun, Sep 30, 2012 at 10:51 AM, Salma Wafi salmawaf...@yahoo.com 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

Re: [R] modified Ward method

2012-09-30 Thread Uwe Ligges
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

Re: [R] norm tables

2012-09-30 Thread John Kane
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] Merging multiple columns into one column

2012-09-30 Thread Stephen Politzer-Ahles
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 mmbal...@mtu.edu To: r-help@r-project.org Cc: Sent:

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

2012-09-30 Thread Bert Gunter
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

Re: [R] Anova and tukey-grouping

2012-09-30 Thread arun
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 2 3 4 5

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

2012-09-30 Thread Mehul Makwana
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

Re: [R] Anova and tukey-grouping

2012-09-30 Thread Landi
typohne.csv http://r.789695.n4.nabble.com/file/n4644635/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

Re: [R] Generating by inverting function

2012-09-30 Thread David Winsemius
On Sep 30, 2012, at 5:34 AM, R. Michael Weylandt wrote: On Sun, Sep 30, 2012 at 10:51 AM, Salma Wafi salmawaf...@yahoo.com 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

Re: [R] Anova and tukey-grouping

2012-09-30 Thread Landi
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

[R] Please help on Sweave question

2012-09-30 Thread Andras Farkas
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

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

2012-09-30 Thread Bert Gunter
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 makwan...@yahoo.com wrote: Hello, I am reaching out for a problem I have seem to hit wall on. I am trying to

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

2012-09-30 Thread Ben Bolker
Bert Gunter gunter.berton at 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

Re: [R] Please help on Sweave question

2012-09-30 Thread Duncan Murdoch
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.

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

2012-09-30 Thread Joshua Wiley
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: [R] Unexpected behavior with weights in binomial glm()

2012-09-30 Thread David Winsemius
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 you will

[R] nlme: spatial autocorrelation on a sphere

2012-09-30 Thread Dan Bebber
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 =

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

2012-09-30 Thread jpm miao
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