Re: [R] Linear models interaction

gauravbhatti wrote: > > My data looks like following: > cera3[i, ] batch lcl29 pdt > Untreated 3.185867 1 0 0 > Untreated.4 3.185867 0 0 0 > LCL29 4.357552 1 1 0 > LCL29.6 3.446256 0 1 0 > PDT 2.765535

Re: [R] two questions for R beginners

Patrick Burns wrote: > > * What were your biggest misconceptions or > stumbling blocks to getting up and running > with R? > > (This derives partly from teaching) The fact that this xapply-stuff was not idempotent (worse: not always) and that you need a monster like do.call() to straighten th

Re: [R] RODBC looping sql script

Hi, You should use the paste command: for (i in 1:length(vessel[,1])) { query<-paste("Select * From vessel Where common.utility(", vessel[i,1],", ", ,vesselyear[i,2], ")", sep="") final<-sqlQuery(channel,paste(query)) } Bart Jason Gasper wrote: > > Dear R users, > > I am querying an Oracle

Re: [R] TukeyHSD troubles

treat_code isn't a factor, but a numeric variable. You should use: summary(aov(EtoH~as.factor(treat_code), mydata)) TukeyHSD(aov(EtoH~as.factor(treat_code), mydata)) Bart -- View this message in context: http://n4.nabble.com/TukeyHSD-troubles-tp1570205p1570228.html Sent from the R help mailin

Re: [R] TukeyHSD troubles

> treat_code is a dummy > variable, but that shouldn't matter.  Any suggestions? It does matter to TukeyHSD. If treat_code is a numeric variable with discrete values 0 and 1, then it does not have class "factor". It is true that aov will give the same ANOVA table for a two-level factor as for a

[R] TukeyHSD troubles

I've tried to run a Tukey post-hoc but keep getting this weird error, whether the aov was significant or not. treat_code is a dummy variable, but that shouldn't matter. Any suggestions? Thanks Amy > summary(aov(EtoH~treat_code, mydata)) Df Sum Sq Mean Sq F value Pr(>F) treat_code

Re: [R] Updating a hexbinplot

