Re: [R] plot in generalized additive model (GAM)

2014-07-01 Thread David Winsemius
On Jul 1, 2014, at 11:33 PM, agostinodicia...@tiscali.it wrote: > Dear David, many thanks for your reply. > The aim was to evaluate the relationship between daily mortality (variable > "mortality") and mean daily values of PM10 levels (variable "PM10") in a > specific geographic area, consideri

Re: [R] Stringr / Regular Expressions advice

2014-07-01 Thread arun
#or res <- mapply(`%in%`, accel_data, v.to.match) res1 <- sapply(seq_len(ncol(accel_data)),function(i) accel_data[i]<=tail(v.to.match[[i]],1) & accel_data[i] >=v.to.match[[i]][1]) all.equal(res, res1,check.attributes=F) #[1] TRUE A.K. On Tuesday, July 1, 2014 10:56 PM, arun wrote: Hi Vincen

Re: [R] x axis labelling

2014-07-01 Thread Duncan Mackay
Hi Michael Dates and times are always a problem as they are irregular not 1,2,3 ..., 100 If you want more fancy formatting of the x axis try this First convert your time to a datetime class # Use a dummy date for datetime as it is easier mmdat$time <- seq(strptime("20140702 14", "%Y%m%d %H"), by

Re: [R] plot in generalized additive model (GAM)

2014-07-01 Thread David Winsemius
On Jul 1, 2014, at 12:02 AM, adc wrote: > I performed the following GAM by the MGCV package: I think it's actually spelled in all lower case. > gam(mortality ~ > (PM10) + (Tmax) + (umidity), data = data, family = quasipoisson) > > How > can I obtain a plot of Log-relative risk of mortality

Re: [R] Data visualization: overlay columns of train/test/validation datasets

2014-07-01 Thread David Winsemius
On Jul 1, 2014, at 3:46 PM, Supriya Jain wrote: > Hello, > > Given two different datasets (having the same number and type of columns, > but different observations, as commonly encountered in data-mining as > train/test/validation datasets), is it possible to overlay plots > (histograms) and com

[R] Data visualization: overlay columns of train/test/validation datasets

2014-07-01 Thread Supriya Jain
Hello, Given two different datasets (having the same number and type of columns, but different observations, as commonly encountered in data-mining as train/test/validation datasets), is it possible to overlay plots (histograms) and compare the different attributes from the separate datasets, in o

[R] Using RCMD INSTALL under Spanish version of windows.

2014-07-01 Thread Jan Graffelman
Using RCMD INSTALL with R-version 3.1.0 under a Spanish Windows 7 gives the following error message: rcmd INSTALL MyPackages Mensajes de aviso perdidos In normalizePath(path.expand(path), winslash, mustWork) : path[1]="c:/ARCHIV~1/R/R-31~1.0/library": Acceso denegado Mensajes de aviso perdidos

Re: [R] x axis labelling

2014-07-01 Thread Jim Lemon
On Tue, 1 Jul 2014 06:41:52 PM Michael Millar wrote: > Hi, > > I am new to R and am trying to create a graph with Time(24hr) along the x > axis. Rather than start at 01.00, I wanted to start at 14.00. > > I tried to use the axis(side=1, at=c( )) function but it continues to put > then in numeri

Re: [R] an incredibly trivial question about nls

2014-07-01 Thread Erin Hodgess
Thank you to all. I had actually found the summary and trotted that out. Just had not gotten back to the list. Thanks again! Sincerely, Erin On Tue, Jul 1, 2014 at 5:46 PM, Bert Gunter wrote: > Beauty -- or obscurity -- is in the eyes of the beholder. But I leave > your objections to stand

Re: [R] an incredibly trivial question about nls

2014-07-01 Thread Rolf Turner
In direct contrast to what Bert says, I think this is a very reasonable (and non-trivial) question. The problem results from Gurus structuring the functions that they write in such a way that they are totally opaque to anyone but the ultra-cognoscenti. What is gained by not having things s

Re: [R] an incredibly trivial question about nls

2014-07-01 Thread Peter Langfelder
On Tue, Jul 1, 2014 at 1:27 PM, Erin Hodgess wrote: > Hello R People: > > I'm having a forest/trees location problem with the output of nls. > > If I save the output to an object, and print the object, it shows, amongst > other things, the residual sum of squares. I would like to get that. > > Ho

Re: [R] an incredibly trivial question about nls

2014-07-01 Thread Bert Gunter
1. Why? What do you think it tells you? (The number of parameters in a NONlinear model is probably not what you think it is). 2. ?deviance 3. You've been posting all this time and still didn't try stats:::print.nls ?? -- which is where you would find the answer. Cheers, Bert Bert Gunter Gene

