Re: [R] Question with uniroot function

2015-04-16 Thread li li
Thank you. 2015-04-16 12:33 GMT-04:00 William Dunlap wdun...@tibco.com: Use optimize() to find the minimum and feed that value into uniroot(). Bill Dunlap TIBCO Software wdunlap tibco.com On Thu, Apr 16, 2015 at 7:47 AM, li li hannah@gmail.com wrote: Hi Jeff, Thanks for the

[R] Question with uniroot function

2015-04-15 Thread li li
Hi all, In the following code, I am trying to use uniroot function to solve for the root (a and b in code below) for function f1. I am not sure why uniroot function does not give the answer since when we look the graph, the function does cross 0 twice. Any suggestion? Thanks. Hanna

Re: [R] Question with uniroot function

2015-04-15 Thread Jeff Newmiller
You really need to read the help page for uniroot. The sign needs to be different at the ends of the starting interval. This is a typical limitation of numerical root finders. --- Jeff NewmillerThe

Re: [R] question

2015-04-13 Thread Boris Steipe
That all sounds so straightforward I wonder why you don't just code it up and try it out. You might profit from the advice of chapter 2 of the R-Inferno for your main problem. If some people who you ask for advice think that arrays in R are intrinsically slow, you might also want to look for

[R] question

2015-04-13 Thread Mahdiyeh Erfaniyan
Hi all, Recently I sent an email and I was asked to provide reproducible code of a simple example of my situation. Instead of providing the code, I decided to describe what I need in my code. I've written a function V, which is a function of (r,s); so I have a function V(r,s) in fact. The output

[R] Question on CCA and RDA analysis

2015-04-10 Thread Luis Fernando García
Dear R experts, I wanted to know if you can suggest me any website or tutorial just to learn about how to make a RDA or CDA in R Thanks in advance! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE

Re: [R] Question on CCA and RDA analysis

2015-04-10 Thread Ben Bolker
Luis Fernando García luysgarcia at gmail.com writes: Dear R experts, I wanted to know if you can suggest me any website or tutorial just to learn about how to make a RDA or CDA in R Thanks in advance! I hate to ask, but did you try Googling canonical correspondence analysis R ...

Re: [R] Question on CCA and RDA analysis

2015-04-10 Thread Luis Fernando García
Yeah, The most useful example I found was this. https://gist.github.com/perrygeo/7572735. I always had the idea of this kind of forums was to provide sources not so obvious in the web. If you have something better it would be great. 2015-04-10 18:36 GMT-03:00 Ben Bolker bbol...@gmail.com:

[R] question on waveletcomp plot image

2015-04-03 Thread Sudheer Joseph
Dear R experts, I have used waveletcomp package of R and was trying to get the dates formatted as month year but get below error while trying it with the example provided in http://www.hs-stat.com/projects/WaveletComp/WaveletComp_guided_tour.pdf Kindly help me with the

[R] question about the geometric distribution, please

2015-03-09 Thread Erin Hodgess
Hello! The function rgeom is based on the geometric distribution such that x = 0, 1, Is there a function which produces the geometric results such that x = 1, 2, ... please? Thought I'd check before I started coding. Thank you! Sincerely, Erin -- Erin Hodgess Associate Professor

Re: [R] question about the geometric distribution, please

2015-03-09 Thread David Winsemius
On Mar 9, 2015, at 1:52 PM, Erin Hodgess wrote: Hello! The function rgeom is based on the geometric distribution such that x = 0, 1, Is there a function which produces the geometric results such that x = 1, 2, ... please? Thought I'd check before I started coding. This seemed

Re: [R] question about the geometric distribution, please

2015-03-09 Thread Huang, JS
Maybe you are looking for rnbinom? JS Huang 636.536.5635 -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Erin Hodgess Sent: Monday, March 09, 2015 3:52 PM To: R help Subject: [R] question about the geometric distribution, please Hello! The function

[R] [Question] about plot.xts

2015-02-14 Thread KWCG HE
Hi All, I am trying to use xts and xtsExtra packages to plot multiple time series on one plot. I got two questions about this package. What's the meaning of *The following object is masked from ‘package:xts’:* when load xts and xtsExtra?* which plot.xts will be available if I local xts first

Re: [R] [Question] about plot.xts

