Re: [R] Dataverse (reading files with .tab and .7z suffixes)

2018-05-13 Thread Thomas Levine
Ilio Fornasero writes: > Yet, I am at this point. > > > > > ## 01. Finding the dataverse server and making a search > Sys.setenv("DATAVERSE_SERVER" =3D "dataverse.harvard.edu") > dataverse_search(".Hunger") > > > ## 02. Loading the dataset (in this example, I have chosen the word ".Hunge= > r" to g

[R] Dataverse (reading files with .tab and .7z suffixes)

2018-05-13 Thread Thomas Levine
Ilio Fornasero writes: > I am trying to find a way to retrieve data from Harvard Dataverse website. > I usually don't have problem in web-scraping data but the problem here is > that there are a bunch of data formats such as .tab, .7z and so and > I just can't find a way to retrieve the data I am

Re: [R] R-help Digest, Vol 183, Issue 13

2018-05-13 Thread Thomas Levine
Ilio Fornasero writes: > Hello. > > I am trying to find a way to retrieve data from Harvard Dataverse website. > I usually don't have problem in web-scraping data but the problem here is > that there are a bunch of data formats such as .tab, .7z and so and I just > can't find a way to retrieve

Re: [R] [FORGED] Rotating points, preferably in grid

2016-08-24 Thread Thomas Levine
Well this is great. Now I have answers for both graphics and grid. The rot argument is exactly what I had wanted, except that I had imagined it also working on points. But I had not thought to use unicode, and that will probably make this plot even easier. Thanks

[R] Rotating points, preferably in grid

2016-08-24 Thread Thomas Levine
I want to make a plot in polar coordinates. I want to use pch with shapes that do not have radial symmetry, so I want to rotate them such that they face inwards. I am using grid for my plotting, but I provide motivating examples in graphics. The following plot almost gets me what I want. theta

Re: [R] Making my own graphics device

2014-08-17 Thread Thomas Levine
Thanks for the further tips! I'll look at the R Internals manual, the other built-in devices, and the RSvgDevice package. And I'm still looking for the S source code, but I have a feeling that I'll wind up writing the device before I find the source code. > It may be useful to look at existing int

Re: [R] Making my own graphics device

2014-08-17 Thread Thomas Levine
of Economics > vox: 217-333-4558University of Illinois > fax: 217-244-6678Urbana, IL 61801 > On Aug 17, 2014, at 2:58 PM, Thomas Levine <_...@thomaslevine.com> wrote: > > > I want to make my own graphics device am thus

[R] Making my own graphics device

2014-08-17 Thread Thomas Levine
I want to make my own graphics device am thus looking for documentation about graphics devices. The only thing I've found so far is these directions for making graphics devices with the RGraphicsDevice package. http://www.omegahat.org/RGraphicsDevice/ Could someone point me to any other resources

Re: [R] HPGL or PCL plotting device? Or otherwise plotting plots

2014-07-09 Thread Thomas Levine
ostscript files for the different colors on the plot, create the corresponding HPGL files, edit the SP command in each of them, and concatenate them. But maybe there's a better way? On 09 Jul 13:32, Thomas Levine wrote: > Oh it was easier than I thought. > > postscript(

Re: [R] HPGL or PCL plotting device? Or otherwise plotting plots

2014-07-09 Thread Thomas Levine
Oh it was easier than I thought. postscript('project-contracts.ps') hist(log(projects$n.contracts)) dev.off() Then run this from the shell. pstoedit -f plot-hpgl project-contracts.ps project-contracts.hpgl And send it to the plotter. On 09 Jul 13:10, Thomas Levine wrote: &g

[R] HPGL or PCL plotting device? Or otherwise plotting plots

2014-07-09 Thread Thomas Levine
Hi, I want to print plots on a Roland DXY-1100 plotter. How can I do this from R? I think the easiest thing would be a graphics device for Printer Command Language or Hewlett-Packard Graphics Language, but I haven't managed to find any of those. Thanks Tom __

[R] Stripchart colors don't vary after I sort a data frame

2012-11-28 Thread Thomas Levine
# Hi, # This plot has two colors. overflow <- read.csv('http://chainsaw.thomaslevine.com/overflow.csv', stringsAsFactors = F) png('original.png') stripchart(overflow$precipi ~ overflow$after.9.am, method='stack', pch = 22, bg = overflow$overflow + 1, vertical = T, col = 0) dev.off() # I wan

