Re: [R] Second-order effect in Parametric Survival Analysis

2011-11-12 Thread ryusuke
Thank you Dr. David. I try to summarize it. Assumes x and z are two covariates: x = dummy variable (1 or 0) z = factors (people name) x*z = x + z + x*z therefore this is not a 2nd-order interactions, it should be (for an exponential survival regression):- h(t|(X=x,Z=z)) = exp(Beta0 + XZBeta1) #-

Re: [R] identify duplicate from more than one column

2011-11-12 Thread Joshua Wiley
Hi Carlos, Here is one option: ## read in your data dat <- read.table(textConnection(" obs unithome zsex age 1 015029 18 11 053 2 015029 18 12 049 3 015029 01 11 038 4

[R] identify duplicate from more than one column

2011-11-12 Thread jour4life
Hi all, I've searched everywhere to try to find out how to do this and have had no luck. I am trying to construct identifiers for couples in a dataset. Essentially, I want to identify couples using more than one column as identifiers. Take for instance: obs unithome zsex

[R] Myriam Saavedra M. Sc. Questions about maximun radius distance

2011-11-12 Thread Myriam
Dear Mr. Baddeley   I just graduated from a Masters in Applied Mathematics on Jun19th. My thesis was about spatial distribution /a nalysis of some trees in a part of  the Congo Basic Forest. In my thesis I used your spatial package in R, and today I'm doing a more deeper study about how we chois

Re: [R] dev.new() within a loop

2011-11-12 Thread David Winsemius
On Nov 12, 2011, at 6:04 PM, Giovanni Azua wrote: Hello, I have a loop where I iterate performance data files within a folder, parse and plot them in one shot (see below). However, when executing plot_raw which invokes dev.new(..) all windows come out blank whereas if I execute each file

Re: [R] With an example - Re: rbind.data.frame drops attributes for factor variables

2011-11-12 Thread David Winsemius
On Nov 12, 2011, at 6:40 PM, Sammy Zee wrote: Thanks David. Besides rbind(), is there any other way to add a row to a data frame so that I do not lose the custom attributes. I have already told you the method that I know of. You don't seem to have taken my poin that it is not a data.frame

[R] Simulation over data repeatedly for four loops

2011-11-12 Thread Francesca
Dear Contributors, I am trying to perform a simulation over sample data, but I need to reproduce the same simulation over 4 groups of data. My ability with for loop is null, in particular related to dimensions as I always get, no matter what I try, "number of items to replace is not a multiple

Re: [R] efpFunctional construction (strucchange package)

2011-11-12 Thread bonda
Unfortunately, I couldn't find neither this source file, nor function wmax(). Was it any old version of strucchange? There are only files strucchange.R, strucchange.rdb and strucchange.rdx in directory \strucchange\R. Best, Julia -- View this message in context: http://r.789695.n4.nabble.com/efp

[R] dev.new() within a loop

2011-11-12 Thread Giovanni Azua
Hello, I have a loop where I iterate performance data files within a folder, parse and plot them in one shot (see below). However, when executing plot_raw which invokes dev.new(..) all windows come out blank whereas if I execute each file outside of a loop then I can see the plots properly. Wh

Re: [R] With an example - Re: rbind.data.frame drops attributes for factor variables

2011-11-12 Thread Sammy Zee
Thanks David. Besides rbind(), is there any other way to add a row to a data frame so that I do not lose the custom attributes. Thanks, Sammy On Sat, Nov 12, 2011 at 5:17 PM, David Winsemius wrote: > > On Nov 12, 2011, at 2:47 PM, Sammy Zee wrote: > > When I use rbind() or rbind.data.frame() to

Re: [R] Second-order effect in Parametric Survival Analysis

2011-11-12 Thread David Winsemius
On Nov 12, 2011, at 7:37 AM, ryusuke wrote: Hi experts, http://r.789695.n4.nabble.com/file/n4034318/Parametric_survival_analysis_2nd-order_efffect.JPG Parametric_survival_analysis_2nd-order_efffect.JPG As we know a normal survival regression is the equation (1) Well, I'ld like to modify it to

Re: [R] With an example - Re: rbind.data.frame drops attributes for factor variables

