[R] How fast is the JSS reviewing process?

2014-01-14 Thread Christophe Genolini
Hi all, I thought I read somewhere that the reviewing process of the Journal of Statistical Software was quite fast (arround 3 month). Is that correct? -- Christophe Genolini Maître de conférences en bio-statistique Vice président Communication interne et animation du campus Université Paris

Re: [R] How fast is the JSS reviewing process?

2014-01-14 Thread Achim Zeileis
Christophe: I thought I read somewhere that the reviewing process of the Journal of Statistical Software was quite fast (arround 3 month). Is that correct? I'm not sure whether this is the right place to ask about this... In any case, the JSS web page publishes the submit/accept/publication

Re: [R] coxph with time-dependent covariates

2014-01-14 Thread Terry Therneau
R throws out this event due to the missing time-dependent covariate at day 6. Is there a way I can keep this event without filling in a covariate value at day 6? No. The Cox model assesses the risk associated with each covariate by comparing, at each event time, the values of the subject who

[R] Problems with table

2014-01-14 Thread silvano
Hi, I have 25 questions divided into 3 modules. I'm building tables of DISCIPLINE with each of the questions (Q1, Q2, Q3, ..., Q25). Each table, of course, has a different title. In the first module (questions Q1, Q2, ..., Q7) I don’t have problems with the titles of tables. DISCIPLINA -

[R] About Assigning Missing Values

2014-01-14 Thread Chia-Chieh Lin
Hi all,   I have a data set, named ‘data1.long’, which has 180 subjects with two variables, ‘time’ and ‘patient’. The ‘time’ variable takes 9 values, i.e., 1, 2, 3,…,9, and the ‘patient’ variable takes two values, 0 and 1. There are 10 subjects for each combination of

Re: [R] About Assigning Missing Values

2014-01-14 Thread Adams, Jean
Your assignment of NAs is fine. When you extract, you have to remember to specify whether you want to include or exclude those NA values. Try this: data1.long[data1.long$time==9 !is.na(data1.long$time), time] Jean On Tue, Jan 14, 2014 at 7:22 AM, Chia-Chieh Lin

[R] about linear model in R

2014-01-14 Thread Solip Park
Hello :) I'm Solip Park in Barcelona :) I have one question about linear model in R. I used generalized linear model (glm) in R with three variables (A,B, and C). So I made a model like this; glm.model = glm (Freq ~ A * B * C, family = poisson) and then, anova (glm.model, test = Chisq). It

Re: [R] about linear model in R

2014-01-14 Thread Bert Gunter
I am not sure what you mean by expected table, but perhaps ?predict.glm is what you are looking for. If not, hopefully someone else will help. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is not knowledge. And knowledge is

Re: [R] Plot only a portion of a record

2014-01-14 Thread Greg Snow
If you use the formula method for plot then you can use the subset argument, for example: plot(Sepal.Width ~ Sepal.Length, data=iris) plot(Sepal.Width ~ Sepal.Length, data=iris, subset=Species=='setosa') On Mon, Jan 13, 2014 at 5:39 PM, Janet Choate jsc@gmail.com wrote: Hi R community, i

Re: [R] about linear model in R

2014-01-14 Thread Bert Gunter
Solip: You should usually cc R-help on your responses to increase the chance that others can help or to correct bad advice. To respond to your explanation: Read ?glm and ?family to read about how to construct glm models and family functions for them. glm models are by default fit by iterated

[R] Bivariate Binary model with conditional synchrony measures (CSM) and Odds Ratio

2014-01-14 Thread Adelino Jose Chingore Juga
Hi, Can you help me. I need to fit Bivariate Binary model with conditional synchrony measures (CSM) and Odds Ratio (OR) as measures of association. Can you give some hint code? Regards Adelino [[alternative HTML version deleted]]

[R] Running multidimensional regressions

2014-01-14 Thread andrews Nikolaiv
Dear R helpers!, I have a question on how to run a regression with many indices. To give you a practical example, let y_{itabp} be an dependent variable (representing prices) indexed by i=country, t=time, a=area, b=brand and p=package size. In particular, we collected prices on a

[R] barchar or barp on multiple data

2014-01-14 Thread Adel ESSAFI
Hello list I have the following data in file in attachment. in want to draw bars for every value of V5 in finction of V3 grouped by V2*V1. however, each unique value of V4, V2 and V1 I want to draw a bar. This example is very close to want I want to do. but here, for each value of V3, barchart

Re: [R] R-help Digest, Vol 131, Issue 15

2014-01-14 Thread Vasco Cadavez
Dear R Users, I'm searching for a Self Start function for Weibull function (3 parameters). Somebody can help? Cheers, Vasco Cadavez On 14/01/14 11:00, r-help-requ...@r-project.org wrote: Send R-help mailing list submissions to r-help@r-project.org To subscribe or unsubscribe

