Re: [R] manova: R vs SAS...need some clarification

Please see the footer of this message: we need to know what you did. Also, SAS may have made some assumptions for you without telling you (for example used a numerically ill-conditioned covariance matrix), and we don't know what you did in SAS, either. On Tue, 12 Aug 2008, Pedro Mardones

Re: [R] Update Package on CRAN

stephen sefick wrote: To update a package on CRAN I just update all of the version information stuff etc. and then upload it to the ftp site? Stephen Sefick Yes, just build the package and submit it to CRAN as before, but with an increased version number. Uwe Ligges

Re: [R] Installing R in Ubuntu

On Tue, Aug 12, 2008 at 9:24 PM, Shreyasee Pradhan [EMAIL PROTECTED] wrote: Hi, I am running Ubuntu on my Windows OS through VMware. I am trying to install R in Ubuntu, but not getting with those commands, which are there on the site. Can anyone please tell me how to install it, stepwise,

Re: [R] aligned memory allocation in C

On Tue, 12 Aug 2008, Jeffrey Horner wrote: Christophe Dutang1 wrote: Hi, I'm currently R porting SF Mersenne Twister algorithm of Matsumoto and Saito. To get the full power of their code, I want to use their fonction fill_array32 which need aligned memory. That is to say I need to use the C

Re: [R] Installing R in Ubuntu

Hi, Thanks for that. the way I tried is as follows: 1) Downloaded the r-base package 2) went in that directory where the r-base package was downloaded from command line 3) entered the command sudo apt-get install r-base But got the error, that Couldn't find r-base command. I don't

Re: [R] Installing R in Ubuntu

Shreyasee Pradhan [EMAIL PROTECTED] writes: Hi, Thanks for that. the way I tried is as follows: 1) Downloaded the r-base package 2) went in that directory where the r-base package was downloaded from command line 3) entered the command sudo apt-get install r-base But got the

Re: [R] sqlQuery with date attribute

Hi Many GetReturn-function(code,date) { db-C:/Test.mdb channel-odbcConnectAccess(db) ssql-paste(select * from tblCalendarDate Where CalendarID =,code,and DateRebal =,date) print(ssql)# so as i can see what ssql contains mydata-sqlQuery(channel,ssql) mydata } [snip] This is the content

Re: [R] sqlQuery with date attribute

Thank you for your answer. Actually, I've tried with this function where I added the # symbol between the date: GetReturn-function(code,date) { db-C:/Test.mdb channel-odbcConnectAccess(db) ssql-paste(select * from tblCalendarDate Where CalendarID =,code,and DateRebal= #,date,#) print(ssql)

Re: [R] summary.manova rank deficiency error + data

Pedro Mardones wrote: Dear R-users; Previously I posted a question about the problem of rank deficiency in summary.manova. As somebody suggested, I'm attaching a small part of the data set. #*** test - structure(.Data = list(structure(.Data =

Re: [R] Senging commands to the GUI in Windows through a script

Cool, many thanks Henrik. Tolga Henrik Bengtsson [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 13/08/2008 02:03 To [EMAIL PROTECTED] cc Prof Brian Ripley [EMAIL PROTECTED], r-help@r-project.org Subject Re: [R] Senging commands to the GUI in Windows through a script With AutoIt

Re: [R] Installing R in Ubuntu

Hi, If you download a package to your harddrive for installation you need to use the dpkg command like: 1) Download pacakge (foo.deb) 2) Go to the directory 3) dpkg -i foo.deb But I would advise against this because it is better to use repositories so R get updated automatically. The

[R] mob(party) formula question

I try tu use mob() with my data.frame ('data.frame':288 obs. of 81 variables; factors, numerics and ordered factors) My response is a binary variable and I should use for modelling a logistic regression (family=binomial). I read in the MOB Vignette that I could use a formula like this if I

[R] Calculating an appropriate error from the NNET package for a continuous target