On Thu, Feb 25, 2010 at 1:24 AM, Marcin Kozak wrote: > Dear all, > > Considering this simple example of hexbinplot: > mixdata <- >     data.frame(x = c(rnorm(5000), rnorm(5000,4,1.5)), >                y = c(rnorm(5000), rnorm(5000,2,3)), >                a = gl(2, 5000)) > fig <- hexbinplot(y ~ x

[R] Legend's attribute

Hi there, I use plot(type="o") to draw the line. And I need put some legend to this line. But how can I let the legend to display the line type in the graph generated from type="o"? I can only find the lty in the legend. But I tried several times and still failed. Also, if I use abline(h=0.08) t

[R] cloud() / wireframe()

Hi All, I need some help with how to define: 1. colors in levelplot() 2. fixing labels in auto.key() using wireframe(). Sometimes I get "points" in legend and sometimes I'm getting different signs for each line like "plus sign", "star", "circle", etc... how can I be consistent with these. I'm aw

Re: [R] Ordering categories on a boxplot - a serious trap??

Hi Wilhelm, I agree it's confusing to have a levels() function that does something so different from the levels argument of the factor function. Personally I use levels() only as an extractor, never to change levels. For that I use factor(), with the levels and labels arguments as needed. Best, Is

[R] bwplot() {lattice}

Hi All, I'm trying to plot boxplot graph. I tried barchart with "groups=" option and it worked fine. But when I try to generate same kind of graph using bwplot(), "groups=" option doesn't seem to work. Though this works, yield ~ variety | site * year I'm thinking why "groups=" doesn't work in th

Re: [R] Plotting 15 million points

Hi All I should have included this first up and I think I understand the problem. The load on the server I was running R was heavy which was causing everything to slow up. >summary(s) Min. 1st Qu. MedianMean 3rd Qu.Max. 2 182 2636086 343 463 > length(s) [1]

Re: [R] S4 programming

Try this: x...@y On Thu, Feb 25, 2010 at 11:12 PM, Zhang,Yanwei wrote: > > Dear all, > > I'm new to S4 classes and have a question on this. I want to use S4 because I > want to define explicitly the slots that the new class can have. But other > than that, the new class behaves exactly like a

[R] S4 programming

Dear all, I'm new to S4 classes and have a question on this. I want to use S4 because I want to define explicitly the slots that the new class can have. But other than that, the new class behaves exactly like a list. But this will not allow me to use the generic functions that are already defi

Re: [R] Including Rcpp in Eclipse C++ project

Sam, On 25 February 2010 at 11:01, Sam Garrett wrote: | Hello. I am trying to include R in an eclipse C++ project so that I can use | Rcpp and Rinside. Right now I have R installed on my Ubuntu machine and both | plugins are installed too. I have the directories of the corresponding .so | files s

Re: [R] two questions for R beginners

On Thu, Feb 25, 2010 at 9:31 AM, Patrick Burns wrote: > * What were your biggest misconceptions or > stumbling blocks to getting up and running > with R? 1- Compared to other programming languages it is hard to learn R by example, because it is hard to find code on the web that will do the exact

Re: [R] Heterogeneous Correlation Matrix with Survey Weights

Dear Chris, As you surmised, the functions in the polycor package make no provision for survey weights. You could get polychoric correlations by using weighted contingency tables as input, but that approach won't work for polyserial correlations. Regards, John ---

Re: [R] Accessing named elements of a vector

You can convert this to list: as.list(vx)$a On Thu, Feb 25, 2010 at 6:53 PM, Matthias Gondan wrote: > Dear R developers, > > A great R feature is that elements of vectors, lists and dataframes can have > names: > > vx = c(a=1, b=2) > lx = list(a=1, b=2) > > Accessing element "a" of vx: vx['a'] >

[R] RODBC looping sql script

Dear R users, I am querying an Oracle database using sqlQuery() from the RODBC package. What I would like to do is embed my R Sql query in a for loop. For example suppose the data.frame "vessel" contains two columns: vessel[,1]= vessel id and vessel[,2]=permit year I am using "vessel" as

Re: [R] Ordering categories on a boxplot - a serious trap??

Phil, That works[*], but I still think there is a big problem given how easy it is to do the wrong thing, and that searches lead to dangerous instructions. Hopefully this will serve to keep others out of trouble, but so might an immutable return value from levels(). [*] I have not yet done an

Re: [R] Building R packages in Windows 7

On 25/02/2010 7:56 PM, Eric Ferreira wrote: Thank you, Sir, but how can I demand it to create HTML files? The tools::Rd2HTML function will do the translation, but it takes a bit of work to prepare input for it. The idea is that you don't need to save the files, R just produces them when the br

Re: [R] R and Wordpress

I have found a partial solution to this. Wordpress will only chew the first <- on a given post. Therefore, I add a comment with a <-. e.g. # <- hello <- c(1,2,3) WP will only break the first <-. On Thu, Feb 25, 2010 at 9:33 PM, Steve Lianoglou wrote: > Hi, > > On Wed, Feb 24, 2010 at 11:40 PM

[R] two questions for R beginners

* What were your biggest misconceptions or stumbling blocks to getting up and running with R? * What documents helped you the most in this initial phase? I especially want to hear from people who are lazy and impatient. Feel free to write to me off-list. Definitely write off-list if you are ju

Re: [R] Plotting 15 million points

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Abhishek Pratap > Sent: Thursday, February 25, 2010 3:12 PM > To: r-help@r-project.org > Subject: [R] Plotting 15 million points > > Hi All > > I have a vector of about 15 millio

Re: [R] Subset Question

Thank you Eric and Petr. It seems to be working now! -- View this message in context: http://n4.nabble.com/Subset-Question-tp1568555p1569461.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://st

Re: [R] Subset Question

Thank you Eric and Petr. It seems to be working now! -- View this message in context: http://n4.nabble.com/Subset-Question-tp1568555p1569464.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://st

Re: [R] Ordering categories on a boxplot - a serious trap??

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Schwab,Wilhelm K > Sent: Thursday, February 25, 2010 3:51 PM > To: r-help@r-project.org > Subject: [R] Ordering categories on a boxplot - a serious trap?? > > Hello all, > > I t

Re: [R] R ANOVA gives diferent results than SPSS

Thanks all of you. the main problem I had then was it was not recognizing one of my variables as a grouping variable, it was treating its numerical values as data, not group designations. Thanks for the help in clarifying. I may be a regualr in calling for help on these R sited until I can get th

[R] Ordering categories on a boxplot - a serious trap??

Hello all, I think I probably did something stupid, and R's part was to allow me to do it. My goal was to control the order of factor levels appearing horizontally on a boxplot. Enter search engines and perhaps some creative stupidity on my part, and I came up with the following: v=r

Re: [R] Plotting 15 million points

On Feb 25, 2010, at 6:11 PM, Abhishek Pratap wrote: Hi All I have a vector of about 15 million numbers which I would like to plot. The goal is the see the distribution. I tired the usual steps. I get that way after a long day myself. 1. Histogram : never gets complete my window freeze

Re: [R] two questions for R beginners

My biggest blocker was my misconception that R is extremely difficult to start with. It is powerful and one can do very complicated things ( that consequently turn things complicated) but it comes with very nice defaults and one can produce great results with standard tasks in very little time - e

Re: [R] Plotting 15 million points

Have you considered taking a random subset and plotting that? I'd bet you can get a really impression of the distribution with a few hundred thousand points at most. Tim Glover Senior Environmental Scientist - Geochemistry Geoscience Department Atlanta Area MACTEC Engineering and Consulting,

Re: [R] Plotting 15 million points

On 2/25/10, Abhishek Pratap wrote: > Any efficient methods to play with > 10 million numbers in a vector. > Did you try rggobi? Liviu __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide htt

Re: [R] two questions for R beginners

On Thu, Feb 25, 2010 at 5:39 PM, Carl Witthoft wrote: > Well, here goes... > > I still wish there were a really good monograph on the use and > implementation of factors. To get a good handle on factors, and the sets of contrasts they encode, it is really necessary to study a good statistics book

[R] Plotting 15 million points

Hi All I have a vector of about 15 million numbers which I would like to plot. The goal is the see the distribution. I tired the usual steps. 1. Histogram : never gets complete my window freezes w/out log base 10 2. Density : I first calculated the kernel density and then plotted it which worke

Re: [R] two questions for R beginners

My biggest stumbling blocks to getting up and running with R was whenever I was lazy and impatient. The more you love R, the more it loves you back. Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 R

Re: [R] two questions for R beginners

Patrick Burns wrote: > > * What were your biggest misconceptions or > stumbling blocks to getting up and running > with R? > R was the first scripting language that I *really* invested time in learning. Prior to R I had a few years experience programming in Fortran and had worked on a few pro

[R] Heterogeneous Correlation Matrix with Survey Weights

Hello, I have a data set containing categorical and ordinal factors, as well as sampling weights (i.e., survey weights reflecting unequal probabilities of selection). I want to fit a structural equation model with sem(). I have run sem() on weighted covariance matrices using advice from John F

Re: [R] two questions for R beginners

> "The best way to spread information is to tell someone that it is a secret, > the best way to keep it secret is to put it in  > a manual." ==> Nice quote. ;-) The problem is not that there's too little information, rather there's so much. That's probably because R is so powerful, but it makes