[R] Subsetting on multiple criteria (AND condition) in R

2014-01-14 Thread Jeff Johnson
I'm running the following to get what I would expect is a subset of countries that are not equal to US AND COUNTRY is not in one of my validcountries values. non_us - subset(mydf, (COUNTRY %in% validcountries) COUNTRY != US, select = COUNTRY, na.rm=TRUE) however, when I then do table(non_us) I

[R] Different output for lm Mac vs PC

2014-01-14 Thread CASENHISER, DEVIN M
I've noticed that I get different output when running a linear model on my Mac versus on my PC. Same effect, but the Mac assumes the predictor as a 0 level whereas the PC uses the first category (alphabetically). So for example (using Bodo Winter's example from his online linear models

Re: [R] barchar or barp on multiple data

2014-01-14 Thread John Kane
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example No data. The list strips off most files. If you want to send data give the file a .txt suffix but it is much better to supply the output of dput() . John Kane Kingston ON Canada -Original

Re: [R] Subsetting on multiple criteria (AND condition) in R

2014-01-14 Thread arun
Hi, Try: table(as.character(non_us[,COUNTRY])) A.K. On Tuesday, January 14, 2014 3:17 PM, Jeff Johnson mrjeffto...@gmail.com wrote: I'm running the following to get what I would expect is a subset of countries that are not equal to US AND COUNTRY is not in one of my validcountries values.

Re: [R] Subsetting on multiple criteria (AND condition) in R

2014-01-14 Thread Marc Schwartz
On Jan 14, 2014, at 1:38 PM, Jeff Johnson mrjeffto...@gmail.com wrote: I'm running the following to get what I would expect is a subset of countries that are not equal to US AND COUNTRY is not in one of my validcountries values. non_us - subset(mydf, (COUNTRY %in% validcountries) COUNTRY

Re: [R] Barplot not showing all labels

