similar to: [NB] lm problems

Hi, I have a SPSS datafile which is used for my textbook in the statistics (and which is available on http://abacon.com/fox/s6720p2.sav, but it is originally from ICPSR). When I opened it with SPSS 10 and run Frequencies on it I have got 979 valid data a 27 missing. However, see below (unfortunately, I have used R in preparation of my homework, which caused me an error on this): >

Hi, how to do quickly equivalent of the following? counter = 0 for(i in 1:length(data$S2)) { if(!is.na(data$S2[i])) { counter = counter + 1 } } I have imagined something like length(x,na.rm=TRUE). How can I get values usually taken from tables like z-score, values of t distribution etc.? I could not find them among values mentioned in info file. Thanks, Matej -- Matej Cepl,

I would like to get a curve of normal distrubtion over the histogram. Something like the following (which obviously doesn't work; see attached example). maluj <- function() { vrhy=read.csv("pennies.csv",head=TRUE) hf=table(vrhy$HEADS) postscript("heads.eps",onefile=FALSE,width=4.134,height=3.445,pointsize=12) plot(hf,main="Frequency distribution of

Hi, I have this table (BTW, published by FBI and representing age distribution of the WTC tragedy victims): "LABEL" "FREQ" "MIDPOINT" "1" "Infant (under 1)" 0 "2" "1 to 4" 5 "3" "5 to 8" 1 "4" "9 to 12" 3 "5" "13 to 16" 0 "6" "17 to

I have probably misunderstood something (again :-(): > CEMPLOY<-TEMPLOY[!is.na(TEMPLOY)] > sum(CEMPLOY^2/length(CEMPLOY))-mean(CEMPLOY)^2 [1] 30877.28 > var(TEMPLOY,na.rm=TRUE) [1] 30981.25 Could anybody tell me, what's going on? I suppose that var=sum(x^2/N)-xbar^2 or what did I enter wrong? Thanks Matej -- Matej Cepl, matej at ceplovi.cz, Finger: 89EF 4BC6 288A

How do I get non-squared correlation coefficient in some more sensible way than sqrt(summary(lm(y~x))$r.squared)? Thanks Matej -- Matej Cepl, matej at ceplovi.cz, Finger: 89EF 4BC6 288A BF43 1BAB 25C3 E09F EF25 D964 84AC 138 Highland Ave. #10, Somerville, Ma 02143, (617) 623-1488 The difference between death and taxes is death doesn't get worse every time Congress meets -- Will

Hi, I have a question, how to most properly select set of cities which would be as similar as possible in some particular variables with the City of Boston (which I use as my base line). I thought about ordering cities by sum of ((differences between value of that particular variable for that particular city and the value of same variable for Boston) divided by the standard deviation of the

Hi, I've heard somewhere that it is possible to store inside .ogg files text streams and thus it should be possible to store lyrics for the sond directly in the file. Is it true? Could anybody point me towards some code how to do it (Python preferable, but not required)? Is there any documentation for text streams in .ogg files? Thanks for any reply, Mat?j -- GPG Finger: 89EF 4BC6 288A

Hi, I have got to my hands an excellent book by Michael J. Crawley ``Statistical Computing: An Introduction to Data Analysis using S-Plus'' (John Wiley & Sons, Ltd, ISBN 0-471-56040-5). Its beauty for me is in the fact, that it is more of ``An Introduction to Data Analysis'' than ``using S-Plus'', but I guess that it may be of interest for many others. Most of the

https://bugs.freedesktop.org/show_bug.cgi?id=28449 Summary: I messed up and now nouveau doesn't work. Product: xorg Version: git Platform: x86 (IA32) OS/Version: Linux (All) Status: NEW Severity: critical Priority: medium Component: Driver/nouveau AssignedTo: nouveau at lists.freedesktop.org

Hi all In two days, I am giving a small, informal workshop in our Department about using R. Initially, it was just for statistician, but surprisingly (to me anyhow) many mathematicians are also coming who have a MATLAB background. They are coming at the Workshop from a teaching perspective. They are considering using R to avoid licensing issues with MATLAB. One thing they were hoping me to

Hi, I understand that I can find t-test for the significance of the correlation coefficient on sample to population, but is there way how to make this test for more than one pair of variables together as easily as I can make a correlation matrix for more variables with cor(cbind(A,B,C,D))? Thanks, Matej -- Matej Cepl, matej at ceplovi.cz, Finger: 89EF 4BC6 288A BF43 1BAB 25C3 E09F EF25 D964

Hi, I would love to make multiple histograms transposed one on another in order to show relation between the sets. I tried to write a function like this, but R tells me, that I cannot use add=FALSE in high-level commands. That's nice but I am supposed to do? rm(list=ls()) # what's wrong with underscore? #getwd("/home/matej/docs/skola/stat\_anal-cj3534/assign01/")

is anybody whoo know something about making frequency tables in R. i'm just beginner in R. thank you. -=-=-=-=-=-=-=-=-=-=- milan toth http://toth.host.sk -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Hi, what is the minimal hardware requirement for run R on windows 9x? and for run it in a linux with X and a light windowmanager? Not for hard use, only for learning. Thanks Ronaldo -- Genius may have its limitations, but stupidity is not thus handicapped. -- Elbert Hubbard -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ]

Hi, when I try to load data from a file _in a script_ (run either with R --vanilla < script.R, or R BATCH --vanilla script.R, or even with source() from inside R), then I got this error message and scripts halts > #data.multiple <- read.table(file="multiple.csv",sep="\t") > load(file="multiple.RData") Error in open.connection(con, "rb") :

Hello List, I am trying to install RODBC package in a debian linux box but getting the following message. Can anyone help me to find what I am doing wrong here: $ R CMD INSTALL RODBC_1.0-4.tar.gz ### * Installing *source* package 'RODBC' ... checking for gcc... gcc checking for C compiler default output... a.out checking whether the C compiler works... yes checking whether we are cross

I have moved the R-announce mailing list from the Majordomo software to a combination of procmail & mailman an hour or so ago. You find a link in the footer of this message giving an URL where you can modify your subscription, look at the archives (everything from the beginning in 1997) and more. CHANGE: Postings to R-announce are now moderated, ie. they have to be approved by a list

I have moved the R-announce mailing list from the Majordomo software to a combination of procmail & mailman an hour or so ago. You find a link in the footer of this message giving an URL where you can modify your subscription, look at the archives (everything from the beginning in 1997) and more. CHANGE: Postings to R-announce are now moderated, ie. they have to be approved by a list

Dear all, my problem is following. I know Stata, but currently I have to use R. Could You please help in finding the analogy to R. (1) creating of City-Dummy. I know in Stata: g byte city=0 replace city=1 if city==12&year==2000 and (2) Create a Time-Dummy-Variable g byte T2000=0 replace T2000=1 if year==2000 (3) I need the City DUmmy for the following combination: I have the