Re: [R] two questions for R beginners

Well, here goes... I still wish there were a really good monograph on the use and implementation of factors. I had to do a certain amount of digging to learn that {assign, get, eval, expression, call, parse, deparse} all existed and how they play together. Sometimes they are look like the C

Re: [R] Restructure some data

Harold - Here's what I came up with: tapply(as.vector(as.matrix(dat[5:7])), + list(rep(dat$id,3),as.vector(as.matrix(dat[2:4]))),I) item1 item10 item2 item3 item4 item5 item7 item9 1NA NA 1NANA 1NA 0 2 0 NANANANA 1 1N

[R] missing values in 'ncdf', integer versus byte

I'm working with the ncdf package, and I'm seeing a behavior which I am unsure whether it is a bug or a feature. I am creating a variable, which I make with the precision of "byte" and also set a missing value, When I put into this variable, all my NA values are put in properly as the missing val

Re: [R] Alternatives to linear regression with multiple variables

On Feb 22, 2010, at 7:46 AM, Guy Green wrote: I wonder if someone can give some pointers on alternatives to linear regression (e.g. Loess) when dealing with multiple variables. Taking any simple table with three variables, you can very easily get the intercept and coefficients with:

Re: [R] two questions for R beginners

Patrick Burns wrote: * What were your biggest misconceptions or stumbling blocks to getting up and running with R? * What documents helped you the most in this initial phase? I especially want to hear from people who are lazy and impatient. Can't be bothered with questionnaires and can't wait