I have been using the NNET package and have successfully run neural networks on both continuous and binary targets. I managed to search the internet and found out how to capture the error resulting from a binary model no problems. My problem now is that I am trying to find how to calculate an

Re: [R] aligned memory allocation in C

Yes, it seems a good idea but your two questions are also good questions! Le 13 août 08 à 04:52, Jeffrey Horner a écrit : Christophe Dutang1 wrote: Hi, I'm currently R porting SF Mersenne Twister algorithm of Matsumoto and Saito. To get the full power of their code, I want to use their

Re: [R] dixon test

Hi, thank you very much for your useful help =). just a question...I don't know what is the distribution of my data (normal, T, etc...). So, how can I set the type parameter? There is a type value to use in case of a distribution-free statistical test? Thank you so much! Fernando

Re: [R] fPortfolio constraints, maxsumW

JPB == John P. Burkett [EMAIL PROTECTED] on Tue, 12 Aug 2008 10:46:28 -0400 JPB Running R version 2.6.1 under Gentoo Linux and using the fPortfolio JPB package, I am having trouble specifying a sector constraint. One of the JPB constraints to be imposed is that assets 1 and 2

[R] problems with packages tseries and robustbase

Dear R Users, Is there a known problem with downloading packages robustbase and tseries from the UK CRAN website ? Thanks in advance, Tolga = R version 2.7.1 (2008-06-23) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN

[R] Help me: nls and try function

Dear All, I have these problems: 1) How can use the function try in nls model: try(nls(...)) 2) I have 100 colun with data and I want ro prepare 99 file with the first colun with the others Time A1 A2 A3 A4 AN. I want to have 99 files with a)Time and A1 b)Time and A2 n) Time AN thanks

Re: [R] problems with *downloading* packages tseries and robustbase