2014-01-14 Thread John Kane
I am not sure that I got the data correctly--it is much better to supply sample data using dput(). See ?dput for more information but I think something like this will work dat1 / - structure(list(cty = structure(1:70, .Label = c(AE, AN, AR, AT, AU, BB, BD, BE, BH, BM, BN, BO, BR, BS, CA, CH,

Re: [R] Different output for lm Mac vs PC

2014-01-14 Thread Marc Schwartz
On Jan 14, 2014, at 2:23 PM, CASENHISER, DEVIN M de...@uthsc.edu wrote: I've noticed that I get different output when running a linear model on my Mac versus on my PC. Same effect, but the Mac assumes the predictor as a 0 level whereas the PC uses the first category (alphabetically). So

Re: [R] Subsetting on multiple criteria (AND condition) in R

2014-01-14 Thread William Dunlap
Here is a reproducible example of your problem where you do not want to see a table entry for Medium. tmp_df - data.frame(Size=factor(rep(c(Small,Medium,Large),1:3), levels=c(Small,Medium,Large))) non_medium - subset(tmp_df, Size != Medium, select=Size) table(non_medium) non_medium

Re: [R] Different output for lm Mac vs PC

2014-01-14 Thread Greg Snow
I would suggest running the code: options('contrasts') on both machines to see if there is a difference. Having the default contrasts set differently would be one explanation. On Tue, Jan 14, 2014 at 2:28 PM, Marc Schwartz marc_schwa...@me.com wrote: On Jan 14, 2014, at 2:23 PM, CASENHISER,

Re: [R] Different output for lm Mac vs PC

2014-01-14 Thread Duncan Mackay
I think this is a case of using dput to make sure that the datasets are the same. I was bitten once by this and now dput between machines. Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England Armidale NSW 2351 Email: home: mac...@northnet.com.au

Re: [R] barchar or barp on multiple data

2014-01-14 Thread Adel ESSAFI
Ok, I put here the data exported with dput(). structure(list(V1 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,

Re: [R] Different output for lm Mac vs PC

2014-01-14 Thread Marc Schwartz
Good catch Greg. The Mac output observed can result from either: options(contrasts = c(contr.helmert, contr.poly)) or options(contrasts = c(contr.sum, contr.poly)) being run first, before calling the model code. I checked the referenced tutorial and did not see any steps pertaining to

Re: [R] R-help Digest, Vol 131, Issue 15

2014-01-14 Thread Rui Barradas
Hello, Maybe function ?SSweibull in package stats. Hope this helps, Rui Barradas Em 14-01-2014 18:29, Vasco Cadavez escreveu: Dear R Users, I'm searching for a Self Start function for Weibull function (3 parameters). Somebody can help? Cheers, Vasco Cadavez On 14/01/14 11:00,

Re: [R] Barplot not showing all labels

2014-01-14 Thread Jeff Johnson
Thanks John (and everyone else as well). John's example got it very close. I can tweak from here. Thanks! On Tue, Jan 14, 2014 at 1:22 PM, John Kane jrkrid...@inbox.com wrote: I am not sure that I got the data correctly--it is much better to supply sample data using dput(). See ?dput for more

Re: [R] Subsetting on multiple criteria (AND condition) in R

2014-01-14 Thread Jeff Johnson
Thanks so much Marc and for those that responded. Mark's suggestion with droplevels gave me the desired result. I'm new to figuring out how to post reproducible code. I'll try using the set.seed and rnorm functions next time and hope that does the trick. Thanks everyone! On Tue, Jan 14, 2014

[R] rJava works on R-32bit but fails in R 64bit

2014-01-14 Thread Hui Du
Hi All, I have R 64bit and R-32 bit installed in my windows 7. For 64 bit, the version info is R.Version() $platform [1] x86_64-w64-mingw32 $arch [1] x86_64 $os [1] mingw32 $system [1] x86_64, mingw32 $status [1] $major [1] 3 $minor [1] 0.2 $year [1] 2013 $month [1] 09 $day [1] 25

[R] Why is kmeans incredibly slow on this Ubuntu machine?

2014-01-14 Thread Ben Harrison
X is a 5000 x 5 matrix (numerical variables) k = 6 My desktop machine is a Xeon-processor Dell Precision T3500. [Ubuntu 12.04.4 LTS (GNU/Linux 3.2.0-58-generic x86_64) with 24 GB RAM] [R version 3.0.2 (2013-09-25) -- Frisbee Sailing Copyright (C) 2013 The R Foundation for Statistical Computing

Re: [R] R-help Digest, Vol 131, Issue 15

2014-01-14 Thread David Winsemius
On Jan 14, 2014, at 10:29 AM, Vasco Cadavez wrote: Dear R Users, I'm searching for a Self Start function for Weibull function (3 parameters). Have you looked at: ?selfStart # be sure to look at the links in the See Also section -- David Somebody can help? Cheers, Vasco

Re: [R] rJava works on R-32bit but fails in R 64bit

2014-01-14 Thread Jeff Newmiller
Post plain text per the posting guide? Install the 64bit version of the Java Runtime? --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#.

[R] (no subject)

2014-01-14 Thread Keniajin Wambui
I am using R3.02 on WIndows 7 64bit. How can I test for trend across ordered groups. I am using `prop.trend.test` but it doesn't order into groups and doesn't give the expected output of `sum of ranks` , `z value`... I have this data set and am doing `prop.trend.test (birthweight2$lbw2

[R] Implement a Fomula

2014-01-14 Thread Norman McBride
Now, if I were to move the data over from an excel doc would that change anything. I do understand that that I would have to change it to a csv document and create a variable for it. Just to give you a more visual understanding of what I am doing I attached a picture of the head of the data. So

[R] Assignment Week 1

2014-01-14 Thread Dorcas Nduati
Hi, I don't want to cheat for Assignment Week 1 I calculated the answers using microsoft excel. Please help me with extracting data, counting missing values in R software will really be grateful. Don't want to be a liar. -- Dorcas Nduati PHPT Dept. Room 111 UoN Office:0202453621 Ext 27107

[R] Model averaging using QAICc

2014-01-14 Thread Diana Virkki
Hi all, I am having some trouble running GLMM's and using model averaging with QAICc. Let me know if you need more detail here: I am trying to run GLMM's on count data in the package glmmADMB with a negative binomial distribution due to overdispersion. The dispersion parameter has now reduced to

Re: [R] Assignment Week 1

2014-01-14 Thread Koushik Saha
show us some of your effort? what have you tried? did u went through the video lectures in the site? On Wed, Jan 15, 2014 at 3:11 AM, Dorcas Nduati dndu...@gmail.com wrote: Hi, I don't want to cheat for Assignment Week 1 I calculated the answers using microsoft excel. Please help me with

Re: [R] Implement a Fomula

2014-01-14 Thread PIKAL Petr
Hi It seems to me that you would do yourself a favor if you go through R-Intro. See in line. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Norman McBride Sent: Tuesday, January 14, 2014 11:56 PM To: r-help@r-project.org

Re: [R] Assignment Week 1

2014-01-14 Thread vikram ranga
Hi Dorcas Nduati, Better place to discuss and ask for help about home work is discussion forum of that site. It is one of recommendations is that: you should not post your homework on Rhelp or stackOverflow. On Wed, Jan 15, 2014 at 12:56 PM, arun smartpink...@yahoo.com wrote: Hi,May be this