[R] Accessing named elements of a vector

Dear R developers, A great R feature is that elements of vectors, lists and dataframes can have names: vx = c(a=1, b=2) lx = list(a=1, b=2) Accessing element "a" of vx: vx['a'] Accessing element "a" of lx: lx[['a']] or lx$a Might be a matter of taste, but I like the $ very much. Unfortunatel

[R] Including Rcpp in Eclipse C++ project

Hello. I am trying to include R in an eclipse C++ project so that I can use Rcpp and Rinside. Right now I have R installed on my Ubuntu machine and both plugins are installed too. I have the directories of the corresponding .so files set in the Project > Properties > Settings Gcc c++ linker under L

Re: [R] taking the median across similar data

Hi Hari, Here is a suggestion: with(x, aggregate(x[,-1], list(probename), median)) where x is your data set. HTH, Jorge On Thu, Feb 25, 2010 at 4:25 PM, Harikrishnadhar <> wrote: > Dear All, > > I am analyzing the miRNA data set in which I have 817 unique probes for > each > they have 20 fe

Re: [R] Normal distribution (Lillie.test())

On Feb 25, 2010, at 12:20 PM, Bosken wrote: Hi, Thanks for your reaction. The purpose of my test is to check if my NoiseGenerators really are Normal Distributed en witch circuit is the best! So I need some good test to do this. But what with: Fortune(), can't find anything about it..

[R] taking the median across similar data

Dear All, I am analyzing the miRNA data set in which I have 817 unique probes for each they have 20 features each . I have to group the similar features and take the median across them so that I have a data with no repeats to perform invariant analysis . My data looks something similar format

Re: [R] How to add a title to represent four different plot in lm function

Hi there, Try this: par(mfrow=c(2,2)) for( i in 1:4 ) plot(1:10) mtext("Title",side=3,outer=TRUE,padj=3) Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biost

[R] Restructure some data