[R] Can I specify POSIX[cl]t column classes inside read.csv?

2012-04-23 Thread Thomas Levine
I'm loading a nicely formatted csv file.     #!/usr/bin/env Rscript     kpi <- read.csv(       # This is a dump of the username, date_joined and last_login columns       # from the auth_user Django table.       'data/2012-04-23.csv',       colClasses = c('character')     )     print(kpi[sample(nro

Re: [R] R development master class: NYC, Dec 12-13

2011-11-15 Thread Thomas Levine
Hmm... you know I'd love to run a study on how software and other information displays affect the speed, accuracy and reliability with which people make insights about data. Tom On Tue, Nov 15, 2011 at 4:44 PM, Stavros Macrakis wrote: > > Last time, I was told that I couldn't list my R package a

Re: [R] squared "pie chart" - is there such a thing?

2011-07-23 Thread Thomas Levine
How about just a stacked bar plot? barplot(matrix(c(3,5,3),3,1),horiz=T,beside=F) Tom On Fri, Jul 22, 2011 at 7:14 AM, Naomi Robbins wrote: > > Hello! > It's a shoot in the dark, but I'll try. If one has a total of 100 > (e.g., %), and three components of the total, e.g., > mytotal=data.frame(x

Re: [R] computer name

2011-06-12 Thread Thomas Levine
Not exactly R, but how about > pcname <- system('uname -n',intern=T) Tom On Sun, Jun 12, 2011 at 11:19 PM, pdb wrote: > > Is there an r function that will be able to identify the computer the code is > running on? > > I have some common code that I run on several computers and each has a > data

Re: [R] Odp: Three sigma rule

2011-05-31 Thread Thomas Levine
I think you really want a normality test. If that's what you want, you have more options than the three-sigma rule. http://en.wikipedia.org/wiki/Normality_test Tom On Tue, May 31, 2011 at 12:31 PM, Bert Gunter wrote: > Folks: > > On Tue, May 31, 2011 at 8:48 AM, Petr PIKAL wrote: >> Hi >> >> r

Re: [R] help with barplot

2011-05-28 Thread Thomas Levine
You can do pretty well without ggplot actually. boxplot(Time~paste(Incidents,Months),data=DF,border=c('grey20','red')) On Sat, May 28, 2011 at 2:55 AM, steven mosher wrote: > Thanks, > >  ggplot is on my list of things to learn before Hadley comes here to the > bay area >  to give a session on i

Re: [R] help with barplot

2011-05-27 Thread Thomas Levine
Does this do it? barplot(t(matrix(pivot$x,4)),beside=T) Tom On Fri, May 27, 2011 at 6:08 PM, steven mosher wrote: > > Hi, > > I'm really struggling with barplot > > I have a data.frame with 3 columns. The first column represents an > "incident" type > The second column represents a "month" > Th

Re: [R] Null

2011-05-05 Thread Thomas Levine
Maybe you were doing something like fcv <- read.csv('fcv.csv') instead of fcv <- read.csv('fcv.csv')[1] (I haven't tested this.) Tom On Thu, May 5, 2011 at 8:48 AM, pcc wrote: > > This is probably a very simple question but I am completely stumped!I am > trying to do shapiro.wilk(x) test on

Re: [R] Generalized Linear Model

2011-04-27 Thread Thomas Levine
Because you have two dependent variables, you'll want to to use a multivariate logit. mlogit does this, but I don't know the syntax off hand. If you just wanted to look at one dependent variable, it would be the following (which Alex said) glm(y~x1*x2,family='binomial') On Mon, Apr 25, 2011 at 3

Re: [R] How to erase (replace) certain elements in the data.frame?

2011-04-24 Thread Thomas Levine
As Joshua said, mine was indeed different from yours. And it didn't work on non-numeric data. But this one seems to work right: random.del_vec <- function (x, n.keeprows, del.percent){ del<-function(notkeep){ k<-floor(length(notkeep)*del.percent/100) notkeep[sample.int(length(notkeep),k)

Re: [R] How to erase (replace) certain elements in the data.frame?

2011-04-24 Thread Thomas Levine
This should do the same thing random.del <- function (x, n.keeprows, del.percent){   del<-function(col){     col[sample.int(length(col),length(col)*del.percent/100)]<-NA     col   }   change<-n.keeprows:nrow(x)   x[change,]<-lapply(x[change,],del)   x } This is faster because it's vectorized. [1

Re: [R] Sample size estimation for sample surveys

2011-04-04 Thread Thomas Levine
hits, many of which could be relevant. > > HTH, > Dennis > > On Mon, Apr 4, 2011 at 6:05 PM, Thomas Levine > wrote: >> >> Hi, >> >> Is there an R package for estimating sample size requirements for >> parameter estimation in sample surveys? In particula

[R] Sample size estimation for sample surveys

2011-04-04 Thread Thomas Levine
Hi, Is there an R package for estimating sample size requirements for parameter estimation in sample surveys? In particular, I'm interested in sample size estimation for stratified and systematic sampling. I have a textbook with appropriate formulae, but it'd be nice if I didn't have to type in al

Re: [R] simple if question

2011-03-26 Thread Thomas Levine
Posting some sample data would help, but I think something like this is what you want data[data$school=='Cornell University',] For example CO2[CO2$Type=='Quebec',] Tom 2011/3/26 Sebastián Daza : > Hi everyone, > I have just got different samples from a dataframe (independent and > exclusive, t

Re: [R] Two matrix loop

2011-03-26 Thread Thomas Levine
This way uses a three-dimensional array instead of the nested apply. It seems to take the same amount of time, even on larger datasets, but it may give you ideas. distance=function(x) daisy(x, metric = 'gower') persons=array(dim=c(2,nrow(donor)*nrow(receiver),ncol(receiver))) persons[1,,]=donor[r

Re: [R] Extracting columns from a class

2011-03-19 Thread Thomas Levine
Hi, Here is the prcomp output. tom=prcomp(matrix(rnorm(25),5,5)) R functions often output lists. To see what's in this one, run "names(tom)" or type "tom$" and use tab completion. Once you do that, the following is more obvious. pc1=tom$rotation[,1] sd1=tom$sdev[1] column=c(sd1,pc1) write.tabl

Re: [R] Easy help with plots, font size

2011-02-03 Thread Thomas Levine
There's also pointsize, depending on your device. On Thu, Feb 3, 2011 at 9:31 PM, gigi1234 wrote: > > Ok thanks. I did have some success using the cex options. Figured out that > I > was putting some of the arguments in the wrong place. > > Thanks again! > > On Thu, Feb 3, 2011 at 8:24 PM, David

[R] Converting numbers into words

2010-12-05 Thread Thomas Levine
Example data desk=data.frame( deskchoice=c('mid','mid','left','bookdrop','mid','bookdrop') ) -- I like doing stuff like the line below, especially when I'm using Sweave. print(paste('Within the observation period,',nrow(desk), 'patrons approached the circulation desk.')) -- But what if I wan

Re: [R] R script on linux?

2010-11-02 Thread Thomas Levine
Open a terminal, then run these two commands. cd /home/the/directory/with/your/script R Then run this in R source('yourscript.R') Tom 2010/11/2 Jonathan P Daily : > What is the error message? > -- > Jonathan P. Daily > Technician - USGS Leetown Science Cente

Re: [R] make many barplot into one plot

2010-10-31 Thread Thomas Levine
hierobarp or barNest from {plotrix} may do this more neatly. 2010/10/31 Sibylle Stöckli > Dear R users > > > I would like to group my barplot graph (see example on the R help > link). The proposed R code, adding individual bars to the plot, looks > really overwhelming. My specific dataset just

Re: [R] R and Matlab

2010-10-28 Thread Thomas Levine
Is there a particular reason you can't use csv? write.csv() in R It seems that you can read csv in Matlab with this http://www.mathworks.com/help/techdoc/ref/importdata.html Tom 2010/10/28 Claudia Beleites : >> I am looking for ways to use R and Matlab. Doing the data transformations >> in >> R

Re: [R] times

2010-10-28 Thread Thomas Levine
I couldn't figure that out, but you can work around it by adding the seconds chron(dates=lwc.file[,1],times=paste(lwc.file[,2],':00',sep='')) Tom 2010/10/28 thoeb : > > Hi! I have an input table with a column "Dates" in the format > "month/day/year" (eg. "5/11/1999" and a column "Times" in the f

Re: [R] Multiple, multiple regressions...

2010-10-27 Thread Thomas Levine
Write a function() that does everything you want to do for one regression, then run that on all of them. It'll look something like this. library(car) myreg=function(X,v1,v2,v3,v4,v5) { reg=lm(X ~ v1 + v2 + v3 + v4 + v5) crPlots(reg) } Then run the function a lot, maybe in a loop or with sapply()

Re: [R] Changing origin of line in radial plot

2010-10-27 Thread Thomas Levine
I guess you have something like this. > testlen <- c(sin(seq(0,1.98*pi,length=100))+2+rnorm(100)/10) > testpos <- seq(0,1.98*pi,length=100) > radial.plot(testlen,testpos,rp.type="p",main="Test Polygon",line.col="blue") (http://addictedtor.free.fr/graphiques/graphcode.php?graph=75) Does this work

Re: [R] Y Axis Labels

2010-10-27 Thread Thomas Levine
Oops. That was backwards plot(scale~time,axes=F) axis(2,at=0:4,labels=c("Never", "Once per month", "A few times per month", "A few times per week","Everyday")) axis(1) Tom 2010/10/27 Thomas Levine : > More specifically > > time=rnorm(2

Re: [R] lineplot.CI {sciplot}: continuous line

2010-10-27 Thread Thomas Levine
0.05, pch = > c(15,15), col=c("grey", "black"), lty=c(1,1), main = "Experiments 2 and 3", > ylim=c(0.3,1), fixed=TRUE, ci.fun= function(x) c(mean(x)-1.96*se(x), > mean(x)+1.96*se(x)),leg.lab=c("cond2", "cond3")) > > Thanks, > Fotis

Re: [R] Embedding graphics in a pdf()

2010-07-14 Thread Thomas Levine
Woah! That's so awesome! And now I've found even more functions of my drawing programs that can be replaced with R. Tom 2010/7/14 Marc Schwartz > On Jul 14, 2010, at 1:38 PM, Thomas Levine wrote: > > > I've had two reasons for wanting to embed graphics in R

[R] Embedding graphics in a pdf()

2010-07-14 Thread Thomas Levine
I've had two reasons for wanting to embed graphics in R pdf output. 1. I am plotting something on top of a surface (It's actually a desk.) of which I have a picture and would like to place a picture underneath. 2. I can produce all of my presentation slides in R without LaTeX but have a few pictur

Re: [R] Passing a LaTeX length to Sweave

2010-07-10 Thread Thomas Levine
n Murdoch : > On 10/07/2010 10:35 AM, Thomas Levine wrote: >> >> I'd really love to do this >> >> \setkeys{Gin}{width=\columnwidth} >> >> <>= >> plot(y~x) >> @ >> >> The \columnwidth macro does not work, of course, in the secon

[R] Passing a LaTeX length to Sweave

2010-07-10 Thread Thomas Levine
I'd really love to do this \setkeys{Gin}{width=\columnwidth} <>= plot(y~x) @ The \columnwidth macro does not work, of course, in the second line. What can I do instead? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

Re: [R] SAS for R-users

2010-05-23 Thread Thomas Levine
Thanks for the suggestions! This will keep me busy for a while. Tom 2010/5/15 Muenchen, Robert A (Bob) : >> Thomas Levine wrote: >>Bob Muenchen says that 'Ralph O’Brien says that >>in a few years there will be so many students >>graduating knowing mainly R that [he]

Re: [R] Fitting a function

2010-05-14 Thread Thomas Levine
be a function that does this already, but I couldn't find it either. Sorry for sending so many emails 2010/5/14 Thomas Levine : > I spoke too soon; that is not linear. > > I've never used the model-fitting functions, but you may check them out. > http://developer.r-project.org/

Re: [R] Fitting a function

2010-05-14 Thread Thomas Levine
I spoke too soon; that is not linear. I've never used the model-fitting functions, but you may check them out. http://developer.r-project.org/model-fitting-functions.txt Tom 2010/5/14 Thomas Levine : > x<-c(0.5,4,6,8,12) > y<-c(0.021,0.021,0.020,0.018,0.012) > lm(y~x

Re: [R] Fitting a function

2010-05-14 Thread Thomas Levine
x<-c(0.5,4,6,8,12) y<-c(0.021,0.021,0.020,0.018,0.012) lm(y~x) 2010/5/14 Dani Valverde : > Hello, > It is a very naive question, but here it is. I have this values: > > x: 0.5 4 6 8 12 > > y: 0.021 0.021 0.020 0.018 0.012 > > I need to fit a function to them. How can I do it with R? > Thank you so

Re: [R] Fwd: nonlinearity and interaction

2010-05-14 Thread Thomas Levine
Creating the 5 indicator variables will be easy if you post your code and sample data. This may also allow people to help with the first problem you were having. Tom 2010/5/14 William Simpson : > [posted this at 9:25 and still hasn't appeared on the list at 13:26] > > > I have the following set-u

[R] SAS for R-users

2010-05-14 Thread Thomas Levine
There are loads of resources for users of any other statistics package who are learning R. For example http://www.google.com/search?q="r+for+sas-users"; The reverse isn't the case http://www.google.com/search?q="sas+for+r-users"; Having heard rumors of how unelegant other statistics packages ar

[R] Bar plots with bars made of stacked text

2010-03-31 Thread Thomas Levine
I would like to make bar plots where the bars are composed of text like this: http://www.thomaslevine.com/lowres/text_bars.png Is there a package that will help me with this? Thanks Tom __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

[R] Error in tapply when reordering levels of a factor

2010-02-27 Thread Thomas Levine
I have this > grades$grade ... [4009] A B A- A- A- B+ A A- B+ B A B B B A A- A A- A- B+ A- A A B+ [4033] A- A- A- A A- B A A A- A Levels: A A- A+ B B- B+ C C+ I want to change the order of the levels > reorder(grades$grade,c('A+','A','A-','B+','B','B-','C+','C')) Error in tapply

Re: [R] Change the scale on a barplot's y axis

2010-02-27 Thread Thomas Levine
s the axes in the > barplot using axes=FALSE, and then add the axes using axis(1) and > axis(2,..) with the ... as above. > > Incidentally, I'm not sure I'd have converted your numbers that way, but > if it's worked it's worked. > > Steve E >>>> Thomas

[R] Change the scale on a barplot's y axis

2010-02-27 Thread Thomas Levine
I have grades data. I read them from a csv in letter-grade format. I then converted them to levels levels(grades$grade)=c('A+','A','A-','B+','B','B-','C+','C','C-','D+','D','D-') And then to numbers grades$gp=grades$grade levels(grades$gp)=c(4.3,4.0,3.7, 3.3,3.0,2.7, 2.3,2.0,1.7, 1.3,1.0,0.7) gr

[R] A graphic art piece I created that only statisticians would understand and the use of R in art

2009-06-29 Thread Thomas Levine
This isn't the typical r-help post, but I think people here will appreciate it. The current Shirt.Woot (http://shirt.woot.com) derby prompt is fake band names, so I chose the name "The Probable Error of a Mean." And I made my submission with R. Here it is: http://www.thomaslevine.org/the_probable_

Re: [R] How to import timestamps from emails into R

2009-06-24 Thread Thomas Levine
One last tiny problem: How do I add months to the scale? It currently just has years http://school.thomaslevine.org/mywall.png Thanks again Tom On Sat, Jun 20, 2009 at 12:14 PM, Thomas Levine wrote: > I wasn't really thinking that far ahead; plot tries to do something, so I > figu

Re: [R] How to import timestamps from emails into R

2009-06-20 Thread Thomas Levine
ng(tt)), putting the time axis along the bottom. Thanks Tom On Sat, Jun 20, 2009 at 11:15 AM, Gabor Grothendieck < ggrothendi...@gmail.com> wrote: > Try this: > > plot(seq_along(tt), tt) > > > On Sat, Jun 20, 2009 at 10:55 AM, Thomas Levine > wrote: > > Here'

Re: [R] How to import timestamps from emails into R

2009-06-20 Thread Thomas Levine
0 > Wed, 17 Jun 2009 21:50:41 -0700" > > # L <- readLines("myfile.txt") > L <- readLines(textConnection(Lines)) > tt <- as.POSIXct(L, format = "%a, %d %b %Y %H:%M:%S") > > > > On Fri, Jun 19, 2009 at 6:06 PM, Thomas Levine > wrote: >

[R] How to import timestamps from emails into R

2009-06-19 Thread Thomas Levine
I am analysing occurrences of a phenomenon by time, and each of these timestamps taken from email headers represents one occurrence. (The last number is the time zone.) I can easily change the format. Sun, 14 Jun 2009 07:33:00 -0700 Sun, 14 Jun 2009 08:35:10 -0700 Sun, 14 Jun 2009 21:26:34 -0700 M

Re: [R] Labeling barplot bars by multiple factors

2009-05-28 Thread Thomas Levine
>> [mailto:r-help-boun...@r-project.org] On Behalf Of Thomas Levine >> Sent: Thursday, May 28, 2009 5:04 AM >> To: Jim Lemon >> Cc: r-help@r-project.org >> Subject: Re: [R] Labeling barplot bars by multiple factors >> >> Both of those worked, but hierobarp l

Re: [R] Labeling barplot bars by multiple factors

2009-05-28 Thread Thomas Levine
Both of those worked, but hierobarp looked a bit easier, so I used that. The one annoying thing is that it sorts alphabetically. Tom On Thu, May 28, 2009 at 6:46 AM, Jim Lemon wrote: > Thomas Levine wrote: > >> I want to plot quantitative data as a function of three two-level fac

[R] Labeling barplot bars by multiple factors

2009-05-27 Thread Thomas Levine
I want to plot quantitative data as a function of three two-level factors. How do I group the bars on a barplot by level through labeling and spacing? Here 's what I'm thinking of. Also, I'm pretty sure that I want a barplot, but there

Re: [R] Error in FUN with tapply and by

2009-05-22 Thread Thomas Levine
onPrefHand": 1 1 1 1 1 1 1 1 1 1 ... > str(Grip) Factor w/ 2 levels "Round","Tripod": 1 1 1 1 1 1 1 1 1 1 ... On Fri, May 22, 2009 at 2:46 PM, jim holtman wrote: > You need to supply str for the original arguments; the error message had a > different set of parameter

[R] Error in FUN with tapply and by

2009-05-22 Thread Thomas Levine
ne.integrated,mean,na.rm=T)[2]-tapply(time,Baseline.integrated,mean,na.rm=T)[1])} To help explain what the truecost function does: > tapply(t_p1,Baseline.integrated,mean,na.rm=T) Baseline Integrated 212.8000 252.8402 > truecost(t_p1) [1] 40.04021 Then I

Re: [R] Plotting multiple ablines

2009-04-03 Thread Thomas Levine
Thu, Apr 2, 2009 at 8:29 AM, r...@quantide.com wrote: > May be: > > plot(c(-1, 1) , c(-1, 1), type = "n") > n = 4 > a = rep(0, n) > b = 1:n/n > > > fun = function(i, a, b, col = 1 , ...) { > abline(a[i], b[i], col = col[i], ...) > } > > lapply(1:

[R] Plotting multiple ablines

2009-04-01 Thread Thomas Levine
I really want to do this: abline( a=tan(-kT*pi/180), b=kY-tan(-kT*pi/180)*kX ) where kX,kY and kT are vectors of equal length. But I can't do that with abline unless I use a loop, and I haven't figured out the least unelegant way of writing the loop yet. So is there a way to do this without a loo

[R] Plot the highest point in a contour plot

2009-03-26 Thread Thomas Levine
I'm plotting contour plots with contourplot. which.max gives me the index of the highest point of a matrix. I can find the point in the matrix from here, convert it to a point on a graph and add it with points, but you'd think someone's already written this. How do I plot the highest point in a c

Re: [R] Contour plots of four two-dimensional matrices

2009-03-15 Thread Thomas Levine
tion will be different for > those two options. > > contourplot is is a lattice plotting function. See Figure 6.10 on Sarkar's > Lattice book pages. levelplot is the closest analog to filled contour in > lattice. > -- > David Winsemius > > > > On Mar 15, 2009,

[R] Contour plots of four two-dimensional matrices

2009-03-15 Thread Thomas Levine
I have four large two-dimensional matrices of which I want to create contour plots. Something like filled.contour() contourplot() works but only gives me one plot at a time. If I combine the four matrices into one three-dimensional matrix, which I'll name "seven", there should be a way of doing s

Re: [R] Inefficiency of SAS Programming

2009-03-02 Thread Thomas Levine
R depends on all of those things to run, but you only have to use those programs through R. The software depends on these other tools, but the human doesn't have to switch interfaces. Tom! On Fri, Feb 27, 2009 at 9:22 PM, Gabor Grothendieck wrote: > On Fri, Feb 27, 2009 at 8:53 AM, Frank E Harr

Re: [R] Inefficiency of SAS Programming

2009-02-27 Thread Thomas Levine
I had enrolled in a statistics course this semester, but after the first class, I dropped it because it uses SAS. This thread makes me quite glad. Tom! On Fri, Feb 27, 2009 at 8:48 AM, Frank E Harrell Jr wrote: > Wensui Liu wrote: >> >> Thanks for pointing me to the SAS code, Dr Harrell >> After

Re: [R] Download daily weather data

2009-02-27 Thread Thomas Levine
ere are two parts to the task: find >>patterns that allow you to pull out the datum/data you're after; and >>then write a program to pull it/them out. Also, of course, download >>the webpage (but that's no issue). >> >>I bet you'd be able to find a comm

Re: [R] Download daily weather data

2009-02-26 Thread Thomas Levine
too much work, so I just go by the current temperature in the morning, which leads to many false positives. Thomas Levine! On Thu, 2009-02-26 at 22:51 +, Barry Rowlingson wrote: > 2009/2/26 Thomas Levine : > > I'm writing a program that will tell me whether I should wear a c

[R] Download daily weather data

2009-02-26 Thread Thomas Levine
a&state=NY&site=BGM&textField1=42.4422&textField2=-76.5002&e=0&FcstType=digital), but I can't figure out how to import them. Someone must have needed to do this before. Suggestions? Thomas Levine! __ R-help@r-project.org ma