2011-11-12 Thread David Winsemius
On Nov 12, 2011, at 2:47 PM, Sammy Zee wrote: When I use rbind() or rbind.data.frame() to add a row to an existing dataframe, it appears that attributes for the column of type "factor" are dropped. See the sample example below to reproduce the problem. Please suggest How I can fix this.

Re: [R] changelog for MASS?

2011-11-12 Thread R. Michael Weylandt
It's on CRAN, where you can also see the record of updates: http://cran.r-project.org/web/packages/MASS/index.html Specifically, http://cran.r-project.org/web/packages/MASS/NEWS Hope this helps On Sat, Nov 12, 2011 at 8:55 AM, Xu Wang wrote: > Does anyone know where I can find a changelog for M

Re: [R] list.dir() function

2011-11-12 Thread Scott
firstly - it looks like you define orgDir but then your code later refers to indexDir? that being said, I don't think that is relevant. You could always try upgrading to the latest version and see if the problem goes away (if you're not already). I can't replicate the issue here but we can't gues

[R] changelog for MASS?

2011-11-12 Thread Xu Wang
Does anyone know where I can find a changelog for MASS? It's difficult to know whether I should ask my company to update the package or not. We are usually required to show the changelog, but I can't find one. Thank you, Xu -- View this message in context: http://r.789695.n4.nabble.com/changel

Re: [R] R v2.13.2 - Cannot find Rcmd on path?

2011-11-12 Thread jack306
Dear all: I could not able to find rcmd using the following path: either R version (2.13.2 or 2.14.0); either rtools version (2.13 or 2.14). My os is winxp. The variables set in path are as follows. c:\Rtools\bin;c:\Rtools\MinGW\bin;C:\Perl\bin; C:\Perl\site\bin; c:\Rtools\bin; C:\Program Files

[R] With an example - Re: rbind.data.frame drops attributes for factor variables

2011-11-12 Thread Sammy Zee
When I use rbind() or rbind.data.frame() to add a row to an existing dataframe, it appears that attributes for the column of type "factor" are dropped. See the sample example below to reproduce the problem. Please suggest How I can fix this. Thanks, Sammy a=c("Male", "Male", "Female", "Male") b=c

Re: [R] Passing parameters

2011-11-12 Thread ftonini
Hi Annie, The error is in the R script...you start reading with Args[1] but instead you should start with Args[5] etc. because (I believe) that R --slave --vanilla --args count as [1] [2][3][4] Therefore, Args[5] it's the first argument you use...

Re: [R] plotting a bivariate quadratic curve

2011-11-12 Thread Barry Rowlingson
On Sat, Nov 12, 2011 at 10:50 AM, Yackov Lubarsky wrote: > Hi, > > I would like to plot a x,y curve described by the equation : > > Ax^2 + Bx + Cy^2 + Dy + E == 0 > > (A,B,C,D,E are constants) > > This sounds like quite the common task but haven't been able to figure > out how to do this. Could yo

Re: [R] plotting a bivariate quadratic curve

2011-11-12 Thread Duncan Murdoch
On 11-11-12 5:50 AM, Yackov Lubarsky wrote: Hi, I would like to plot a x,y curve described by the equation : Ax^2 + Bx + Cy^2 + Dy + E == 0 (A,B,C,D,E are constants) This sounds like quite the common task but haven't been able to figure out how to do this. Could you please help ? I am new to

Re: [R] logistic Regression using lmer

2011-11-12 Thread Ben Bolker
arunkumar gmail.com> writes: > > Hi > > can we perform logistic regression using lmer function. Please help me with > this > Yes. library(lme4) glmer([reponse]~[fixed effects (covariates)]+(1|[grouping variable]), data=[data frame], family=binomial) Further questions on this

[R] State space model

2011-11-12 Thread Kristian Lind
Hi, I'm trying to estimate the parameters of a state space model of the following form measurement eq: z_t = a + b*y_t + eps_t transition eq y_t+h = (I -exp(-hL))theta + exp(-hL)y_t+ eta_{t+h}. The problem is that the distribution of the innovations of the transition equation depend on the pr

Re: [R] Please Help