Suppose I have a data frame like "dat" below. For some context, this is the format that represents student's taking a computer adaptive test. first.item is the first item that student was administered and then score.1 is the student's response to that item and so forth. item.pool <- paste("item

Re: [R] read.table (trying one more time)

On Feb 25, 2010, at 3:10 PM, Mestat wrote: Hey listers, I've posted this doubt other day... But I am still having a hard time... I have a MAC and I am not getting how to use the read.table command... At the bottom of the FINDER application I have the following path: Macintosh HD>Usuários>U

Re: [R] How to fill in a region with different patterns?

Hi, If you are curious you might like to try a highly experimental Grid function I wrote some time ago, library(grid) source("http://gridextra.googlecode.com/svn/trunk/R/patternGrob.r";) grid.newpage() grid.pattern(x=seq(1/6, 5/6, length=6), width=unit(1/8,"npc"), height=unit(0.5,"npc"),

[R] Linear models interaction

My data looks like following: cera3[i, ] batch lcl29 pdt Untreated 3.185867 1 0 0 Untreated.4 3.185867 0 0 0 LCL29 4.357552 1 1 0 LCL29.6 3.446256 0 1 0 PDT 2.765535 1 0 1 PDT.5 3.584963

Re: [R] Alternatives to linear regression with multiple variables

Well, the help page for the loess function says that the formula can include up to 4 predictor variables. There are also additive models (mgcv or gam (or other) package). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -

Re: [R] Odp: How to read percentage and currency data?

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Hrishi Mittal > Sent: Wednesday, February 24, 2010 8:56 AM > To: r-help@r-project.org > Subject: Re: [R] Odp: How to read percentage and currency data? [snip] > It's a shame tha

[R] Help with simple bootstrap test

Hi all Forgive me, I'm a total R newbie, and this seems to be a straightforward simple bootstrap problem, but after a whole day of trying to figure out how to do it I'm ready to give up. Part of the problem is that every example and every help page seems to be about doing something more far more c

Re: [R] Restructure some data

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Doran, Harold > Sent: Thursday, February 25, 2010 10:35 AM > To: r-help@r-project.org > Subject: [R] Restructure some data > > Suppose I have a data frame like "dat" below. For s

Re: [R] Normal distribution (Lillie.test())

Hi, Thanks for your reaction. The purpose of my test is to check if my NoiseGenerators really are Normal Distributed en witch circuit is the best! So I need some good test to do this. But what with: Fortune(117) and fortune(234), can't find anything about it.. Thanks for the help! Bosken --

Re: [R] Inverse function

?uniroot -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of li li > Sent: Wednesday, February 24, 2010 2:31 PM >

Re: [R] Single character input without Enter

Greg Snow-2 wrote: > > The playSudoku function in the sudoku package has 2 examples/methods of > responding to single key strokes, whether those methods will work for your > application or not depends on what you are trying to do. > > Thanks, Greg, it boils down to the rather basic getGraphi

Re: [R] How to fill in a region with different patterns?

The rect function will draw a rectangle and can fill it with diagonal lines at specified angle and density. But before you do that, really consider if and why you want to. Using the diagonal lines became popular when the only way to get quality graphics was with a pen plotter (mechanical arm a

Re: [R] How to add a title to represent four different plot in lm function

Try something like: > par(oma=c(0,0,3,0)) > mtext("your text here", outer=TRUE Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-b

Re: [R] two questions for R beginners

Patrick, I would add one more question: * where did you look for help expecting answers, but did not find them? If you add hubris to laziness and impatience, you have Larry Wall's 3 virtues of a programmer. To new users of R who may not understand why Patrick is asking: Patrick Burns is the a

[R] Minimum Spanning Trees

Hi, I need to find all minimum spanning trees of an unweighted graph. Is there a way in R to do that? Thanks -- View this message in context: http://n4.nabble.com/Minimum-Spanning-Trees-tp1569351p1569351.html Sent from the R help mailing list archive at Nabble.com. __

[R] nested design

This is a nested design where team is nested in group. Method is fixed, while team and group are random. Is this the correct R code for this? gunload.aov <- aov(rounds ~ method+group + Error((team %in% group)/method),data=gunload, qr=T) -- View this message in context: http://n4.nabble.com/nest

Re: [R] two questions for R beginners

On 2/25/10, Patrick Burns wrote: > * What were your biggest misconceptions or > stumbling blocks to getting up and running > with R? > > * What documents helped you the most in this > initial phase? > > I especially want to hear from people who are > lazy and impatient. > I'm quite resilient

Re: [R] Help with simple bootstrap test

The boot() function in the 'boot' package expects to find a function for the statistic with two arguments: The data object plus a row index object. You don't indicate enough to see how you will be resampling. It you sum all elements in your table, resampling would have to be one of: 1. A samp

[R] read.table (trying one more time)

Hey listers, I've posted this doubt other day... But I am still having a hard time... I have a MAC and I am not getting how to use the read.table command... At the bottom of the FINDER application I have the following path: Macintosh HD>Usuários>User>Marcio>UdeM>Travail Dirigé>Data>MU284 Populati

Re: [R] Normal distribution (Lillie.test())

May be you should have said: "Normality tests are standard farce in a host of statistical texts." Ravi. --- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine

Re: [R] Normal distribution (Lillie.test())

Install and load the fortunes package first, then run fortune(117), etc. Then run fortune() quite a few times for possible enlightenment (or at least mild entertainment). Do your NoiseGenerotors need to generate exactly normal data (they don't, see SnowsPenultimateNormalityTest), or is there a

Re: [R] Single character input without Enter

The getGraphicsEvent function only works on windows (at least last time I checked), so if you are working only on windows, you can use that. If you want something that works cross platform, then the tcltk solution is a possibility. Your original question did not indicate if you would be willing

Re: [R] two questions for R beginners

Apparently I need to explain the "lazy and impatient" comment. No offence was intended (quite the contrary). The meaning of it is that the higher your level of frustration, the more valuable your comments are likely to be to me. On 25/02/2010 17:31, Patrick Burns wrote: * What were your bigges

Re: [R] angle phase mapping

Look at the my.symbols function in the TeachingDemos package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of

Re: [R] two questions for R beginners

I started using statistical software with the commercial product S+ when I obtained a new HP735 workstation. We kept the S+ license going for a number of years until I heard about R. It was an easy transition and because I have been proficient in fortran and perl, the scripting came naturally

Re: [R] Bimodal distribution

You could use: library(TeachingDemos) x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) TeachingDemos:::vis.test( x, TeachingDemos:::vt.normhist ) Then click on the plot that looks most different. Sorry about needing TeachingDemos::: for now, the functions were accidentally left out of the NAMESPACE fil

Re: [R] Minimum Spanning Trees

Hi, On Thu, Feb 25, 2010 at 12:16 PM, harunpirim wrote: > > Hi, > > I need to find all minimum spanning trees of an unweighted graph. > Is there a way in R to do that? If you're working with graphs in R, you'll probably want to see what the igraph library has to offer. For instance, it has a "mi

[R] error in lmLists in lme4 package (bug?)

Hello, I am trying to use lmLists in the lme4 package and copying over very standard code from the nlme package given in 'Mixed-Effects Models in S and S-Plus'. It appears to not accept an 'I(age-11)' in the formula, though it will accept the formula with out the subtraction of 11 from age. Th

Re: [R] Reading a password-protected Excel workbook

Use RDCOMClient or rcom to open it and then use RODBC to read it. On Thu, Feb 25, 2010 at 7:26 AM, Dimitri Shvorob wrote: > > I try reading a password-protected spreadsheet with RODBC, passing a file or > an ODBC source to odbcConnect, e.g. > > con = odbcConnect(dsn = "mywkbk", uid = "", pwd = "m

Re: [R] Building R packages in Windows 7

On 25/02/2010 11:49 AM, Eric Ferreira wrote: Ok, I'm working under: Windows 7 Professional 32bits, 4 GB RAM, 320 GB HD, Intel Core 2 Duo processor R 2.10.1 I've installed: Rtools211 MikteX 2.8 HTML Help Workshop Setting my PATH to: c:\Rtools\bin;c:\Rtools\perl\bin;c:\Rtools\MinGW\bin;c:\Arquiv

Re: [R] Normal distribution (Lillie.test())

You should probably read fortune(117) and fortune(234) (and possibly some of the original discussions that lead to the fortunes). Reading the help page for the SnowsPenultimateNormalityTest function (TeachingDemos package) may also help. If you are happy with the plots, but still feel the need

Re: [R] Use 2 "ifelse" to sort data

There is the %<% operator in the TeachingDemos package that allows comparisons like: 5 %<% a %<% 10 But in this case it is not needed (and in this case it takes more key strokes than 5 < a & a < 5, so may not be worth it). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

Re: [R] Odp: export results

If your data can fit in memory, consider creating a list of the intermediate values and then 'cbind'ing the result into a matrix that you want to write out. result <- lapply(1:100, function(x) your function ) # if everything is the same size result <- do.call(cbind, result) On Thu, Feb 2

Re: [R] data mining

It depends on what you want to do. Start by looking at: http://cran.r-project.org/web/views/ There are also some nice startup tutorials on google, start from that... Tal Contact Details:---

Re: [R] How to extract one of four plots in a linear regression model

Hi, You can extract a plot in using the option which in specifying the number of the plot (from 1 to 6). For example: > plot(lm.D9, which=1) Regards, Alain Guillet On 25-Feb-10 16:50, FMH wrote: Dear All, A linear regression model could be fitted by using lm function and the plot functio

Re: [R] Building R packages in Windows 7

On 25/02/2010 10:56 AM, Eric Ferreira wrote: This is my first package. I'm just getting started doing that, following the steps described on you website... I really don't know how I asking for CHMs to be produced, sorry. All I can suggest is that you need to be less stingy with information.

[R] How to add a title to represent four different plot in lm function

Dear All, A linear regression model could be fitted by using lm function and the plot function can be used to check the assumptions of the model. The example is as followed. require(graphics) ## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data.

[R] How to do: Correlation with "blocks" (or - "repeated measures" ?!) ?

Hello dear R help group, I have the following setup to analyse: We have about 150 subjects, and for each subject we performed a pair of tests (under different conditions) 18 times. The 18 different conditions of the test are complementary, in such a way so that if we where to average over the test

Re: [R] Building R packages in Windows 7

On 25/02/2010 10:40 AM, Eric Ferreira wrote: Dear Duncan Thank so much for your reply. Actually, I'm using the latest version of R and the problem persists. What do you use instead of HTML Help Workshop for newer R versions? We just produce text and HTML help pages on demand, and LaTeX ones

[R] How to extract one of four plots in a linear regression model

Dear All, A linear regression model could be fitted by using lm function and the plot function can be used to check the assumption of the model. The help menu shows few instances on suitable coding for fitting such a linear model. In addition, four different plots could be extracted simultaneo

[R] Odp: export results

Hi r-help-boun...@r-project.org napsal dne 25.02.2010 12:42:46: > > Hi all, > > I am looping through a function for 100 time in the middle of my code, and I > want to output the results from this function. Is there a way to write the > results into a txt or csv fil? For example, I write the r

Re: [R] mtext with at=NULL

On 25/02/2010 8:58 AM, Ulrike Feldmann wrote: Hello, in the help of mtext I found "at: If ‘length(at)==0’ (the default), the location will be determined by ‘adj’" But if I use mtext( "hello world", at=NULL) there comes the following error message: Fehler in mtext("hallo Welt", side = 3, adj =

Re: [R] error using pvcm() on unbalanced panel data

On 2/25/10, Liviu Andronic wrote: > I can reproduce the error on a modified example from the vignette: > > require(plm) > > data("Hedonic") > > Hed <- pvcm(mv ~ crim + zn + indus + chas + nox + rm + age + dis +rad + > tax + ptratio + blacks + lstat, Hedonic, model = "within",index = "townid")

[R] export results

Hi all, I am looping through a function for 100 time in the middle of my code, and I want to output the results from this function. Is there a way to write the results into a txt or csv fil? For example, I write the results from the first loop to the first column of a spreadsheet and the results

[R] mtext with at=NULL

Hello, in the help of mtext I found "at: If ‘length(at)==0’ (the default), the location will be determined by ‘adj’" But if I use mtext( "hello world", at=NULL) there comes the following error message: Fehler in mtext("hallo Welt", side = 3, adj = NA, at = NULL) : 'at' mit Länge 0 angegeben

[R] Reading a password-protected Excel workbook

I try reading a password-protected spreadsheet with RODBC, passing a file or an ODBC source to odbcConnect, e.g. con = odbcConnect(dsn = "mywkbk", uid = "", pwd = "mypwd") but get "Could not decrypt file" pop-up error message. Can anyone help? Thank you. -- View this message in context: htt

[R] error using pvcm() on unbalanced panel data

Dear all I am trying to fit Variable Coefficients Models on Unbalanced Panel Data. I managed to fit such models on balanced panel data (the example from the "plm" vignette), but I failed to do so on my real, unbalanced panel data. I can reproduce the error on a modified example from the vignette:

Re: [R] behavior of seq_along

On Feb 25, 2010, at 9:02 AM, Dale Steele wrote: I'm trying to understand the behavior of seq_along in the following example: x <- 1:5; sum(x) y <- 6:10; sum(y) data <- c(x,y) S <- sum( data[seq_along(x)] ) S T <- sum( data[seq_along(y)] ) T Why is T != sum(y) ? Look at seq_along(y) seq

Re: [R] Building R packages in Windows 7

On 25/02/2010 9:06 AM, Eric Ferreira wrote: Dear useRs, I'm having trouble building R packages in Windows 7 regarding HTML help Workshop. Pointing PATH to c:\Program Files\HTML help Workshop does work in Windows (e.g. Vista) and does not in Windows 7. Some tips?? We don't use the HTML Help

Re: [R] reducing data.frame

Perhaps the reshape package? It's just about impossible to read your data layout. Could you resubmit the example using dput()? Thanks --- On Thu, 2/25/10, AC Del Re wrote: > From: AC Del Re > Subject: [R] reducing data.frame > To: r-help@r-project.org > Received: Thursday, February 25, 2010

  1   2   >