Re: [R] Fwd: combining data from multiple read.delim() invocations.

2014-07-01 Thread Bert Gunter
On Tue, Jul 1, 2014 at 12:03 PM, John McKown wrote: > On Tue, Jul 1, 2014 at 11:31 AM, David L Carlson wrote: > >> There is a better way. First we need some data. This creates three files >> in your home directory, each with five rows: >> >> write.table(data.frame(rep("A", 5), Sys.time(), Sys.tim

[R] an incredibly trivial question about nls

2014-07-01 Thread Erin Hodgess
Hello R People: I'm having a forest/trees location problem with the output of nls. If I save the output to an object, and print the object, it shows, amongst other things, the residual sum of squares. I would like to get that. However, when I look at names or str of the object, I can't find the

[R] Fwd: combining data from multiple read.delim() invocations.

2014-07-01 Thread John McKown
On Tue, Jul 1, 2014 at 11:31 AM, David L Carlson wrote: > There is a better way. First we need some data. This creates three files > in your home directory, each with five rows: > > write.table(data.frame(rep("A", 5), Sys.time(), Sys.time()), > "A.tab", sep="\t", row.names=FALSE, col.name

Re: [R] Generating Patient Data

2014-07-01 Thread David Winsemius
On Jun 25, 2014, at 1:49 PM, David Winsemius wrote: > > On Jun 24, 2014, at 11:18 PM, Abhinaba Roy wrote: > >> Hi David, >> >> I was thinking something like this: >> >> ID Disease >> 1 A >> 2 B >> 3 A >> 1C >> 2D >> 5A >> 4B >> 3D >> 2A >> .... >> >>

Re: [R] x axis labelling

2014-07-01 Thread Michael Peng
Did you add xaxt = "n" in the plot function? Try the following: plot(x,y, xaxt = "n") axis(1, at = c(14, 20),labels = c("14h", "20h") ) 2014-07-01 12:41 GMT-05:00 Michael Millar : > Hi, > > I am new to R and am trying to create a graph with Time(24hr) along the x > axis. Rather than start at 01

[R] Using CSS package to extract text from html

2014-07-01 Thread dstrick1
This being my first post, I'm sure I'll do something discordant with convention, so forgive me in advance. Basically, I am trying to extract text from an html file using the CSS package in R. However, I am unable to do so because it seems that the text itself is not identified with any class and t

[R] R help

2014-07-01 Thread Andre Weeks
To whom it may concern: I installed R 3.1 and I get this. In normalizePath(path.expand(path), winslash, mustWork) : path[1]="\\network\users\aweeks\My Documents/R/win-library/3.1": Access is denied Is there any way to change this path? I have looked it up on the internet but cannot seem to fin

Re: [R] 1-dinemsional point process

2014-07-01 Thread David Winsemius
It's unclear why density estimates are not being mentioned. Also suggest you search: install.packages("sos") require(sos) findFn("scan statistic") On Jun 30, 2014, at 7:35 PM, Doobs wrote: > Hi, > As a new user, is it possible to look at clustering/dispersion processes of > a 1D point process

[R] x axis labelling

2014-07-01 Thread Michael Millar
Hi, I am new to R and am trying to create a graph with Time(24hr) along the x axis. Rather than start at 01.00, I wanted to start at 14.00. I tried to use the axis(side=1, at=c( )) function but it continues to put then in numeric order. Is there another way I can add labels to the x axis? Tha

Re: [R] combining data from multiple read.delim() invocations.

2014-07-01 Thread David L Carlson
I agree it is not necessarily faster, but the code is more compact since we don't have to initialize the variable or explicitly refer to the index. For big data it has the disadvantage of storing the data twice. For speed, this is faster and does not store the data twice, but is system depende

Re: [R] Stringr / Regular Expressions advice

2014-07-01 Thread VINCENT DEAN BOYCE
Sara, Yes, I modified the code that you provided and it worked quite well. Here is the revised code: . accel_data <- data *# pattern to be identified* v.to.match <- c(438, 454, 459) # call the below function anytime the "v.to.match" criteria changes to ensure match is updated v.matches <- ap

Re: [R] combining data from multiple read.delim() invocations.

2014-07-01 Thread Bert Gunter
Maybe, David, but this isn't really it. Your code just basically reproduces the explicit for() loop with the lapply. Maybe there might be some advantage in rbinding the list over incrementally adding rows to the data frame, but I would be surprised if it made much of a difference either way. Of c

Re: [R] combining data from multiple read.delim() invocations.

2014-07-01 Thread David L Carlson
There is a better way. First we need some data. This creates three files in your home directory, each with five rows: write.table(data.frame(rep("A", 5), Sys.time(), Sys.time()), "A.tab", sep="\t", row.names=FALSE, col.names=FALSE) write.table(data.frame(rep("B", 5), Sys.time(), Sys.time(

Re: [R] logistic regression for data with repeated measures

2014-07-01 Thread Mitchell Maltenfort
http://stats.stackexchange.com/questions/62225/conditional-logistic-regression-vs-glmm-in-r might be a good start Ersatzistician and Chutzpahthologist I can answer any question. "I don't know" is an answer. "I don't know yet" is a better answer. On Tue, Jul 1, 2014 a

[R] logistic regression for data with repeated measures

2014-07-01 Thread Suzon Sepp
Hi, It seems that I'm quite lost in this wide and powerful R's universe, so I permit myself to ask your help about issues with which I'm struggling. Thank you, I would like to know if the answer’s accuracy (correct = 1; incorrect = 0) varies depending on 2 categorical variables which are the gr

[R] combining data from multiple read.delim() invocations.

2014-07-01 Thread John McKown
Is there a better way to do the following? I have data in a number of tab delimited files. I am using read.delim() to read them, in a loop. I am invoking my code on Linux Fedora 20, from the BASH command line, using Rscript. The code I'm using looks like: arguments <- commandArgs(trailingOnly=TRUE

[R] Socket Connection in R

2014-07-01 Thread Param Jeet
I am trying to create socket connection in R. socket <- make.socket("localhost",2099,T,T) msg2<-'function=subscribe|item=MI.EQCON.1|schema=last_price;ask;bid' write.socket(socket,msg2) read.socket(socket,252,FALSE) When I run the read.socket line, I get error: Error in read.s

[R] Order Book details in R Interactive Brokers Package

2014-07-01 Thread jeeth ghambole
Hello All, I am working on project of Automated trade execution using R and Interactive Brokers Package. I have successfully implemented and tested the connection of R with Interactive Brokers API. Implementing orders, placing orders too are working fine. The only problem is that while executing

[R] plot in generalized additive model (GAM)

2014-07-01 Thread adc
I performed the following GAM by the MGCV package: gam(mortality ~ (PM10) + (Tmax) + (umidity), data = data, family = quasipoisson) How can I obtain a plot of Log-relative risk of mortality vs. PM10 ? thanks agostino Scopri istella, il nuovo motore per il web italiano. Istella garantisc

[R] (PLM- package) Residual-Plotting and missing Values

2014-07-01 Thread Katharina Mersmann
Dear R-Community, I tried plotting the residuals of an FE-model estimated via plm . And detected that there are no residuals in the plot for the last two countries. I guess this happens because for some countries values are missing and R gives me the following for > fixed.reg1.new$resid[1]

Re: [R] c() with POSIXlt objects and their timezone is lost

2014-07-01 Thread Uwe Ligges
On 23.06.2014 23:52, Marc Girondot wrote: When two POSIXlt objects are combine with c(), they lost their tzone attribute, even if they are the same. I don't know if it is a feature, but I don't like it ! Marc > es <- strptime("2010-02-03 10:20:30", format="%Y-%m-%d %H:%M:%S", tz="UTC") > es

Re: [R] Dead link in the help page of as.Date()

2014-07-01 Thread Uwe Ligges
On 23.06.2014 21:41, Christofer Bogaso wrote: Hi, I was reading the help page for as.Date() function for some reason, and noticed a Matlab link: http://www.mathworks.com/help/techdoc/matlab_prog/bspgcx2-1.html Thanks, updated now. Best, Uwe Ligges It looks like this link is dead. So ma

Re: [R] How to plot individual pdf files for each wrapped plot with ggplot2?

2014-07-01 Thread Bea GD
Just solved the first problem! I had to generate a plot and then plotted. Now it's saved into pdf. Only the second issue: *2. I don't know how to write the code to get the 8 sex/day combinations.* Thanks! On 01/07/2014 12:59, Bea GD wrote: > Thanks a lot for your reply Trevor! > > I've been wo

Re: [R] How to plot individual pdf files for each wrapped plot with ggplot2?

2014-07-01 Thread Bea GD
Thanks a lot for your reply Trevor! I've been working with the code but I cannot make it work. I have 2 main problems: 1. From running the loop I get pdf files with no pages generated. 2. I don't know how to write the code to get the 8 sex/day combinations. library(ggplot2) library(resha

[R] Using external SQLite installation for RSQLite in Windows?

2014-07-01 Thread Nick Eubank
Hello All, I'm trying to figure out how to link RSQLite to an external sqlite3 installation compiled for a 64bit platform. I see from the CRAN installation instructions (see: http://cran.r-project.org/web/packages/RSQLite/INSTALL) that on a unix machine there's a way to set the configuration to ac