Re: [R] sub-setting rows based on dates in R

2017-01-31 Thread Jim Lemon
Hi Md, This kind of clunky, but it might do what you want. df1<-read.table(text="DateRainfall_Duration 6/14/2016 10 6/15/2016 20 6/17/2016 10 8/16/2016 30 8/19/2016 40", header=TRUE,stringsAsFactors=FALSE) df1$Date<-as.Date(df1$Date,"%m/%d/%Y")

Re: [R] HELP with GLM

2017-01-31 Thread Bert Gunter
Dear Xavier: It sounds like you have a right mess! Perhaps others cleverer and more diligent than I can sort through it and diagnose the problem. However, it really *does* sound like you need to step back, take a deep breath, and spend some time with an R tutorial or two (there are many good ones

Re: [R] HELP with GLM

2017-01-31 Thread David Winsemius
> On Jan 31, 2017, at 6:35 PM, CHIRIBOGA Xavier > wrote: > > Dear colleagues, > > > I am trying to perform GLM ..but I got some objects masked: and an error > message below > > > a <- read.table(file.choose(), h<-T) >> head(a) > time treatment transinduc > 1

Re: [R] Failure to understand namespaces in XML::getNodeSet

2017-01-31 Thread Mark Sharp
Hadley, It’s sometimes amazing the mistakes I can make. No, it did not do what I wanted, which was read_xml(str_c(with_ns_xml, collapse = “") Reproducible example follows: library(stringr) library(xml2) ## Given the correct argument value for collapse, the next two lines work no_ns <-

[R] Error in mcp2matrix(model, linfct = linfct)

2017-01-31 Thread Karen Castillioni
Problem: Error in mcp2matrix(model, linfct = linfct) Post hoc test is not working. Does anyone know what I did wrong? modmisto<-lme(Cobertura~Tratamento, random=~1|Parcela, data=Cover_BraquiT3) summary(modmisto) tukey<-glht(modmisto, mcp(Tratamento="Tukey")) Error in mcp2matrix(model, linfct

Re: [R] caculate correlation

2017-01-31 Thread Jim Lemon
Hi Elham, It looks to me as though you are looking for matches in the name field and then asking for all columns except the first (-1). If you only have two columns in "coding.rpkm", and "name" is the first column, you will get whatever is in the second column that has a match in the "name"

[R] need help in trying out sparklyr - spark_connect will not work

2017-01-31 Thread Taylor, Ronald C
Hello R-help list, I am a new list member. My first question: I was trying out sparklyr (in R ver 3.3.2) on my Red Hat Linux workstation, following the instructions at spark.rstudio.com as to how to download and use a local copy of Spark. The Spark download appears to work. However, when I try

[R] HELP with GLM

2017-01-31 Thread CHIRIBOGA Xavier
Dear colleagues, I am trying to perform GLM ..but I got some objects masked: and an error message below a <- read.table(file.choose(), h<-T) > head(a) time treatment transinduc 11 CHA0+Db 1,0768 21 CHA0+Db 1,0706 31 CHA0+Db 1,0752 41 CHA0+Db 1,0689 5

Re: [R] Failure to understand namespaces in XML::getNodeSet

2017-01-31 Thread Hadley Wickham
I think you want x <- read_xml(' http://labkey.org/etl/xml;> MFIA 9-Plex (CharlesRiver) ') The collapse argument do what you think it does. Hadley On Tue, Jan 31, 2017 at 5:36 PM, Mark Sharp wrote: > Hadley, > > Thank you. I am able to get the xml_ns_strip() function

Re: [R] Failure to understand namespaces in XML::getNodeSet

2017-01-31 Thread Mark Sharp
Hadley, Thank you. I am able to get the xml_ns_strip() function to work with my file directly so I will likely be able to reach my immediate goal. However, I still have had no success with understanding the namespace problem. I am not able to use read_xml() using the object I generated for the

Re: [R] Challenge extracting months

2017-01-31 Thread Jim Lemon
Hi Kwesi, I worked through your code below, and I think that when you have the two variables "mon.t1" and "seas.t1" you can select a "rolling quarter" like this: # the file name in your example is different from the one you sent era<-read.table(file="SAfr_700hpa_7x5II.txt",header=FALSE,sep=" ",

Re: [R] [FORGED] filter correlation data

2017-01-31 Thread Rolf Turner
On 01/02/17 10:34, Elham - via R-help wrote: hello everybody,I have a very very huge table in R from calculating correlation,how can I filter it per spearman correlation and p-value before export it,I mean what is the function that I use?I want to select the pairs for value (r), , greater than

Re: [R] filter correlation data

2017-01-31 Thread Elham - via R-help
actually,First I searched in net,after that I sent my question, because I couldn't find the function. On Wednesday, February 1, 2017 1:34 AM, Bert Gunter wrote: ... And why did you not do a web search on "correlation coefficient in R", which would have led you

Re: [R] Failure to understand namespaces in XML::getNodeSet

2017-01-31 Thread Hadley Wickham
See the last example in ?xml2::xml_find_all or use xml2::xml2::xml_ns_strip() Hadley On Tue, Jan 31, 2017 at 9:43 AM, Mark Sharp wrote: > I am trying to read a series of XML files that use a namespace and I have > failed, thus far, to discover the proper syntax. I have a

Re: [R] filter correlation data

2017-01-31 Thread Bert Gunter
... And why did you not do a web search on "correlation coefficient in R", which would have led you almost imediately to ?cor and friends? -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in

[R] filter correlation data

2017-01-31 Thread Elham - via R-help
hello everybody,I have a very very huge table in R from calculating correlation,how can I filter it per spearman correlation and p-value before export it,I mean what is the function that I use?I want to select the pairs for value (r), , greater than 0.9 (directly correlated) and less than -0.9

[R] sub-setting rows based on dates in R

2017-01-31 Thread Md Sami Bin Shokrana
Hello guys, I am trying to solve a problem in R. I have 2 data frames which look like this: df1 <- DateRainfall_Duration 6/14/2016 10 6/15/2016 20 6/17/2016 10 8/16/2016 30 8/19/2016 40 df2 <- DateRemoval.Rate 6/17/201664.7 6/30/201622.63

Re: [R] caculate correlation

2017-01-31 Thread Elham -
Hi Dear Jim, I did it, both return a vector of name of the genes with different length,as I said before I have list of coding and noncoding so the length are not same. where is number?! and at the end of print there is this error : <0 rows> (or 0-length row.names) On Tuesday, January

[R] Failure to understand namespaces in XML::getNodeSet

2017-01-31 Thread Mark Sharp
I am trying to read a series of XML files that use a namespace and I have failed, thus far, to discover the proper syntax. I have a reproducible example below. I have two XML character strings defined: one without a namespace and one with. I show that I can successfully extract the node using

Re: [R] (no subject)

2017-01-31 Thread PIKAL Petr
Hi Of course you have fewer values than in original data. > x<-1:21 > length(rowSums(embed(x,3))) [1] 19 > length(x) [1] 21 > embed(x,3) [,1] [,2] [,3] [1,]321 [2,]432 [3,]543 [4,]654 [5,]765 [6,]876 [7,]98

Re: [R] g parameter for deltaMethod() as a function

2017-01-31 Thread Marc Girondot via R-help
Dear John and list members, I have found a solution using the package nlWaldTest. I post the solution in case someone else will have this problem. Here is a summary of the problem: I would like use the delta method for a function for which no derivative using D() can be calculated. I would

Re: [R] caculate correlation

2017-01-31 Thread Jim Lemon
Hi Elham, On Tue, Jan 31, 2017 at 7:28 PM, Elham - wrote: > Hi Dear Jim, > > I did it, both return a vector of name of the genes with different length,as > I said before I have list of coding and noncoding so the length are not > same. > > where is number?! > Not in the