Re: [R] Putting all independent variables in one variable so performing many similar tests is more convenient

2009-02-01 Thread Thomas Levine
oops, I mean aov On Sun, Feb 1, 2009 at 9:27 AM, Thomas Levine wrote: > I want to do something like this. > > avo(q~a+b+c+d+e+f+g+h+i+j+k+l) > avo(r~a+b+c+d+e+f+g+h+i+j+k+l) > avo(s~a+b+c+d+e+f+g+h+i+j+k+l) > > (There's likely a better way to do this actually, but I th

[R] Putting all independent variables in one variable so performing many similar tests is more convenient

2009-02-01 Thread Thomas Levine
I want to do something like this. avo(q~a+b+c+d+e+f+g+h+i+j+k+l) avo(r~a+b+c+d+e+f+g+h+i+j+k+l) avo(s~a+b+c+d+e+f+g+h+i+j+k+l) (There's likely a better way to do this actually, but I think this'll work.) How do I define e=a+b+c+d+e+f+g+h+i+j+k+l such that the following works? avo(q~e) avo(r~e)

Re: [R] How should I organize data to compare differences in matchedpairs?

2008-01-25 Thread Thomas Levine
I have everything in the data frame now. The calculations I need to preform for each round are slightly different. How do I specify to make one calculation for one round and one for the other? (One round is called "C", and the other is "F") Thomas Levine On Thu, 2008-01-24 at

Re: [R] How should I organize data to compare differences in matchedpairs?

2008-01-24 Thread Thomas Levine
By accident, I didn't send this to the list. On Thu, 2008-01-24 at 17:54 -0500, Thomas Levine wrote: > Oh, right, I don't need the differences. I only needed to get the > differences before because I was doing them sloppily in a spreadsheet > and needed to do a t-test m

[R] How should I organize data to compare differences in matched pairs?

2008-01-24 Thread Thomas Levine
ll of the duplicated information that the previous method has. Anyway, I'm sure someone's done this before, so I'd like to see what other people have done for data like these. Thomas Levine [[alternative HTML version deleted]] _