2015-02-14 Thread Jeff Newmiller
It means that xtsExtra has a different version if plot.xts, and that will be the one that gets run when you call that function. If you want to use the one that is defined in xts, you can specify it... xts::plot.xts(...) In the future, don't post in HTML, since this is a plain text only list.

Re: [R] Question

2015-01-20 Thread Uwe Ligges
On 19.01.2015 23:58, Glenda Palacios wrote: I have any problems in R: The downloaded binary packages are in /var/folders/s9/kr21631n3nj30hb_4td3pv_8gn/T//RtmpDDAFmJ/downloaded_packages Warning messages: 1: In download.file(url, destfile, method, mode = wb, ...) : downloaded length

[R] Question

2015-01-19 Thread Glenda Palacios
I have any problems in R: The downloaded binary packages are in /var/folders/s9/kr21631n3nj30hb_4td3pv_8gn/T//RtmpDDAFmJ/downloaded_packages Warning messages: 1: In download.file(url, destfile, method, mode = wb, ...) : downloaded length 2122881 != reported length 2235972 2: 'tar' returned

Re: [R] Question abt 'factanal' function

2015-01-18 Thread PIKAL Petr
you mean by couldn't get any output Cheers Petr -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of ??? Sent: Sunday, January 18, 2015 9:22 AM To: r-help Subject: [R] Question abt 'factanal' function Hi, I tried to run 'factanal' function for Factor

[R] Question abt 'factanal' function

2015-01-18 Thread 오건희
Hi, I tried to run 'factanal' function for Factor Analysis, but failed to get results. The error message was following, Error in factanal(data, factors = 3, rotation = varimax) : unable to optimize from this starting value I searched for the solution of this error, and one said putting

Re: [R] Question about package principal