It works from my home ISP (Virgin Media) and from .ox.ac.uk, so I think the problem is local to you, perhaps your DNS server. Ask you IT support for erm ... support. Please do try to use an accurate subject line (see the posting guide). (Why don't peoople just try a different mirror or some

Re: [R] problems with *downloading* packages tseries and robustbase

Hi, It may be a problem with the mirror and my location. However, I should have added that the following packages had no problems installing from that mirror and location, which is why I suspected it was something more. RBloomberg strucchange car lmtest nlme corrgram RODBC MSBVAR xtable vars

Re: [R] Help me: nls and try function

For the first question, you have provided the answer -- try(nls(...)). Was there something else you wanted? For part 2, this should work: for (i in names(myData)[-1]){ # skip first column with Time write.table(myData[, c(Time, i)], file=i) } 1) How can use the function try in nls model:

Re: [R] Location of HTML help files [on Firefox 3]

I've added a couple of workarounds to this in R-patched and hence the upcoming R 2.7.2. 1) There is a new menu item on Rgui to go directly to SearchEngine.html. 2) help.start() has a new argument searchEngine=TRUE to do the same. R 2.7.2 is 12 days' away, so it would be appreciated if Firefox

Re: [R] dixon test

giov [EMAIL PROTECTED] 13/08/2008 10:59:32 just a question...I don't know what is the distribution of my data (normal, T, etc...). So, how can I set the type parameter? You must assume an underlying distribution or you can't do an outlier test. Outliers are just unusually extreme data

Re: [R] Threshold vector error correction models

Hello I worked on threshold cointegration for my master thesis and wrote code for R. This code will be published in a next release of package tsDyn, when I will have time to finish it (there is a difference between using its own code and making it avalaible on R... I did'nt realize there is

[R] rgl/compiz problem

I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom in and out with the mouse wheel but the mouse buttons do nothing. I can't rotate the plot. This has been mentioned and discussed here and in other lists before, and the solution is to

Re: [R] mob(party) formula question (example)

Here is an example that produces the same error: Read in the following as textfile (save as DFExample.txt): 1 2 3 4 7 8 9 10 12 13 14 15 16 17 18 19 21 22 23 25 27 28

Re: [R] rgl/compiz problem

Two days ago I installed compiz on my Debian laptop. It plays fine with the OpenGL games that I also have on that computer. (My son plays the games in the brief interludes between my intense R hacking sessions. I, of course, have no time for such frivolity. The production cycle is sleep - eat -

Re: [R] mob(party) formula question

On Wed, 13 Aug 2008, Birgitle wrote: I try tu use mob() with my data.frame ('data.frame':288 obs. of 81 variables; factors, numerics and ordered factors) My response is a binary variable and I should use for modelling a logistic regression (family=binomial). I read in the MOB Vignette

Re: [R] aligned memory allocation in C

On Wed, 13 Aug 2008, Christophe Dutang1 wrote: Hi, I'm currently R porting SF Mersenne Twister algorithm of Matsumoto and Saito. To get the full power of their code, I want to use their fonction fill_array32 which need aligned memory. That is to say I need to use the C function memalign on

[R] R freezes on Xeon multiprocessor and Win 64-bit

Hi everybody, Performing a stepAIC on a glm.nb object, from a database of more than 10,000 records and about 50 independent variables, on a 64-bit workstation with two Intel Xeon 3.20Ghz processors (keeping the HyperThreading option disabled in the BIOS), using 4 out of 7Gb available RAM,

[R] merging data sets to match data to date

Hi everyone, I want to extract data from a data set according to dates specified in a vector. I have created a blank matrix with row names (dates) that I want to extract from the full data set. I have then performed a merge to try to o/p rows corresponding to common dates to a results matrix,

[R] Comination of two barcharts and one xyplot

Hi Rhelpers, I would like to have some help with a plot which is beyond my capabilities. This plot that I am seeking involves an overlay of two different barcharts and one xyplot. The code that I have used is the following : #save(df1,file=M:\\KBR\\df1.RData) load(file=M:\\KBR\\df1.RData) #

Re: [R] merging data sets to match data to date

rcoder wrote: Hi everyone, I want to extract data from a data set according to dates specified in a vector. I have created a blank matrix with row names (dates) that I want to extract from the full data set. I have then performed a merge to try to o/p rows corresponding to common dates to a

Re: [R] merging data sets to match data to date

Try this: x - data.frame(Dates = seq(as.Date('2008-01-01'), as.Date('2008-01-31'), by = 'days'), Values = sample(31)) subset(x, Dates %in% as.Date(c('2008-01-05', '2008-01-20'))) On 8/13/08, rcoder [EMAIL PROTECTED] wrote: Hi

Re: [R] merging data sets to match data to date

zoo and merge.zoo- read the help files. Use chron to generate a list of dates that correspond to the one that you want, and then merge away. This should get you started Stephen Sefick On Wed, Aug 13, 2008 at 9:08 AM, Erik Iverson [EMAIL PROTECTED] wrote: rcoder wrote: Hi everyone, I want

Re: [R] Comination of two barcharts and one xyplot

not reproducible On Wed, Aug 13, 2008 at 9:07 AM, ravi [EMAIL PROTECTED] wrote: Hi Rhelpers, I would like to have some help with a plot which is beyond my capabilities. This plot that I am seeking involves an overlay of two different barcharts and one xyplot. The code that I have used is

Re: [R] rgl/compiz problem

Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom in and out with the mouse wheel but the mouse buttons do nothing. I can't rotate the plot. This has been mentioned and

Re: [R] mob(party) formula question

Many thanks for your answer and the code that you offered me. I get this error message after calling mob (look at my given example). I guess it has something to do with the missings? The iris example works also fine for me. Sorry that I am not enough into statistics to really understand the

[R] Tiny help for tiny function

I just started to write tiny functions and therefore I appologise in advance if I am asking stupid question. I wrote a tiny function to give me back from the original matrix, a matrix showing only the values smaller -0.8 and bigger 0.8. y-c(0.1,0.2,0.3,-0.8,-0.4,0.9)

[R] The standard deviation of measurement 1 with respect to measurement 2

Hi, I have two (different types of) measurements, say X and Y, resulting from the same set of experiments. So X and Y are paired: (x_1, y_1), (x_2, y_2), ... I am trying to calculate the standard deviation of Y with respect to X. In other words, in terms of the scatter plot of X and Y, I would

Re: [R] Tiny help for tiny function

You can do this: ifelse(XY = 0.8 | XY = -0.8 | is.na(XY), XY, low corr) On 8/13/08, Birgitle [EMAIL PROTECTED] wrote: I just started to write tiny functions and therefore I appologise in advance if I am asking stupid question. I wrote a tiny function to give me back from the original

Re: [R] Tiny help for tiny function

Many thanks. Much easier than my solution B. Birgitle wrote: I just started to write tiny functions and therefore I appologise in advance if I am asking stupid question. I wrote a tiny function to give me back from the original matrix, a matrix showing only the values smaller -0.8 and

Re: [R] mob(party) formula question

On Wed, 13 Aug 2008, Birgitle wrote: Many thanks for your answer and the code that you offered me. I get this error message after calling mob (look at my given example). I guess it has something to do with the missings? Yes, you have to handle NAs in advance if you want to fit that model.

Re: [R] summary.manova rank deficiency error + data

Thanks for the reply. The SAS output is attached but seems to me that doesn't correspond to the wihtin-row contrasts as you suggested. By the way, yes the data are highly correlated, in fact each row correspond to the first part of a signal vector. Thanks anyway PM

Re: [R] mob(party) formula question

Thanks again. Unfortunately I have always this missing values problem. But the missings have also a meaning and its impossible to code it differently or impute. Also thanks for the explanation. Now I understand. B. Achim Zeileis wrote: On Wed, 13 Aug 2008, Birgitle wrote: Many thanks

Re: [R] mob(party) formula question

On Wed, 13 Aug 2008, Birgitle wrote: Thanks again. Unfortunately I have always this missing values problem. But the missings have also a meaning and its impossible to code it differently or impute. That's ok. Just to clarify: NAs are not allowed in the response or the modeling variables. In

Re: [R] Comination of two barcharts and one xyplot

Hi Rhelpers, Thanks a lot, Stephen, for showing me the way to get a data frame into a pasteable format with the dput command. My code is given below with the new correction. This should work, as Stephen says, right off the bat :-) ## df1 is the first data frame df1 -structure(list(Year =

Re: [R] ignoring zeros or converting to NA

Hi, since many suggestions are following the form of x[x==0] (or similar) I would like to ask if this is really recommended? What I have learned (the hard way) is that one should not test for equality of floating point numbers (which is the default for R's numeric values, right?) since the

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

Dear Peter and Henrik, Thanks for your replies - this helps speed up a bit, but I thought there would be something much faster. What I mean is that I thought that a particular value of a level could be accessed instantly, similarly to a hash key. Since I've got about 6000 levels in that data

[R] bmp header

Hi R users, I have a xml file. A value of one of the nodes of the xml file is a bmp image(RAW format) encoded in base64. I would like to read this image by R. I think I should do the following steps: 1. Decoding it from base64 to binary. 2. Removing the header of the image file 3. building the

Re: [R] rgl/compiz problem

2008/8/13 Ben Bolker [EMAIL PROTECTED]: Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom in and out with the mouse wheel but the mouse buttons do nothing. I can't rotate the

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

I still don't understand what you are doing. Can you make a small example that shows what you have and what you want? Is ?split what you are after? Emmanuel Levy wrote: Dear Peter and Henrik, Thanks for your replies - this helps speed up a bit, but I thought there would be something much

Re: [R] ignoring zeros or converting to NA

Integers (up to a fairly high limit) are represented exactly, as are fractions whose denominator is a power of two (again up to a fairly high limit), so x==0 is fine in that sense. If x is computed by floating point operations you do have to worry whether these are exact, eg, with

Re: [R] rgl/compiz problem

Barry Rowlingson wrote: 2008/8/13 Ben Bolker [EMAIL PROTECTED]: Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: oe So he tried mouse buttons in combination. Holding B1 and then B3 and moving the mouse resulted in a zoom operation. Holding first B3 and then B1 resulted in rotation

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

Sorry for being unclear, I thought the example above was clear enough. I have a data frame of the form: name info 1 YAL001C 1 2 YAL001C 1 3 YAL001C 1 4 YAL001C 1 5 YAL001C 0 6 YAL001C 1 7 YAL001C 1 8 YAL001C 1 9 YAL001C 1 10 YAL001C 1

[R] Flag variable

Hi All, I have 4000 case which have string variables in them, i want to do some fuzzy matching and create a new variable that is of the same length with 0 or 1's if i use the code test- agrep(web Klick,ETC$Exposure.Type , max = 2, ignore.case = TRUE) it works but i get length(test) [1] 3127

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

Wow great! Split was exactly what was needed. It takes about 1 second for the whole operation :D Thanks again - I can't believe I never used this function in the past. All the best, Emmanuel 2008/8/13 Erik Iverson [EMAIL PROTECTED]: I still don't understand what you are doing. Can you make

Re: [R] summary.manova rank deficiency error + data

Pedro Mardones wrote: Thanks for the reply. The SAS output is attached but seems to me that doesn't correspond to the wihtin-row contrasts as you suggested. By the way, yes the data are highly correlated, in fact each row correspond to the first part of a signal vector. Thanks anyway PM

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

split if probably what you are after. Here is an example: n - 270 x - data.frame(name=sample(1:6000,n,TRUE), value=runif(n)) # split it into 6000 lists system.time(y - split(x$value, x$name)) user system elapsed 0.800.201.07 str(y[1:10]) List of 10 $ 1 : num [1:454]

Re: [R] which(df$name==A) takes ~1 second! (df is very large), but can it be speeded up?

If you want the index, then use: system.time(y - split(seq(nrow(x)), x$name)) user system elapsed 0.810.060.88 str(y[1:10]) List of 10 $ 1 : int [1:454] 6924 17503 26880 39197 42881 50835 57896 62624 65767 75359 ... $ 2 : int [1:440] 9954 25619 25761 33776 56651 60372 61042

Re: [R] ignoring zeros or converting to NA

The help page on binary operators (see ?==) confirms that binary representation of fractional representation is not catered for and points to all.equal as a more suitable test method for those cases. Steve E Thomas Lumley [EMAIL PROTECTED] 13/08/2008 16:47 Integers (up to a fairly high limit)

[R] tcl/tk example in batch

The example for learning tcl/tk under R at http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/OKtoplevel.html suggests running it from batch - but when I do, the window flashes by and the example ends. I'm under XP pro. Is there a workaround? Should I create a modal window instead so it

[R] Arguments to lm() within a function - object not found

Hi all, I'm having some difficulty passing arguments into lm() from within a function, and I was hoping someone wiser in the ways of R could tell me what I'm doing wrong. I have the following: lmwrap - function(...) { wts - somefunction() print(wts) # This works, wts has the values I expect

[R] reverse orientation of text in plot margins

Dear R users, I am trying to reverse the orientation of axis labels and title in the right margin of a plot, so that they read from top to bottom. I know that this can be done using text() as follows: par(mar=c(5,4,4,4)+.1) plot(1:4,las=0) par(new=T) y - rnorm(4)

[R] conditional IF with AND

Hi everyone, I'm trying to create an if conditional statement with two conditions, whereby the statement is true when condition 1 AND condition 2 are met: code structure: if ?AND? (a[x,y] condition1, a[x,y] condition2) I've trawled through the help files, but I cannot find an example of the

Re: [R] dixon test

Thank you so much, I have not much experience on outliers =), I thought that there were nonparametric distribution-free outliers test =(. What is the most general distribution I can use? I did histogram of my data set and sometimes normal distribution seems to occur, sometimes an uniform

Re: [R] subsetting matrix according to columns with character index

Try this: x V1 V2 V3 1 a1 c1 1 2 a1 c1 2 3 a2 c1 1 4 a1 c2 1 5 a1 c2 2 lis - split(x, list(x$V1, x$V2), drop = TRUE) do.call(rbind, unname(lis[sapply(lis, function(x)all(1:2 %in% x[,3]))])) On Wed, Aug 13, 2008 at 3:00 PM, Ralph S. [EMAIL PROTECTED] wrote: Hi, I have a long matrix

Re: [R] conditional IF with AND

if(cond1 cond2) { ... } rcoder wrote: Hi everyone, I'm trying to create an if conditional statement with two conditions, whereby the statement is true when condition 1 AND condition 2 are met: code structure: if ?AND? (a[x,y] condition1, a[x,y] condition2) I've trawled through the help

[R] Arguments to lm() within a function - object not found

Hi all, I'm having some difficulty passing arguments into lm() from within a function, and I was hoping someone wiser in the ways of R could tell me what I'm doing wrong. I have the following: lmwrap - function(...) { wts - somefunction() print(wts) # This works, wts has the values I expect

Re: [R] conditional IF with AND

See: ?`` On Wed, Aug 13, 2008 at 1:45 PM, rcoder [EMAIL PROTECTED] wrote: Hi everyone, I'm trying to create an if conditional statement with two conditions, whereby the statement is true when condition 1 AND condition 2 are met: code structure: if ?AND? (a[x,y] condition1, a[x,y]

[R] How to test for a random effect in a repeated measu res analysis using anova.mlm ?

Dear “R” masters, I am trying to conduct an ANOVA with repeated measures using the command anova.mlm for data structured according to a Randomized Block Design. I would like to account for a random effect but cannot find a way to incorporate it in the analysis. NB. I tried using the argument “M”

Re: [R] re placing default labels in lattice

You can see the source code of demo script: file.show(system.file(demo/labels.R, package = lattice)) On Wed, Aug 13, 2008 at 11:20 AM, Andrewjohnclose [EMAIL PROTECTED] wrote: Dear all, I am having a little trouble deciphering how to change the default x-axis labels in a lattice xyplot (or

Re: [R] ignoring zeros or converting to NA

FYI, there is an isZero() in the R.utils package that allows you to specify the precision. It looks like this: isZero - function (x, neps=1, eps=.Machine$double.eps, ...) { (abs(x) neps*eps); } /Henrik On Wed, Aug 13, 2008 at 8:23 AM, Roland Rau [EMAIL PROTECTED] wrote: Hi, since many

Re: [R] conditional IF with AND

On 13-Aug-08 16:45:27, rcoder wrote: Hi everyone, I'm trying to create an if conditional statement with two conditions, whereby the statement is true when condition 1 AND condition 2 are met: code structure: if ?AND? (a[x,y] condition1, a[x,y] condition2) I've trawled through the help

Re: [R] which alternative tests instead of AIC/BIC for choosing models

Many thanks John, appreciate the advice, Tolga John C Frain [EMAIL PROTECTED] 13/08/2008 18:51 To [EMAIL PROTECTED] cc r-help@r-project.org Subject Re: [R] which alternative tests instead of AIC/BIC for choosing models My initial idea would be to forget about AIC and BIC, ask the

Re: [R] Arguments to lm() within a function - object not found

On Wed, 13 Aug 2008, Pete Berlin wrote: Hi all, I'm having some difficulty passing arguments into lm() from within a function, and I was hoping someone wiser in the ways of R could tell me what I'm doing wrong. I have the following: lmwrap - function(...) { wts - somefunction() print(wts)

[R] need help with stat functions(like adaboost, random forests and glm)

Ok, so basically I have a dataframe named data_frame data_frame contains: startdate startprice endpricethreshold1 endpricethreshold2 endpricethreshold3 all of these endpricethresholds are true/false binary vectors. They are true or false depending on whether the endprice was above or below

Re: [R] subsetting matrix according to columns with character index

I tried this - I get an empty set: 0 rows (or 0-length row.names) I guess this happens because the z variable takes only one value per row?? What works is: DFsub-DF[DF$z == 1 | DF$z == 2,] but then, I do not eliminate the entries where there is only one entry for z given an a and c

Re: [R] subsetting matrix according to columns with character index

i don't think i understood what you were trying to do, atleast based on Henrique's solution which I haven't cut and pasted yet in order to understand. Did Henrique's solution do what you wanted ? On Wed, Aug 13, 2008 at 2:45 PM, Ralph S. wrote: I tried this - I get an empty set: 0 rows (or

Re: [R] which alternative tests instead of AIC/BIC for choosing models

By way of partial follow-up to my own question, and on the odd chance anyone else wonders about this issue, some alternatives to this appear to be in the leaps package, which implements the leaps routine (Mallows Cp) and regsubsets. In my case Mallows' Cp does not work either (see below), so I

Re: [R] subsetting matrix according to columns with character index

sorry ralph. i meant the OR instead of the AND so that was my bad mistake. the subset function should also work with the OR. i think i understand better what you want now also. the approach below for doing what you want assumes that , if there are 2 rows associated with the values in the

Re: [R] subsetting matrix according to columns with character index

Ralph: I looked at Henrique's solution and he does 2 things which make it better than mine. 1) He splits based off the first two columns where I just split based on the second. So, my split assumes that the same rows are next to each other which is an unnecessary assumption. 2) He actually

Re: [R] merging data sets to match data to date

Dear Henrique, This is exactly what I need. Thank you very much for your help! rcoder Henrique Dallazuanna wrote: Try this: x - data.frame(Dates = seq(as.Date('2008-01-01'), as.Date('2008-01-31'), by = 'days'), Values

Re: [R] Arguments to lm() within a function - object not found

Thanks very much for the quick reply. I had looked at the help for lm, but I clearly skimmed over the critical part explaining where weights is evaluated. Thanks, Pete On 13/8/2008, Prof Brian Ripley wrote: On Wed, 13 Aug 2008, Pete Berlin wrote: Hi all, I'm having some difficulty passing

Re: [R] conditional IF with AND

Thank you all for your replies. This is all very useful information for me! Ted, thank you very much for the extra explanation and example. Many thanks, rcoder Ted.Harding-2 wrote: On 13-Aug-08 16:45:27, rcoder wrote: Hi everyone, I'm trying to create an if conditional statement with

Re: [R] which alternative tests instead of AIC/BIC for choosingmodels

your model 3 is the unrestricted model and your models 1 and 2 are restricted models. you can test model 1 and 2 against model 3 using the anova function, e.g. anova(model2,model3), which, for the case of OLS estimation, are compared with an F-test. If the test is insignificant, the simpler model

Re: [R] reverse orientation of text in plot margins

On Wed, 13-Aug-2008 at 06:32PM +0200, Karel Van den Meersche wrote: | | Dear R users, | | I am trying to reverse the orientation of axis labels and title in | the right margin of a plot, so that they read from top to bottom. I | know that this can be done using text() as follows: |

Re: [R] which alternative tests instead of AIC/BIC for choosing models

Cp is either the same thing as AIC, or an approximation to it. So it is not an 'alternative'. See e.g. the discussion in MASS or ?add1. On Wed, 13 Aug 2008, [EMAIL PROTECTED] wrote: By way of partial follow-up to my own question, and on the odd chance anyone else wonders about this issue,

Re: [R] The standard deviation of measurement 1 with respec t to measurement 2

Firas Swidan frsswdn at gmail.com writes: Hi, I have two (different types of) measurements, say X and Y, resulting from the same set of experiments. So X and Y are paired: (x_1, y_1), (x_2, y_2), ... I am trying to calculate the standard deviation of Y with respect to X. In other

Re: [R] which alternative tests instead of AIC/BIC for choosing models

Dear R Users, I am looking for an alternative to AIC or BIC to choose model parameters. This is somewhat of a general statistics question, but I ask it in this forum as I am looking for a R solution. Suppose I have one dependent variable, y, and two independent variables, x1 an

[R] change 3x3 cell size in 1x1 cell size

Hi All, I wish to change 3x3 pixel size in 1x1 pixel size my grid. I have this fuction: dem.area - ([EMAIL PROTECTED],[EMAIL PROTECTED],1])*([EMAIL PROTECTED],[EMAIL PROTECTED],1]) dem.pixelsize - round(5*sqrt(dem.area/length(ground$Z)),0) dem.pixelsize where is the input to change?

[R] Conditional statement used in sapply()

Hi, I have data stored in a list that I would like to aggregate and perform some basic stats. However, I would like to apply conditional statements so that not all the data are used. Basically, I want to get a specific variable, do some basic functions (such as a mean), but only get the data

Re: [R] Conditional statement used in sapply()

Hi Mark, How about this? result - sapply(split(res, res$Volume0)$`TRUE`, mean) There is one thing I'm not sure: is res$Volume a vector or single numeric? -gary -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Altaweel, Mark R. Sent: Wednesday, August

Re: [R] Conditional statement used in sapply()

-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Altaweel, Mark R. Sent: Wednesday, August 13, 2008 3:03 PM To: r-help@r-project.org Subject: [R] Conditional statement used in sapply() Hi, I have data stored in a list that I would like to

Re: [R] Conditional statement used in sapply()

Hi, Yes, that's it. I got the correct results. Thanks everyone for their help once again. This is a great help board. Mark -Original Message- From: Steven McKinney [mailto:[EMAIL PROTECTED] Sent: Wed 8/13/2008 5:29 PM To: Altaweel, Mark R.; r-help@r-project.org Subject: RE: [R]

Re: [R] rgl/compiz problem

Barry Rowlingson wrote: I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom in and out with the mouse wheel but the mouse buttons do nothing. I can't rotate the plot. This has been mentioned and discussed here and in other lists before,

Re: [R] rgl/compiz problem

on 08/13/2008 06:03 PM Duncan Murdoch wrote: Barry Rowlingson wrote: I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom in and out with the mouse wheel but the mouse buttons do nothing. I can't rotate the plot. This has been mentioned

Re: [R] Conditional statement used in sapply()

Hello - Altaweel, Mark R. wrote: Hi, I have data stored in a list that I would like to aggregate and perform some basic stats. However, I would like to apply conditional statements so that not all the data are used. Basically, I want to get a specific variable, do some basic functions (such

Re: [R] Comination of two barcharts and one xyplot

At 01:17 14/08/2008, you wrote: Hi Rhelpers, Thanks a lot, Stephen, for showing me the way to get a data frame into a pasteable format with the dput command. My code is given below with the new correction. This should work, as Stephen says, right off the bat :-) ## df1 is the first data frame

Re: [R] rgl/compiz problem

My laptop has an nVidia card. Maybe that's why it works? Simon. On Wed, 2008-08-13 at 13:17 +, Ben Bolker wrote: Barry Rowlingson b.rowlingson at lancaster.ac.uk writes: I have just encountered the problem with rgl where plot3d figures don't interact with the mouse. My plots zoom

  1   2   >