2011-11-12 Thread David Winsemius
Two problems with your post: a) you used a format that is not accepted by the mail server and resulted in all of your linefeeds being removed, so it is now basically unreadable. As the Posting Guide says, please use plain text. (And please read the Posting Guide for yourself ... and do it

[R] processing data of hierarchical river network

2011-11-12 Thread tlange
Hi, I have some hydrochemical data from observation stations of a river network. The data set provides information about the river hierarchy, i.e. the order of the rivers (main river is 1st order, contributing is 2nd order, and contributing to 2nd is 3rd order... etc). Now I would like to use the

Re: [R] Subsetting data leads to funky plots

2011-11-12 Thread Sarah Goslee
Hi, On Sat, Nov 12, 2011 at 2:04 AM, Vinny Moriarty wrote: > I'm trying out a basic plot, but something about the way I subset my data > leads to problems with the plot. > > > Here is the first bit of my data set > > year,date,location,quadrat_juvenile,photo_location,photo_exists,genus,count,dive

[R] Please Help

2011-11-12 Thread avish kamdar
HiI want to construct a logliikelood function in RHere is the situationy=number of particles emitted in 1 hr period~pois(30)p=probability of detection of radiation particlesx=number of particles detected by a radiation detector~pois(30p)where p~beta(a,1)I have to calculate the loglikehood for a

[R] logistic Regression using lmer

2011-11-12 Thread arunkumar1111
Hi can we perform logistic regression using lmer function. Please help me with this -- View this message in context: http://r.789695.n4.nabble.com/logistic-Regression-using-lmer-tp4034056p4034056.html Sent from the R help mailing list archive at Nabble.com. _

[R] Second-order effect in Parametric Survival Analysis

2011-11-12 Thread ryusuke
Hi experts, http://r.789695.n4.nabble.com/file/n4034318/Parametric_survival_analysis_2nd-order_efffect.JPG Parametric_survival_analysis_2nd-order_efffect.JPG As we know a normal survival regression is the equation (1) Well, I'ld like to modify it to be 2nd-order interaction model as shown in equa

[R] plotting a bivariate quadratic curve

2011-11-12 Thread Yackov Lubarsky
Hi, I would like to plot a x,y curve described by the equation : Ax^2 + Bx + Cy^2 + Dy + E == 0 (A,B,C,D,E are constants) This sounds like quite the common task but haven't been able to figure out how to do this. Could you please help ? I am new to R so probably missing something basic. Thanks

Re: [R] need help

2011-11-12 Thread S Ellison
Supreet kaur [ksupre...@gmail.com] > how do I calculate the reliability between the two groups > using the ICCs? Looking for the irr (inter-rater reliability package) might help; as the name suggests, it includes a variety of measures of rater reliabilit

Re: [R] Estimating IRT models by using nlme() function

2011-11-12 Thread Cengiz Zopluoğlu
Thanks for your suggestion, Dimitris. I know (and use) all of these R packages for IRT (I am actually writer of one of them, EstCRM). I also use specialized commercial software such as BILOG and MULTILOG. My goal is to be able to estimate IRT model parameters in a nonlinear mixed effects regressi

Re: [R] Help

2011-11-12 Thread Patrick Burns
You've been told how to do what you ask. But I'm not convinced that you really want to do what you asked. It might be better to do whatever you want with the data leaving it all in one object. There are many ways of doing that, the 'by' function is one of them. On 11/11/2011 20:24, Francesca wr

Re: [R] Estimating IRT models by using nlme() function

2011-11-12 Thread Dimitris Rizopoulos
Note that there are several IRT dedicated packages in R, such as the ltm and eRm packages. For more info you may check the Psychometrics Task View at: http://cran.r-project.org/web/views/Psychometrics.html I hope it helps. Best, Dimitris On 11/11/2011 9:38 PM, Cengiz Zopluoğlu wrote: Hi,

Re: [R] Win upgrade pb (virus)

2011-11-12 Thread Mario Valle
Thanks to all that responded privately! If you have AVG and needs the colorspace (or ggplot2) package do the following from a suitably privileged account: 1) Open AVG gui and select Components->Anti virus from the menu 2) Click "Advanced settings..." and then go to "Temporary disable AVG prote