2015-01-12 Thread William Revelle
As David suggests, look at your data. For instance, there seems to be only 1 case (#2) for seat.width.club with non-zero data. I find it hard to believe that the other planes have seat.widths of 0! I think you probably want to code the 0s as missing, rather than 0. You also want to rethink

Re: [R] Question about package principal

2015-01-11 Thread David Winsemius
On Jan 11, 2015, at 3:55 AM, 오건희 wrote: Hi, I tried to run principal function in the 'psych' package, but it failed to do.. here is both my code and error message. I searched on the web, but couldn't find the exact answer I wanted. data-read.csv(

[R] Question about package principal

2015-01-11 Thread 오건희
Hi, I tried to run principal function in the 'psych' package, but it failed to do.. here is both my code and error message. I searched on the web, but couldn't find the exact answer I wanted. data-read.csv( https://raw.githubusercontent.com/mylesmharrison/delta_PCA_kmeans/master/delta.csv

[R] Question

2014-12-28 Thread Iskender Karagul
Dear All, I am a very fresh user of R platform. I completed the download of R for Windows, from the http://cran.r-project.org/ Web address. After opening the program I wanted to Go to Packages Tab and firstly Set CRAN mirror and then Install Package(s). However, I saw a statement in the

Re: [R] Question

2014-12-28 Thread Duncan Murdoch
On 28/12/2014 11:12 AM, Iskender Karagul wrote: Dear All, I am a very fresh user of R platform. I completed the download of R for Windows, from the http://cran.r-project.org/ Web address. After opening the program I wanted to Go to Packages Tab and firstly Set CRAN mirror and then

Re: [R] Question

2014-12-28 Thread Jeff Newmiller
This is (mostly) normal. The R software comes with an initial set of packages that are kept in the Program Files directory, which is only modifiable using administrator privileges. If you are running your machine just for you (as most people do) then I highly advise accepting the option during

[R] Question

2014-12-25 Thread Ankita Shukla
Hi I was fitting this model with two matrices. But this bifit function didnt run in my R studio and i got error message. Could anyone guide me how to install this command. I am giving that command below to have better understanding. x.f - cbind(log(Q5.all.pf), log(Q5.all.pf)^2) y.f -

[R] Question about CHAID

2014-12-17 Thread Francisco M. da Rocha
Hallo there, I would like to work with CHAID, but the newest version of R does have it. So I thought I could use an older version of R which accepts or has the library CHAID. Could you tell me which version it is and where to download it? Thanks a lot in advance. Francisco

Re: [R] Question about CHAID

2014-12-17 Thread Achim Zeileis
On Wed, 17 Dec 2014, Francisco M. da Rocha wrote: Hallo there, I would like to work with CHAID, but the newest version of R does have it. So I thought I could use an older version of R which accepts or has the library CHAID. Could you tell me which version it is and where to download it?

[R] Question

2014-12-04 Thread Ghasemzadeh-Barvarz, Massoud
Could you please let me know whether R software has any advantage over MATLAB to analyze very big datasets? Thanks you so much [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

[R] Question on LIMMA analysis with covariates and some missing data

2014-12-03 Thread Bertrand
Hello, I have a dataset of asthma patients for which white blood cells gene expression was measured with one-color Affymetrix microarrays (N~500, asthma is a factor with 4 levels: control, moderate, severe, severe smokers). I also have an extensive clinical dataset related, but with many

[R] question regarding rmvDAG in package pcalg

2014-11-25 Thread Jack Luo
Hi, I am trying to use rmvDAG in pcalg package to generate data from DAG structure. One thing I found is that when the number of variables gets large, there can be really large numbers in the data matrix. I played around with different parameters and it looks like the same case. library(pcalg)

[R] Question on R lattice graphics

2014-11-19 Thread Julie Hope
Hi, I am currently using the lattice library and trellis graphs to explore my data. Without going into too much detail of the data I just have a question about how dot plot in the lattice library scales things... I have normalised data that I am exploring as a function of time partitioned by

Re: [R] Question on R lattice graphics

2014-11-19 Thread Julie Hope
Hi, Ignore the last email - R has plotted the number of data points (100) I have and using that as a scale rather than the normalised data values Sorry to be a pest!! Julie Julie Hope (NERC PhD Student) Sediment Ecology Research Group, University of St Andrews School of Geoscience, Bangor

[R] Question about range of letters

2014-10-04 Thread Nia Gupta
Hello, I have a column with a bunch of letters. I would like to keep some of these letters (A,C,D,L) and turn the rest into 'X'. I have tried using ifelse with '|' in between the argument but it didn't work nor did 4 separate ifelse statements. Example, I currently have: Letters    A    B  

Re: [R] Question about range of letters

2014-10-04 Thread Ben Bolker
Nia Gupta nia_gupta at yahoo.com writes: Hello, I have a column with a bunch of letters. I would like to keep some of these letters (A,C,D,L) and turn the rest into 'X'. I have tried using ifelse with '|' in between the argument but it didn't work nor did 4 separate ifelse statements.

Re: [R] Question about range of letters

2014-10-04 Thread Robert Baer
On 10/4/2014 8:21 AM, Nia Gupta wrote: Hello, I have a column with a bunch of letters. I would like to keep some of these letters (A,C,D,L) and turn the rest into 'X'. I have tried using ifelse with '|' in between the argument but it didn't work nor did 4 separate ifelse statements.

[R] question regarding pcalg package

2014-09-26 Thread Jack Luo
Hi, I am trying to use the pcalg package to do some causal inference on some high dimensional omics data (~ 30k variables). It takes forever to run and my machine get stuck. I just realized that conditional independence test can be calculated without using the correlation matrix (I think it's

Re: [R] Question about searchTwitter{twitteR}

2014-09-11 Thread Cory N
You are only able to search twitter history for a short period of time. gnip.com and similar companies offer historical tweets for sale. cn On Sunday, September 7, 2014 9:21:34 AM UTC-5, Axel Urbiz wrote: Hello, The function searchTwitter() with the arguments supplied as below would

[R] Question about searchTwitter{twitteR}

2014-09-07 Thread Axel Urbiz
Hello, The function searchTwitter() with the arguments supplied as below would give me a different number of results on different days I run this code. Maybe it is my lack of understanding about what the date arguments are supposed to do in this function, but I would think I should be getting the

Re: [R] Question about searchTwitter{twitteR}

2014-09-07 Thread Amos B. Elberg
Twitter tweets aren't a stable database. I wouldn't expect the search results to stay stable, as tweets are retweeted, deleted, accounts are closed, privacy settings adjusted, etc. And if there are more than 1000 results, I don't know that twitter is internally ordered so you'd get the same

Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-29 Thread peter dalgaard
I'm no expert on hurdle models, but it seems that you are unaware that the negative binomial and the truncated negative binomial are quite different things. -pd On 29 Aug 2014, at 05:57 , Nick Livingston nlivings...@ymail.com wrote: I have sought consultation online and in person, to no

Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-29 Thread Achim Zeileis
On Fri, 29 Aug 2014, peter dalgaard wrote: I'm no expert on hurdle models, but it seems that you are unaware that the negative binomial and the truncated negative binomial are quite different things. Yes. You can replicate the truncated count part of the hurdle model with the zerotrunc()

Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-29 Thread Nick Livingston
: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS To: peter dalgaard pda...@gmail.com Cc: Nick Livingston nlivings...@ymail.com, r-help@r-project.org Date: Friday, August 29, 2014, 5:26 AM On Fri, 29 Aug 2014, peter dalgaard wrote: I'm

Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-29 Thread David Winsemius
...@uibk.ac.at wrote: Subject: Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS To: peter dalgaard pda...@gmail.com Cc: Nick Livingston nlivings...@ymail.com, r-help@r-project.org Date: Friday, August 29, 2014, 5:26 AM On Fri, 29 Aug 2014

Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-29 Thread Prof Brian Ripley
: just untick 'binary'.) Thank you again. I appreciate your input. -Nick On Fri, 8/29/14, Achim Zeileis achim.zeil...@uibk.ac.at wrote: Subject: Re: [R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

[R] Question regarding the discrepancy between count model parameter estimates between pscl and MASS

2014-08-28 Thread Nick Livingston
I have sought consultation online and in person, to no avail. I hope someone on here might have some insight. Any feedback would be most welcome. I am attempting to plot predicted values from a two-component hurdle model (logistic [suicide attempt yes/no] and negative binomial count [number of

Re: [R] Question

2014-07-30 Thread arun
Hi Farnoosh, Regarding the first question: dat2 - dat1 dat1$Mean - setNames(unsplit(sapply(split(dat1[,-1], dat1[,1]),rowMeans, na.rm=T),dat1[,1]),NULL) dat1   Unit q1 q2 q3 Mean 1    A  3  1  2 2.00 2    A  2 NA  1 1.50 3    B  2  2  4 2.67 4    B NA  2  5 3.50 5    C  3  2

Re: [R] Question regarding lattice::levelplot and distribution of colors

2014-07-07 Thread Sunny Srivastava
Science University of New England Armidale NSW 2351 Email: home: mac...@northnet.com.au -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Sunny Srivastava Sent: Monday, 7 July 2014 09:35 To: R mailing list Subject: [R

Re: [R] Question regarding lattice::levelplot and distribution of colors

2014-07-07 Thread Duncan Mackay
for certain ranges Duncan From: Sunny Srivastava [mailto:research.b...@gmail.com] Sent: Monday, 7 July 2014 18:57 To: Duncan Mackay Cc: R Subject: Re: [R] Question regarding lattice::levelplot and distribution of colors I am sorry but I think I have been unclear. The problem is not color theme

[R] Question regarding lattice::levelplot and distribution of colors

2014-07-06 Thread Sunny Srivastava
Hello R-helpers: I think there is some problem with my code, but I would like to seek you help because I can't spot it. I have a data.frame defined as follows: testdf - structure(list(yy = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c(R, L), class =

Re: [R] Question regarding lattice::levelplot and distribution of colors

2014-07-06 Thread David Winsemius
On Jul 6, 2014, at 4:35 PM, Sunny Srivastava wrote: Hello R-helpers: I think there is some problem with my code, but I would like to seek you help because I can't spot it. I have a data.frame defined as follows: testdf - structure(list(yy = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,

Re: [R] Question regarding lattice::levelplot and distribution of colors

2014-07-06 Thread Duncan Mackay
: home: mac...@northnet.com.au -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Sunny Srivastava Sent: Monday, 7 July 2014 09:35 To: R mailing list Subject: [R] Question regarding lattice::levelplot and distribution of colors Hello R

Re: [R] Question regarding lattice::levelplot and distribution of colors

2014-07-06 Thread Sunny Srivastava
Hello David: Thanks for your response. I am not sure, but isn't an ordered categorical variable integer for all practical purposes? (R L) Further, the problem persists if I change 'at' from seq(-0.3, 0.3, length = 20) to seq(-0.3, 0.3, length = 50; I think this argument is for the color key

[R] Question on JAVA and R

2014-06-16 Thread Luis Eduardo Castillo Méndez
As I discuss today we are trying to do with my partner Steven Penaloza and a small program of supervised classification of images with unconventional using R methods, connecting libraries and own code of R to an environment of JAVA, in our case we are using NETBEANS to interface, the connection

[R] question about chi values GLM

2014-06-14 Thread Luis Fernando García
Dear all, I am making an analysis using a GLM using three explanatory variables and a response variable. I need to obtain a table similar to this one, http://postimg.org/image/5sau79wlt/r nevertheless, I have not been able to do it. I am having a hard time specially getting the chi square

Re: [R] question about chi values GLM

2014-06-14 Thread peter dalgaard
The column labeled Deviance pretty much _is_ the chi-square, specifically the likelihood ratio test statistic, which has an asymptotic chi-square distribution. (Using test=Rao gives you the alternative Rao efficient score test, which in your case doesn't make much of a difference.) Notice

Re: [R] question about chi values GLM

2014-06-14 Thread Prof Brian Ripley
On 14/06/2014 09:45, peter dalgaard wrote: The column labeled Deviance pretty much _is_ the chi-square, specifically the likelihood ratio test statistic, which has an asymptotic chi-square distribution. (Using test=Rao gives you the alternative Rao efficient score test, which in your case

[R] question about how the mantel test is run in the ade4 package

2014-06-04 Thread Monaly Mistry
Hi, I'm running a mantel test of the location of individuals in space against the behavioural score. I wanted to see whether those with similar behavioural scores are closer together in space relative to those that are further apart. This is the code I have written:

Re: [R] Question about setdiff()

2014-06-03 Thread Σταῦρος Μακράκης
There are two kinds of set difference: the usual set difference https://en.wikipedia.org/wiki/Set_difference#Relative_complement, also called the asymmetric difference or relative complement, and the symmetric difference https://en.wikipedia.org/wiki/Symmetric_difference. In R, setdiff(a,b) is

[R] Question about setdiff()

2014-06-02 Thread Raphael Päbst
Hello everyone, I have a question which is probably rooted in my lack of understanding when it comes to math. I just did the following: v - c(1:20) w - c(11:30) setdiff(v, w) and got: 1 2 3 4 5 6 7 8 9 10 Then I did the following: setdiff(w, v) and got, not surprisingly: 21 22 23 24 25 26 27

Re: [R] Question about setdiff()

2014-06-02 Thread Jeff Newmiller
About time for you to adjust your expectations... looks right to me from both a mathematical sense and as the functions are designed. --- Jeff NewmillerThe . . Go Live...

Re: [R] Question about setdiff()

2014-06-02 Thread Pascal Oettli
Hello, From the help page: Performs *set* union, intersection, (asymmetric!) difference, equality and membership on two vectors. Hope this helps, Pascal On Mon, Jun 2, 2014 at 3:57 PM, Raphael Päbst raphael.pae...@gmail.com wrote: Hello everyone, I have a question which is probably rooted in

Re: [R] Question about setdiff()

2014-06-02 Thread Raphael Päbst
Thanks everyone! It is just as I expected, I just didn't understand how setdiff() works. Raphael On 6/2/14, Pascal Oettli kri...@ymail.com wrote: Hello, From the help page: Performs *set* union, intersection, (asymmetric!) difference, equality and membership on two vectors. Hope this

Re: [R] Question about setdiff()

2014-06-02 Thread arun
Hi, Please check this link: http://r.789695.n4.nabble.com/meaning-of-asymmetric-on-help-page-for-intersect-td877408.html union(setdiff(v,w), setdiff(w,v)) #or in this case setdiff(union(v,w),intersect(v,w)) #or  setdiff(c(v,w),c(v,w)[duplicated(c(v,w))]) A.K. .pae...@gmail.com wrote: Hello

[R] question about wfct in minpack.lm package

2014-05-15 Thread Andras Farkas
Dear All,   made some headways with my nls model, thanks for the help... I would like to ask if you could provide some insights to the following:   I am trying to use wfct() from minpack.lm so that I could do the weighting using one of the predictors, but getting an error message:  

Re: [R] Question on Cluster Package, agnes() function

2014-05-03 Thread Martin Maechler
Anna F... on Thu, 1 May 2014 22:09:28 + writes: Hi Martin, I am a statistician at National Jewish Health in Colorado, and I have been working on clustering a dataset using Ward's minimum variance. When plotting the dendrogram, the y-axis is labeled as 'height'. Can you

[R] Question about rgamma when the shape and the scale are vectors, not scalars

2014-04-30 Thread Stefano Sofia
Dear list users, could somebody explain to me which is the meaning of rgamma(3, a, b) when both a and b are two vectors of length 3? I know that a is the shape and b is the scale, but I thought that they could only be scalars. Thank you for your help Stefano

Re: [R] Question about rgamma when the shape and the scale are vectors, not scalars

2014-04-30 Thread Greg Snow
Your command will generate 3 random values from gamma distributions, the first will be from a gamma with shape a[1] and scale b[1], then the 2nd will come from a gamma with shape a[2] and scale b[2] and the 3rd will have shape a[3] and scale b[3]. On Wed, Apr 30, 2014 at 3:00 PM, Stefano Sofia

Re: [R] Question: Do I need to set refit=FALSE when testing for random effects with anova()?

2014-04-07 Thread Christian Brauner
The question has been answered here: http://stackoverflow.com/questions/22892063/do-i-need-to-set-refit-false-when-testing-for-random-effects-in-lmer-models-wi [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Question: Do I need to set refit=FALSE when testing for random effects with anova()?

2014-04-06 Thread Christian Brauner
Hello, I am currently testing whether I should include certain random effects in my lmer model or not. I use the anova function for that. My procedure so far is to fit the model with a function call to lmer() with REML=TRUE (the default option). Then I call anova() on the two models where one of

[R] Question about c.trellis (latticeExtra)

2014-03-26 Thread Pascal Oettli
Dear r-help mailing list readers, I am facing a problem using the c.trellis function from the latticeExtra package, to merge 2 trellis objects (1 levelplot and 1 xyplot). Using the following example, it works well without customizing the y-axis of levObj1. But when the y-axis is customized

[R] R question about intraindividual variables

2014-03-13 Thread Nastassia Hajal
I am creating variables representing intraindividual means and standard deviations for longitudinal data with the following code: data$TraitHAPPYmean - with(data, ave(Happy, ID, FUN=function(x) mean(x, na.rm=TRUE) ) ) data$TraitHAPPYsd - with(data, ave(Happy, ID, FUN=function(x) sd(x,

Re: [R] question on more efficient data block-match processing

2014-03-07 Thread Rainer Schuermann
What I would do: # read in your sample data mbr - read.table( clipboard, header = TRUE, stringsAsFactors = FALSE ) # create a vector with the codes you want to consider code.list - c(A,B,C,D,E) # reduce the data accordingly mbr - mbr[ mbr$code %in% code.list, ] # get your model matrix using

[R] Question: How Kolmogorov Test is computed by R function ks.test, using scaled difference or unscaled

2014-03-07 Thread chee chen
Dear All,   I would like to ask for your help on the following when testing uniformity by the KS test:   1.  We know that the Kolmogorov test is based on the supremum D of the difference between the empirical distribution and the uniform distribution 2.  We also know that sqrt(n)D has a limiting

[R] Question on 'get'

2014-03-06 Thread Brian Smith
Hi, I was trying to get at some values from 'data' using the get function. Here is my code: class(data) [1] data.frame class(data$gender.factor) [1] factor head(data$gender.factor) [1] Male Female Male Female Male Female Levels: Male Female xx - get(data$gender.factor) Error in

Re: [R] Question on 'get'

2014-03-06 Thread Sarah Goslee
Why do you need to use the get() command? If you want to access a column of a known data frame by name: R fakedata - data.frame(A=1:3, B=letters[8:10], C=runif(3)) R fakedata[A] A 1 1 2 2 3 3 If you're trying to access a data frame by name, then you do need get, but can then subset normally.

[R] question on more efficient data block-match processing

2014-03-06 Thread Mckinstry, Craig
I have a medical insurance claims datafile divided into blocks by member, with multiple lines per member. I am process these into a one line per member model matrix. Member block sizes vary from 1 to 50+. I am match attributes in claims data to columns in the model matrix and have been

[R] R question

2014-02-02 Thread dvf208
I have generated a L1 penalized Cox model using the penalized package in R. I used the optL1() function to generate the Breslow object (see below): fit - optL1(surv.obj, penalized = ..., etc) In the reference manual, it says the fit$predictions are the cross-validated predictions for the left

Re: [R] Question on density values obtained from kde2d() from package MASS

2014-01-23 Thread Peter Löwenberg
Hello Andy, Do you already understand kde2d output? if I run max(tmp$z) do I take the maximum value of the 2d density? Cheers, Peter Em quinta-feira, 24 de novembro de 2011 16h12min19s UTC-2, Andreas Klein escreveu: Hello, I am a little bit confused regarding the density values

[R] Question on adding a p-value in bwplot

2014-01-16 Thread Hae-Young Kim
Hi, I am using bwplot to depict the box plots for two group by 6 time points. I need to add 6 p-values in each time point to compare two group at each time point. P-values are (0.0020, 0.0204, 0.3361, 0.0185, 0.1981, and 0.6677). I could depict the two box plots per each time point using the

[R] Question regarding the package rugarch

2014-01-07 Thread Grünenfelder Corina
Dear all, I have a question regarding the package rugarch. I would like to fit a garch model with exogenous variables using rugarch. My code is as follows, where data.reg is a time series object where the first column corresponds to the response variable and the remainder columns are the

[R] Question: reproducibility of random sampling with replacement

2014-01-01 Thread Chee Chen
Dear All, I would like to ask for your help on reproducibility of random sampling with replacement. For example, one re-samples the rows with replacement of a residual matrix and uses the new residual matrix thus obtained to produce a statistic ; repeat this for a certain number of times. My

Re: [R] Question: reproducibility of random sampling with replacement

2014-01-01 Thread Rui Barradas
Hello, Inline. Em 01-01-2014 22:12, Chee Chen escreveu: Dear All, I would like to ask for your help on reproducibility of random sampling with replacement. For example, one re-samples the rows with replacement of a residual matrix and uses the new residual matrix thus obtained to produce a

Re: [R] Question: reproducibility of random sampling with replacement

2014-01-01 Thread Rolf Turner
Why on earth would you expect S and T to be the same given what you have done. I am unable to rightly apprehend the confusion of ideas that could provoke such a question, (Charles Babbage). You have to set the *same* seed before each construction. I.e. do set.seed(123) before creating S; then

Re: [R] Question: reproducibility of random sampling with replacement

2014-01-01 Thread Jeff Newmiller
If you want to reproduce the same sequence twice, then you need to set the seed at the beginning of each calculation. You are only doing it for the second calculation below. --- Jeff NewmillerThe

Re: [R] Question: reproducibility of random sampling with replacement

2014-01-01 Thread Bert Gunter
You have to set the same seed before each random number generation! You did not do this. Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is not knowledge. And knowledge is certainly not wisdom. H. Gilbert Welch On Wed, Jan 1, 2014 at 2:12

[R] Question

2013-12-16 Thread Rasool Asoodeh
Hello Dears. I have a 2 questions about Discriminant Analysis in R. 1- In ade4 package I perform this analysis for both quantitative and qualitative variables using discrimin function. R gives me Canonical weights/ Loadings and Canonical scores but doesn't give *F* or *t *or *X-square

[R] question about how to install package rjava in R 3.02

2013-12-04 Thread tlw1987 [via R]
Hello, everybody , recently , I want to install the package rjava in windows 7 64bt. But it can not success.This error messages is as follow: install.packages(rjava) Warning in install.packages : package ‘rjava’ is not available (for R version 3.0.2) Warning in install.packages : Perhaps

Re: [R] question about how to install package rjava in R 3.02

2013-12-04 Thread Pascal Oettli
Hello, Please read carefully the message. It says you three times that the name is rJava, not rjava. Regards, Pascal On 5 December 2013 13:19, tlw1987 [via R] ml-node+s789695n4681662...@n4.nabble.com wrote: Hello, everybody , recently , I want to install the package rjava in windows 7 64bt.

[R] Question about ifelse() XXXX

2013-12-02 Thread Dan Abner
Hi all, Can anyone explain what is happening with element 4,4 of c1? ifelse() is not recongizing it as value 1: c1 q1q2q3q4 q1 1.000 0.6668711 0.6948419 0.5758860 q2 0.6668711 1.000 0.6040746 0.4917447 q3 0.6948419 0.6040746 1.000 0.4730732 q4

Re: [R] Question about ifelse() XXXX

2013-12-02 Thread Duncan Murdoch
On 02/12/2013 2:08 PM, Dan Abner wrote: Hi all, Can anyone explain what is happening with element 4,4 of c1? ifelse() is not recongizing it as value 1: FAQ 7.31. Duncan Murdoch c1 q1q2q3q4 q1 1.000 0.6668711 0.6948419 0.5758860 q2 0.6668711

Re: [R] Question about ifelse() XXXX

2013-12-02 Thread William Dunlap
:09 AM To: r-help@r-project.org Subject: [R] Question about ifelse() Hi all, Can anyone explain what is happening with element 4,4 of c1? ifelse() is not recongizing it as value 1: c1 q1q2q3q4 q1 1.000 0.6668711 0.6948419 0.5758860 q2

[R] Question about compatibility

2013-11-23 Thread Juan Manuel Reyes S
Dear R- project I am beginning to work in R. When I was trying to read data for external files with command read.table in R, R was reporting me: example=read.table(file=example.text, header=TRUE,sep=,) Error in file(file, rt) : It is not posible to open connection Además: Mensajes de aviso

Re: [R] Question about compatibility

2013-11-23 Thread Patrick Burns
If you are anything like me, then the main thing that you could do to win the game is to type the name of the file correctly. Given that is a near impossibility for me, when I'm on Windows I use 'file.choose' to get the correct name. Your command would look like:

[R] Question about error of non-numeric argument to binary operator

2013-11-22 Thread Noor Aziani Bt Harun
Hi, I'm going to run my thesis that use R language for Lee-Carter method but some error occur when I'm running this code. Please guide me Sir. male.lca-lca(malaysia.male) Error in pop * mx : non-numeric argument to binary operator Regards [[alternative HTML version deleted]]

Re: [R] Question about error of non-numeric argument to binary operator

2013-11-22 Thread David Winsemius
On Nov 22, 2013, at 1:28 AM, Noor Aziani Bt Harun wrote: Hi, I'm going to run my thesis that use R language for Lee-Carter method but some error occur when I'm running this code. Please guide me Sir. male.lca-lca(malaysia.male) Error in pop * mx : non-numeric argument to binary

Re: [R] Question on xyplot

2013-11-21 Thread Duncan Mackay
-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Carl Witthoft Sent: Thursday, 21 November 2013 23:22 To: r-help@r-project.org Subject: Re: [R] Question on xyplot you didn't show us the code you used to generate the legend. I'm guessing you want to add to the legend list

Re: [R] Question on xyplot

2013-11-21 Thread Carl Witthoft
you didn't show us the code you used to generate the legend. I'm guessing you want to add to the legend list something like lty=0:7 . KB wrote I recently started using R, so I'm not really experienced with it. My question is on adjusting xyplots to get lty lines instead of coloured lines.

[R] question regarding cv.glmnet in glmnet package

2013-10-23 Thread Jack Luo
Hi, I am using glmnet for my data and have questions regarding cv.glmnet: 1. Is the 10 fold CV stratified cross validation for binary classification problem? 2. I am doing binary classification (family = binomial), the plot from cv.glmnet gives the average auc as well as the error bar with

[R] Question ANACOR-HELP

2013-10-01 Thread Matilde de Araújo e Sá Valente Rosa
Dear all, Sorry for bothering you but I am having a doubt for a long time and I cannot find a solution for it. I don't know if you can help me but I am doing a research study on Correspondence Analysis and I have been using the package ANACOR. After performing the Correspondence Analysis, I

[R] question regarding cast function in reshape package with mixture of numeric and character

2013-09-26 Thread Jack Luo
Hi, I am trying to process some data frame with cast function to unwind the stacked variables. I have no problem using cast when the values are all numeric based on the following format: cast(df,subject~v1+v2+v3,value = value,fun.aggregate = mean, fill = NA) However, I am getting some

<    1   2   3   4   5   6   7   8   9   10   >