[R] identical() mystery

2010-03-01 Thread Dario Strbenac
Hello, I have 2 vectors of the same mode and the same contents but I still get FALSE. Any ideas ? > reference <- c(11, 14, 16, 5, 4, 2, 0, 15, 9, 0) > reference [1] 11 14 16 5 4 2 0 15 9 0 > cpgDensity [1] 11 14 16 5 4 2 0 15 9 0 > identical(cpgDensity, reference) [1] FALSE > mode

[R] remove from mailing list

2010-03-01 Thread Bogdan Mirauta
Hello, I receive a lot of emails from the forum members.Is there any way to make the option for a daily summary?If not can may email be removed from the list?I want to continue using the forum but the numebr of emails I receive dialy is too big :) Thanks,Bogdan [[alternative HTM

[R] Smoothing Spline Basis Matrix

2010-03-01 Thread pinkdd
Anybody knows how to generate the basis matrix for smoothing spline? And how about the smoother matrix? I tried to use the following code, but there exist replicated data in X, and the length of smooth.spline(X, S[,i])$y is smaller than X, and then there is error in the last step. spl <- smoot

Re: [R] Calling SAS from R

2010-03-01 Thread Yen Lee
Thank you for your reply. Because my R program is Chinese version, I would try my best to translate and specify my question more precisely. When I type system('"c:\\Program Files\\SAS\\SAS 9.1\\sas.exe","c:\\entropy\\output7\\scale\\syntax.sas"') The warning message is as follow, In sys

Re: [R] identical() mystery

2010-03-01 Thread Søren Højsgaard
It might have to do with the storage.mode(): > reference <- c(11, 14, 16, 5, 4, 2, 0, 15, 9, 0) > storage.mode(reference) [1] "double" > cpgDensity <- as.integer(reference) > storage.mode(cpgDensity) [1] "integer" > identical(reference,cpgDensity) [1] FALSE regards Søren ___

Re: [R] remove from mailing list

2010-03-01 Thread Sharpie
mirauta wrote: > > Hello, > I receive a lot of emails from the forum members.Is there any way to make > the option for a daily summary?If not can may email be removed from the > list?I want to continue using the forum but the numebr of emails I receive > dialy is too big :) > Thanks,Bogdan > T

[R] p-values from bootstrapping of time series (tsboot)

2010-03-01 Thread Markus Troendle
Does anyone know how p-values can be generated if tsboot (stationary bootstrap) for time series is performed? That would be of great help. Thanks a lot for your comments. Markus [[alternative HTML version deleted]] __ R-help@r-project.

[R] [R-pkgs] RQDA 0.1-9 is on CRAN now

2010-03-01 Thread Wincent
Dear all, RQDA is a R package for computer-aided qualitative data analysis (CAQDA). It is an easy-to-use tool to assist in the analysis of textual data. RQDA 0.1-9 is a major upgrade with the following changes: * widgets are enabled only if they are usable. * Better code categories management sys

Re: [R] Calling SAS from R

2010-03-01 Thread Guo-Hao Huang
Did you read the help(system)? You should use ``/'', not ''\'' for you path delimiter in R. The following is a simple example for opening a pdf reader. system(paste('"C:/Program Files/Foxit Software/Foxit Reader/Foxit Reader.exe"'), wait = FALSE) It' also easy to modify the example shown in he

Re: [R] (no subject)

2010-03-01 Thread Jim Lemon
On 03/01/2010 05:00 PM, Twaha Mlwilo wrote: Hello all, Iam newbie here and to R too.Is it possible to customize R to be used to people who has limited knowledge of statistics to perform simple data analysis?.Which language to use?. Any alternative suggestion?. Hi Twaha, You may be looking

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Thu, 25 Feb 2010 17:31:19 + Patrick Burns wrote: > * What were your biggest misconceptions or > stumbling blocks to getting up and running > with R? I didn't have any major stumbling blocks, but even after years of using R I didn't have a clear concept of what exactly a vector, a list an

[R] Customizing R

2010-03-01 Thread Twaha Mlwilo
Hello all, > Iam newbie here and to R too.Is it possible to customize R to be used to > people who has limited knowledge of statistics to perform simple data > analysis?.Which language to use?. Any alternative suggestion?. > > Thank you for your help. > > udd >

Re: [R] Odp: two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Fri, 26 Feb 2010 06:44:44 -0800 (PST) Dieter Menne wrote: > Add me for reshape (the base one); Me too. Each time I try to use it, I invariably spend *hours*, and only end up with error messages ... The 'reshape' and 'plyr' packages are great, though; I use 'plyr' all the time (though someti

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Fri, 26 Feb 2010 11:56:10 -0800 (PST) Jack Siegrist wrote: > What I think would be very helpful is an introduction to programming using > R Here you are: A First Course in Statistical Programming with R http://www.cambridge.org/uk/catalogue/catalogue.asp?isbn=9780521694247 -- Karl Ove Hufth

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Mon, 1 Mar 2010 11:02:59 +0100 Karl Ove Hufthammer wrote: > > * What were your biggest misconceptions or > > stumbling blocks to getting up and running > > with R? > > Also I found it quite confusing that One more thing that still trips me up sometimes. '$' works on data frames but not on m

Re: [R] Loop overwrite and data output problems

2010-03-01 Thread RCulloch
HI Ivan, thanks for your post, I really appreciate the time you've taken over my problem! if (I==1) Sample.dat<-tmp[sample(1:max,1),] else { Sample.dat<-rbind(Sample.dat,tmp[sample(1:max,1),]) This part of the script works - I appreciate that it may not be the best option and I'm perhaps paper

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread RCulloch
Hi Petr, Thank you for your post - I really appreciate you taking the time over my problem. Apologies for not posting more data, it is just that the data set is rather large, and I don't like posting the whole thing on the website for that reason. I have managed to random sample the 98 indivi

Re: [R] two questions for R beginners

2010-03-01 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 01.03.2010 11:26:40: > On Mon, 1 Mar 2010 11:02:59 +0100 Karl Ove Hufthammer > wrote: > > > * What were your biggest misconceptions or > > > stumbling blocks to getting up and running > > > with R? > > > > Also I found it quite confusing that > > One

[R] setting the steps for x axis labels on plot

2010-03-01 Thread Gonzalo Garcia-Perate
Hello, I'm new to R, I've been working with it for the last 2 weeks. I am plotting some data and not getting the labels on the x axis I am expecting on my plot. my code reads #hours in the day h <- c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23) #hp is a data frame with a p

Re: [R] identical() mystery

2010-03-01 Thread Duncan Murdoch
Dario Strbenac wrote: Hello, I have 2 vectors of the same mode and the same contents but I still get FALSE. Any ideas ? reference <- c(11, 14, 16, 5, 4, 2, 0, 15, 9, 0) reference [1] 11 14 16 5 4 2 0 15 9 0 cpgDensity [1] 11 14 16 5 4 2 0 15 9 0 identical

Re: [R] two questions for R beginners

2010-03-01 Thread Duncan Murdoch
Karl Ove Hufthammer wrote: On Fri, 26 Feb 2010 11:56:10 -0800 (PST) Jack Siegrist wrote: What I think would be very helpful is an introduction to programming using R Here you are: A First Course in Statistical Programming with R http://www.cambridge.org/uk/catalogue/catalogue.asp?isbn

Re: [R] setting the steps for x axis labels on plot

2010-03-01 Thread Jim Lemon
On 03/01/2010 10:16 PM, Gonzalo Garcia-Perate wrote: Hello, I'm new to R, I've been working with it for the last 2 weeks. I am plotting some data and not getting the labels on the x axis I am expecting on my plot. my code reads #hours in the day h <- c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread Petr PIKAL
Hi again I believe you are fishing a wrong pond. r-help-boun...@r-project.org napsal dne 01.03.2010 11:46:12: > > Hi Petr, > > Thank you for your post - I really appreciate you taking the time over my > problem. > > Apologies for not posting more data, it is just that the data set is rathe

Re: [R] setting the steps for x axis labels on plot

2010-03-01 Thread Chuck Cleland
On 3/1/2010 6:16 AM, Gonzalo Garcia-Perate wrote: > Hello, I'm new to R, I've been working with it for the last 2 weeks. I > am plotting some data and not getting the labels on the x axis I am > expecting on my plot. > > > my code reads > > #hours in the day > h <- c(0,1,2,3,4,5,6,7,8,9,10,11,12

Re: [R] two questions for R beginners

2010-03-01 Thread Duncan Murdoch
Karl Ove Hufthammer wrote: On Mon, 1 Mar 2010 11:02:59 +0100 Karl Ove Hufthammer wrote: * What were your biggest misconceptions or stumbling blocks to getting up and running with R? Also I found it quite confusing that One more thing that still trips me up sometimes. '$' works

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread RCulloch
Hi Pter, No doubt! I have put a very short form of the data set on the email - it is basically 2 data points from each individual, which should be enough to get an idea of where I'm going wrong.hopefully! I can send this as a .csv if you prefer? Cheers, Ross SITE_NAME SITEYEAR

Re: [R] two questions for R beginners

2010-03-01 Thread Ted Harding
On 01-Mar-10 11:09:51, Petr PIKAL wrote: > Hi > r-help-boun...@r-project.org napsal dne 01.03.2010 11:26:40: >> On Mon, 1 Mar 2010 11:02:59 +0100 Karl Ove Hufthammer >> >> wrote: >> > > * What were your biggest misconceptions or >> > > stumbling blocks to getting up and running >> > > with R? >>

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Mon, 01 Mar 2010 06:37:30 -0500 Duncan Murdoch wrote: > > Some functions output matrices where you would expect them to output > > data frames, and then this problem occurs. (Is there a reason why '$' > > could/should not be made to 'work' on matrices too?) > > > The reason for the differ

Re: [R] R Experts

2010-03-01 Thread Paul Hiemstra
Hi, Maybe the R version Ryan is using is very old? cheers, Paul Erik Iverson wrote: Ryan Kinzer wrote: Erik Thanks for helping. Both of them are factors. That's the problem, they need to be of class Date. See the R NEWS article about Date classes in Volume 4/1. http://cran.r-projec

Re: [R] identical() mystery

2010-03-01 Thread jim holtman
If you use 'dput' you could provide a reproducible example, plus it would probably explain what was happening. On Mon, Mar 1, 2010 at 1:24 AM, Dario Strbenac wrote: > Hello, > > I have 2 vectors of the same mode and the same contents but I still get > FALSE. Any ideas ? > >> reference <- c(11, 1

Re: [R] Calling SAS from R

2010-03-01 Thread Uwe Ligges
On 01.03.2010 08:42, Yen Lee wrote: Thank you for your reply. Because my R program is Chinese version, I would try my best to translate and specify my question more precisely. When I type system('"c:\\Program Files\\SAS\\SAS 9.1\\sas.exe","c:\\entropy\\output7\\scale\\syntax.sas"') The

Re: [R] two questions for R beginners

2010-03-01 Thread Paul Hiemstra
Jack Siegrist wrote: My biggest impediment, as a scientist without previous programming experience, is that the R help is not beginner-friendly. I think it is probably great for experienced programmers and for the people who helped to create the software, to help them remember what they did, but

Re: [R] two questions for R beginners

2010-03-01 Thread Ted Harding
On 01-Mar-10 12:07:52, Karl Ove Hufthammer wrote: > On Mon, 01 Mar 2010 06:37:30 -0500 Duncan Murdoch > > wrote: >> > Some functions output matrices where you would expect them to output >> > data frames, and then this problem occurs. (Is there a reason why >> > '$' >> > could/should not be made

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread Petr PIKAL
The data you put before was just those 98 rows. > data=read.table("clipboard", header=T) > data$ID2 [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 [26] 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 [51] 51 52 53 54 55 56 57 58 59 60 6

[R] help with zoo

2010-03-01 Thread ravi
Hi, I am interested in decomposing a time series and getting the trend, seasonal and irregular variations, as one can get with the "stl" command. My time series is fairly regular, but it has some breaks. From the zoo manual, I gather that it should be possible to convert it to a regular time ser

Re: [R] (no subject)

2010-03-01 Thread Twaha Mlwilo
Jim, Thanks a lot, regards > Date: Mon, 1 Mar 2010 21:02:27 +1100 > From: j...@bitwrit.com.au > To: uddessy2...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] (no subject) > > On 03/01/2010 05:00 PM, Twaha Mlwilo wrote: > > > > > > > > Hello all, > > > > Iam newbie here and to R too.I

[R] MASS::loglm - exploring a collection of models with add1, drop1

2010-03-01 Thread Michael Friendly
I'd like to fit and explore a collection of hierarchical loglinear models that might range from the independence model, ~ 1 + 2 + 3 + 4 to the saturated model, ~ 1 * 2 * 3 * 4 I can use add1 starting with a baseline model or drop1 starting with the saturated model, but I can't see how to get th

[R] Random Forest prediction questions

2010-03-01 Thread Dror
Hi, I need help with the randomForest prediction. i run the folowing code: > iris.rf <- randomForest(Species ~ ., data=iris, > importance=TRUE,keep.forest=TRUE, proximity=TRUE) > pr<-predict(iris.rf,iris,predict.all=T) > iris.rf$votes[53,] setosa versicolor virginica 0.000 0.8074866

Re: [R] variable substitution in for loops

2010-03-01 Thread Jon Erik Ween
Friends First, thanks to all for great feed-back. Open-source rocks! I have a workable solution to my question, attached below in case it might be of any use to anyone. I'm sure there are more elegant ways of doing this, so any further feedback is welcome! Things I've learned (for other noobs

Re: [R] Customizing R

2010-03-01 Thread Liviu Andronic
On 3/1/10, Twaha Mlwilo wrote: > > Iam newbie here and to R too.Is it possible to customize R to be used to > > people who has limited knowledge of statistics to perform simple data > > analysis? > You might want to try one of the several R GUIs available [1]. [1] http://en.wikipedia.org/wiki

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread Petr PIKAL
Hi I finally managed to get your data from CRAN help list Here is slightly modified solutions # just for initiation value<- matrix(NA, 98, 10) # slower way rad <- nrow(data) for( i in 1:10) { # sampling row indices ss <- sample(1:rad) # ordering data according to sampled row indices Sample

[R] How to read.graph a file with strings of characters

2010-03-01 Thread Sara Khalatbari
Hi Say I have the following input *Edges a a b a cc a d a a b b b cc b d b a cc b cc cc cc d cc a d b d cc d d d in a .txt file and I want to read this graph. Which fo

[R] help

2010-03-01 Thread mahalakshmi sivamani
While loading files i am getting this type of error mentioned below. Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'D:/samplelist': No such file or directory could you please help me to resolve this problem.

[R] Public R training

2010-03-01 Thread Sarah Lewis
Mango Solutions provide high quality training courses to customers as diverse as Google and Novartis. We can provide training on specific software products or more varied training on statistics and general data analysis. As part of the course development we produce tailored documentation for

Re: [R] two questions for R beginners

2010-03-01 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 01.03.2010 13:03:24: < snip> > > > > I understand that 2 dimensional rectangular matrix looks quite > > similar to data frame however it is only a vector with dimensions. > > As such it can have items of only one type (numeric, character, ...). > > And

Re: [R] A slight trap in read.table/read.csv.

2010-03-01 Thread Peter Ehlers
On 2010-02-28 14:55, Rolf Turner wrote: I had occasion recently to read in a one-line *.csv file that looked like: "CandidateName","NSN","Ethnicity","dob","gender" "Smith, Mary Jane",111222333,"E","2/25/1989","F" That "F" (for female) in the last field got transformed to FALSE. Apparently rea

Re: [R] identical() mystery

2010-03-01 Thread Duncan Murdoch
Duncan Murdoch wrote: Dario Strbenac wrote: Hello, I have 2 vectors of the same mode and the same contents but I still get FALSE. Any ideas ? reference <- c(11, 14, 16, 5, 4, 2, 0, 15, 9, 0) reference [1] 11 14 16 5 4 2 0 15 9 0 cpgDensity

Re: [R] help with zoo

2010-03-01 Thread Gabor Grothendieck
This will give you a regular series with NAs: z1reg <- as.zooreg(as.ts(z1)) On Mon, Mar 1, 2010 at 8:03 AM, ravi wrote: > Hi, > I am interested in decomposing a time series and getting the trend, seasonal > and irregular variations, as one can get with the "stl" command. My time > series is

Re: [R] two questions for R beginners

2010-03-01 Thread Duncan Murdoch
Karl Ove Hufthammer wrote: On Mon, 01 Mar 2010 06:37:30 -0500 Duncan Murdoch wrote: Some functions output matrices where you would expect them to output data frames, and then this problem occurs. (Is there a reason why '$' could/should not be made to 'work' on matrices too?) The r

[R] Fwd: Re: Kohonen Package

2010-03-01 Thread Peter Ehlers
This was sent to me personally but was probably meant for R-help. Original Message Subject: Re: [R] Kohonen Package Date: Mon, 1 Mar 2010 14:45:13 +1000 From: Martin To: ehl...@ucalgary.ca Hi Any idea if the kohonen package can produce umatrices with hexagons and component p

Re: [R] two questions for R beginners

2010-03-01 Thread John Sorkin
If it looks like a duck and quacks like a duck, it ought to behave like a duck. To the user a matrix and a dataframe look alike . . . except a dataframe can hold non-numeric values. Thus to the users, a matrix looks like a special case of a DF, or perhaps conversely. If you can address elements

[R] Bartlett Test

2010-03-01 Thread Amit Patel
Hi I am trying to conduct a Bartlett test between two groups Samp 1 and Samp 2, both of which are vectors of equal length. I cant find any information on how to do this. Does the data need to be in a structured list. Thanks in advance __ R-help

Re: [R] Smoothing Spline Basis Matrix

2010-03-01 Thread D Sonderegger
I believe that smooth.spline fits a cubic B-spline to the data. So you just need to know the knot points (which are returned by smooth.spline as spl$fit$knot) and then use the bs() function in the splines library. pinkdd wrote: > > Anybody knows how to generate the basis matrix for smoothing s

Re: [R] Bartlett Test

2010-03-01 Thread Richardson, Patrick
?bartlett.test >From the help page "If x is a list, its elements are taken as the samples or >fitted linear models to be compared for homogeneity of variances. In this >case, the elements must either all be numeric data vectors or fitted linear >model objects, g is ignored, and one can simply u

[R] file reading /problems with encoding

2010-03-01 Thread T . Wunder
Hello, I'm a little frightened because of a problem that occured lately as I tried to read in a xml-file (for replacing some variables in the string with values from a data frame). The almost biggest problem is the encoding of the xml-file. Since it is generated by Word 2007 its encoding

Re: [R] two questions for R beginners

2010-03-01 Thread Ted Harding
On 01-Mar-10 13:57:08, Petr PIKAL wrote: > Hi > r-help-boun...@r-project.org napsal dne 01.03.2010 13:03:24: > < snip> >> > I understand that 2 dimensional rectangular matrix looks quite >> > similar to data frame however it is only a vector with dimensions. >> > As such it can have items of only o

[R] why a text editor?

2010-03-01 Thread Ivan Calandra
Dear users, From the recent discussion, I've wondered whether a text editor would be useful. However, I couldn't find a good explanation of what it is used for. For now, under Windows XP, I use the standard R Editor to write and run scripts. What can I do more with a text editor? It is wort

Re: [R] two questions for R beginners

2010-03-01 Thread Duncan Murdoch
On 01/03/2010 9:19 AM, John Sorkin wrote: If it looks like a duck and quacks like a duck, it ought to behave like a duck. To the user a matrix and a dataframe look alike . . . except a dataframe can hold non-numeric values. Thus to the users, a matrix looks like a special case of a DF, or perh

[R] Giovanna Jonalasinio è fuori ufficio, I' m away

2010-03-01 Thread Giovanna . Jonalasinio
Risposta automatica dal 1/3/10 fino al 21/3/10 I'm not in Rome Non sono a Roma [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Mon, 01 Mar 2010 09:09:11 -0500 Duncan Murdoch wrote: > >> The reason for the difference is that data.frames are lists organized > >> into columns (so the $ handling comes from the list, where it means > >> "extract the component") whereas a matrix is a single vector displayed > >> in colum

Re: [R] why a text editor?

2010-03-01 Thread Duncan Murdoch
On 01/03/2010 9:53 AM, Ivan Calandra wrote: Dear users, From the recent discussion, I've wondered whether a text editor would be useful. However, I couldn't find a good explanation of what it is used for. For now, under Windows XP, I use the standard R Editor to write and run scripts. Wha

Re: [R] why a text editor?

2010-03-01 Thread jim holtman
I use Tinn-R under Windows. It lets me keep multiple scripts open and selectively change and execute code. I usually structure a large application as a number of individual scripts that I will run from a master script using "source()". In Tinn-R, I can just highlight the statements I want execut

Re: [R] Random Forest prediction questions

2010-03-01 Thread Liaw, Andy
From: Dror > > Hi, > I need help with the randomForest prediction. i run the folowing code: > > > iris.rf <- randomForest(Species ~ ., data=iris, > > importance=TRUE,keep.forest=TRUE, proximity=TRUE) > > pr<-predict(iris.rf,iris,predict.all=T) > > iris.rf$votes[53,] > setosa versicolor virg

Re: [R] help

2010-03-01 Thread jim holtman
It can not find the file you are trying to read. You might try "choose.files()" which will bet you navigate (at least under Windows) to choose the file that you want. Check that the file is there. Use 'file.exists' to see if it exists where you think it does. On Mon, Mar 1, 2010 at 5:39 AM, mah

Re: [R] Random Forest

2010-03-01 Thread Liaw, Andy
From: Dror > > Hi, > I'm working with randomForest package and i have 2 questions: > 1. how can i drop a specific tree from the forest? Answered in another post. > 2. i'm trying to get the voting of each tree in a prediction > datum using the > folowing code > > pr<-predict(RF,NewData,type="

[R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

2010-03-01 Thread Ravi Kulkarni
Hello, I believe the aov() function in R uses a "Type-I sum-of-squares" by default as against "Type-III". This is relevant for me because I am trying to understand ANOVA in R using my knowledge of ANOVA in SPSS. I can only reproduce the results of an ANOVA done using R through SPSS if I specif

Re: [R] two questions for R beginners

2010-03-01 Thread hadley wickham
> One of the things about R which many (and that certainly includes > me) have to find out the hard way is that you have to *learn* > what to expect! You can't just import it from prior experience in > other contexts. So, by the time you have learned that a matrix > is such that all its elements mu

Re: [R] help

2010-03-01 Thread Uwe Ligges
On 01.03.2010 16:13, jim holtman wrote: It can not find the file you are trying to read. You might try "choose.files()" which will bet you navigate (at least under Windows) to choose the file that you want. Check that the file is there. Use 'file.exists' to see if it exists where you think i

[R] Fitting chi-squared distribution

2010-03-01 Thread Trafim Vanishek
Dear all, I have a question regarding performing test if the data fits chi-squared distribution. For example, using ks.test() I found in the examples how to fit it to gamma or weibull x<-rnorm(100) ks.test(x, "pweibull", shape=2,scale=1) for the gamma, pgamma can be used But I cannot find the v

Re: [R] p-values from bootstrapping of time series (tsboot)

2010-03-01 Thread Uwe Ligges
On 01.03.2010 09:59, Markus Troendle wrote: Does anyone know how p-values can be generated if tsboot (stationary bootstrap) for time series is performed? Well, under H0 we could generate n p-values simply by runif(n, 0, 1) but if you want to apply some specific test, it might make sense to

Re: [R] two questions for R beginners

2010-03-01 Thread Gustaf Rydevik
On Mon, Mar 1, 2010 at 4:02 PM, Karl Ove Hufthammer wrote: > On Mon, 01 Mar 2010 09:09:11 -0500 Duncan Murdoch > wrote: >> >> The reason for the difference is that data.frames are lists organized >> >> into columns (so the $ handling comes from the list, where it means >> >> "extract the componen

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

2010-03-01 Thread Frank E Harrell Jr
I continue to be surprised that the majority of users do not search the help archives linked from www.r-project.org. This is a much discussed topic, including vast amounts of information on the silliness of type III tests. Frank Ravi Kulkarni wrote: Hello, I believe the aov() function in

Re: [R] Fitting chi-squared distribution

2010-03-01 Thread Uwe Ligges
On 01.03.2010 16:15, Trafim Vanishek wrote: Dear all, I have a question regarding performing test if the data fits chi-squared distribution. For example, using ks.test() I found in the examples how to fit it to gamma or weibull x<-rnorm(100) ks.test(x, "pweibull", shape=2,scale=1) for the ga

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Mon, 01 Mar 2010 12:25:20 - (GMT) ted.hard...@manchester.ac.uk wrote: > > A similar type of overloading is used in the 'sp' class functions, > > where you can basically treat a 'SpatialPointsDataFrame', a > > 'SpatialLinesDataFrame' or a 'SpatialPolygonsDataFrame' as a data > > frame, >

Re: [R] file reading /problems with encoding

2010-03-01 Thread Uwe Ligges
On 01.03.2010 15:45, t.wun...@stud.uni-heidelberg.de wrote: Hello, I'm a little frightened because of a problem that occured lately as I tried to read in a xml-file (for replacing some variables in the string with values from a data frame). The almost biggest problem is the encoding of the xml

Re: [R] Odp: Loop overwrite and data output problems

2010-03-01 Thread RCulloch
Hi Petr, Thanks again for trying again with these data, I really appreciate it. Your script works perfectly, but the problem I'm having is how to store the model results so after your script I would do: m1.R<-glm(cbind(res$BEH_T, res$BEH_F) ~ res$SITE + res$YEAR + res$PRECIP_MM

Re: [R] two questions for R beginners

2010-03-01 Thread Karl Ove Hufthammer
On Mon, 01 Mar 2010 14:50:57 - (GMT) ted.hard...@manchester.ac.uk wrote: > as.character(pi) > # [1] "3.14159265358979" > > That raises a few questions about "expectations" too! Expectations can indeed be dangerous. I have been bitten by this one: as.numeric(as.character(pi)) It works

Re: [R] help with zoo

2010-03-01 Thread ravi
Hi Gabor, Thanks for your quick reply. Yes, your solution works. But not quite in the way that I wanted. For example, > length(z1) [1] 174 > length(index(z1reg)) [1] 708  Instead of 52 weeks, the indexing continues all the way from 53 to the next 1000 level, leading to lots of NA's. What I am lo

Re: [R] text editors

2010-03-01 Thread Dwayne Blind
Thanks to all of you !!! 2010/2/26 > I also agree, Emacs without question. The learning curve is a bit steep > but once you know it you can use it > for just about anything, but cleaning the kitchen sink! > > Gérald Jean > Conseiller senior en statistiques, > VP Planification et Développemen

[R] repeated measures anova, car package

2010-03-01 Thread Kay Cichini
Hello list, I' d very much appreciate some help with a two sample repeated measures ANOVA. I did the analysis yielding sign. main effects (between subj.=site, within subj.=cover) and a sign. interaction: Univariate Type II Repeated-Measures ANOVA Assuming Sphericity SS num Df Er

Re: [R] Dock graphs when using R in Eclipse with StatET plug In

2010-03-01 Thread Karl Ove Hufthammer
On Wed, 17 Feb 2010 10:32:35 -0600 Gene Leynes wrote: > By the way, if you have not tried the StatET thing, you should. It's really > nice, and I'm sure I'm not even using all the bells and whistles. The > "automatic variable highlighting" alone is amazing. I very much agree. Note that you can

[R] function odiag(): assigning values to off-diagonal

2010-03-01 Thread Benedikt Gehr
hi I'm trying to use the function odiag(x) for matrix calculations. I would like to assign new values to an off-diagonal in a matrix. When I use the diag (x) function I could write something like p<-matrix(seq(1:12),ncol=4) p.new<-matrix(rep(0,12),ncol=4) diag(p.new)<-diag(p) p.new But this w

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

2010-03-01 Thread Ravi Kulkarni
Apologies. I should have checked the archives first... I did not know that it was so widely discussed... Ravi -- View this message in context: http://n4.nabble.com/Type-I-v-s-Type-III-Sum-Of-Squares-in-ANOVA-tp1573657p1573734.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] why a text editor?

2010-03-01 Thread Joshua Wiley
Dear Ivan, One thing I would add to the current discussion is the ability of some text editors (most?) to check for matched delimiters "(...)". I also find the highlighting of syntax makes the code easier to read. In WinEdt and Tinn-R the basic features are very user friendly (I'm sure others ar

Re: [R] Bartlett Test

2010-03-01 Thread Peter Ehlers
On 2010-03-01 7:25, Amit Patel wrote: Hi I am trying to conduct a Bartlett test between two groups Samp 1 and Samp 2, both of which are vectors of equal length. I cant find any information on how to do this. Does the data need to be in a structured list. I guess you haven't looked very hard.

Re: [R] why a text editor?

2010-03-01 Thread Ivan Calandra
Thanks all for your advice. The highlighting is indeed really interesting. I'll try starting with Emacs (because I also plan to be using MacOS X) Regards, Ivan Le 3/1/2010 16:53, Joshua Wiley a écrit : > Dear Ivan, > One thing I would add to the current discussion is the ability of some > text ed

Re: [R] two questions for R beginners

2010-03-01 Thread Patrick Burns
If it looks like a duck and quacks like a duck, you ought to treat it like a duck. That is, use two subscripts: x[i, j] If you are an ornithologist, then you will know more precisely what can be done. Pat On 01/03/2010 14:19, John Sorkin wrote: If it looks like a duck and quacks like a duck

[R] help with lattice boxplots...

2010-03-01 Thread Kim Jung Hwa
Hi All, I need a small help with following code: I'm trying to convert "dashed lines" to regular ones; and changing default "blue" border color to say "black"... but I'm doing it wrong and its not working. Can anyone help please. Thanks, Code: require(lattice) ex <- data.frame(x=1:10, y=rep(c("A"

Re: [R] file reading /problems with encoding

2010-03-01 Thread T . Wunder
Hello Uwe, Quoting Uwe Ligges : You can arrange the former without a loop by xml <- paste(t, collapse="") Thank you, another trick to improve my R-skills :) For the other problem you are reporting: Can you make (the relevbant part of) your file available (say on some web site) so that we can

[R] integration function

2010-03-01 Thread li li
Hi all, I have some problem regarding the integration function. Can any one take a look at the following and give me some help? Thank you in advance! Hannah ## f0 is a function of r,x, and y f0 <- function(r, x, y){ pnorm((x-sqrt(r)*y)/sqrt(1-r), mean=0,sd=1, lower.tail=

Re: [R] Type-I v/s Type-III Sum-Of-Squares in ANOVA

2010-03-01 Thread Ista Zahn
Hi Ravi, I've posted a short explanation about this at http://yourpsyche.org/miscellaneous that you might find helpful. I'm a psychology student, not a statistician, so caveat emptor. Best, Ista On Mon, Mar 1, 2010 at 10:52 AM, Ravi Kulkarni wrote: > > Apologies. I should have checked the archiv

Re: [R] why a text editor?

2010-03-01 Thread Ted Harding
On 01-Mar-10 15:53:58, Joshua Wiley wrote: > Dear Ivan, > > One thing I would add to the current discussion is the ability > of some text editors (most?) to check for matched delimiters "(...)". > I also find the highlighting of syntax makes the code easier > to read. In WinEdt and Tinn-R the bas

Re: [R] function odiag(): assigning values to off-diagonal

2010-03-01 Thread David Winsemius
On Mar 1, 2010, at 10:53 AM, Benedikt Gehr wrote: hi I'm trying to use the function odiag(x) for matrix calculations. You assume that everyone is using the package that contains odiag? I would like to assign new values to an off-diagonal in a matrix. When I use the diag (x) function I coul

[R] question on DPpackage

2010-03-01 Thread Francesca Ieva
Hi to everyone, I'm a PhD student and I'm involved in non parametric analyses of hierarchical models. I tried to use package DPpackage on my data, but I encountered some problems in interpreting ouputs. Can anybody help me? The problem can be remued as follows: I have a logit hierarchical mode

[R] Advice wanted on using optim with both continuous and discrete par arguments...

2010-03-01 Thread Jeffrey Racine
Dear R users, I have a problem for which my objective function depends on both discrete and continuous arguments. The problem is that the number of combinations for the (multivariate) discrete arguments can become overwhelming (when it is univariate this is not an issue) hence search over the

Re: [R] Expanding a data structure

2010-03-01 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Worik R > Sent: Sunday, February 28, 2010 9:37 PM > To: r-help > Subject: Re: [R] Expanding a data structure > > > change <- c(NA, diff(q1$..1)) > > will be 1 when ..1 goes fro

Re: [R] need help to resolve RODBC error

2010-03-01 Thread Yan Zhang
Thanks for your quick response. I did install the ODBC driver. Initially the machine had unixODBC-devel-2.2.11-7.1 unixODBC-2.2.11-7.1 (both 32 and 64 bit) which I believe were loaded with RPM as part of the Linux OS, but I had problem with the ODBC drive setup, so I downloaded unixODBC-2.2.14

Re: [R] two questions for R beginners

2010-03-01 Thread hadley wickham
> Suppose X is a dataframe or a matrix.  What would you expect to get from > X[1]?  What about as.vector(X), or as.numeric(X)? > > The point is that a dataframe is a list, and a matrix isn't.  If users don't > understand that, then they'll be confused somewhere.  Making matrices more > list-like in

[R] Mainland "Equatorial Guinea" missing?

2010-03-01 Thread mooooore
In the World map, mainland "Equatorial Guinea" seems to be missing although two offshore parts, "Equatorial Guinea:Annobon" and "Equatorial Guinea:Fernando Po", are present. I'd like to be able to fill the mainland piece. How can I correct this? Thanks. -- View this message in context: http://n4

Re: [R] two questions for R beginners

2010-03-01 Thread Albert-Jan Roskam
I though duck-typing was about type-independency? I could feed the bird object bread() or carrots(), or any other method, and that's okay as long as the bird doesn't die. And since ducks don't like carrots [at least, afaik] Quaaack! ;-) Albert-Jan ~~~

  1   2   3   >