[R] Support Counting

2011-04-03 Thread psombe
Hi, I'm new to R and trying to some simple analysis. I have a data set with about 88000 transactions and i want to perform a simple support count analysis of an itemset which is say not a complete transaction but a subset of a transaction. say {A,B,D} is a transaction and i want to find support

Re: [R] loading R object files on an RWeb server

2011-04-03 Thread Dieter Menne
jose romero-3 wrote: > > For one thing, i am using google docs to host the R object file and google > docs has secure https URL's, which apparently cannot be handled by R's > url(). So my questions are these: > > Try ?getURL in the RCurl package Dieter -- View this message in context: htt

[R] RmetricsTools.R

2011-04-03 Thread Luis Felipe Parra
Hello I read on the Rmetrics webpage that all the development packages could be intalled using the following command source("RmetricsTools.R") install.RmetricsDev("[Rmetrics package to be installed] I would like to know where I could get this RmetricsTools.R . I suppose it might be somewhere

Re: [R] detect filetype (as in unix 'file')

2011-04-03 Thread Prof Brian Ripley
On Sun, 3 Apr 2011, Jeroen Ooms wrote: Is there a way in R (in Linux) to detect the type of a file without invoking a shell? E.g to do this: system("file density.plot") density.plot: PDF document, version 1.4 but without using system()? I tried file() and file.info(), but both do display the

Re: [R] Structural equation modeling in R(lavaan,sem)

2011-04-03 Thread Jeremy Miles
On 3 April 2011 12:38, jouba wrote: > > Daer all, > I have a question concerning longitudinal data: > When we have a longitudinal data and we have to do sem analysis there is in > the package lavaan some functions,options in this package that help to do > this or we can treat these data like non

Re: [R] replace last 3 characters of string

2011-04-03 Thread jim holtman
Will this do it for you: > x <- c('asdfasdf', 'sadfasdf/:/', 'sadf', 'asdf/:/') > sub("/:/$", '', x) [1] "asdfasdf" "sadfasdf" "sadf" "asdf" On Sun, Apr 3, 2011 at 7:39 PM, Bert Jacobs wrote: > Hi, > > > > I would like to replace the last tree characters of the values of a certain > column i

Re: [R] Ordering every row of a matrix while ignoring off diagonal elements

2011-04-03 Thread jim holtman
I assume that this is what you did, and I would not call that cheating; it is just a reasonable way to solve the problem: > x <- as.matrix(read.table(textConnection(" 0.00 2.384158 2.0065682 > 2.2998856 + 2.384158 0.00 1.4599928 2.4333213 + 2.006568 1.459993 0.000 0.9733285 + 2.299886

Re: [R] another question on shapefiles and geom_point in ggplot2

2011-04-03 Thread Pierre Roudier
Hi all, 2011/4/4 Felipe Carrillo : > Manuel: > As far as I know one needs gpclibPermit() in order to fortify > see this: > Note: polygon geometry computations in maptools >     depend on the package gpclib, which has a >     restricted licence. It is disabled by default; >     to enabl

[R] Ordering every row of a matrix while ignoring off diagonal elements

2011-04-03 Thread Vivian Shih
Sorry if this is a stupid question but I've been stuck on how to code so that I can order rows of a matrix without paying attention to the diagonal elements. Say, for example, I have a matrix d: d [,1] [,2] [,3] [,4] [1,] 0.00 2.384158 2.0065682 2.2998856 [2,] 2.3

[R] detect filetype (as in unix 'file')

2011-04-03 Thread Jeroen Ooms
Is there a way in R (in Linux) to detect the type of a file without invoking a shell? E.g to do this: > system("file density.plot") density.plot: PDF document, version 1.4 but without using system()? I tried file() and file.info(), but both do display the information I am looking for. -- View

[R] loading R object files on an RWeb server

2011-04-03 Thread jose romero
Hello list: I have some R code/data sets that i'd like to make available to others through RWeb server use.  The code/data were saved via "save.image" and should be made available by the following command sequence: con <- url("https:// ") load(con) close(con) However, this is where the pro

[R] RGtk2: How to populate an GtkListStore data model?

2011-04-03 Thread Cleber N. Borges
hello all I am trying to learn how to use the RGtk2 package... so, my first problem is: I don't get the right way for populate my gtkListStore object! any help is welcome... because I am trying several day to mount the code... Thanks in advanced Cleber N. Borges --- # my

Re: [R] How do I modify uniroot function to return .0001 if error ?

2011-04-03 Thread Hans W Borchers
eric aol.com> writes: > > I am calling the uniroot function from inside another function using these > lines (last two lines of the function) : > > d <- uniroot(k, c(.001, 250), tol=.05) > return(d$root) > > The problem is that on occasion there's a problem with the values I'm > passing to uni

Re: [R] Graph many points without hiding some

2011-04-03 Thread Samuel Dennis
Thanks heaps for all your help, sorry for the late reply, I'm a bit overwhelmed by all the possibilities here! My variables have different lengths which makes automated randomisation difficult Peter, although not impossible. Thanks for the suggestion of 3D graphs Nick, I'm not sure how I can make

[R] LIVE, ONLINE COURSE: Using R Software for Academic Research Analyses

2011-04-03 Thread Geoffrey Hubona
The Information Institute (http://www.information-institute.org), a charitable, non-profit, educational and scientific organization, in conjunction with faculty from Virginia Commonwealth University (VCU), is offering live, interactive, synchronous online classes (one AM and one PM to reach all tim

Re: [R] Linear Model with curve fitting parameter?

2011-04-03 Thread stephen sefick
Steven: You are exactly right sorry I was confused. ### so log(y-intercept)+log(K) is a constant called b0 (is this right?) lm(log(Q)~log(A)+log(R)+log(S)-1) is fitting the model log(Q)=a*log(A)+r*log(R)+s*log(S) (no beta 0) and lm(log(Q)~

Re: [R] replace last 3 characters of string

2011-04-03 Thread Bert Jacobs
Thx I could imagine it was so simple.:) Bert -Original Message- From: Dirk Eddelbuettel [mailto:e...@master.debian.org] On Behalf Of Dirk Eddelbuettel Sent: 04 April 2011 01:39 To: Bert Jacobs Cc: r-help@r-project.org Subject: Re: [R] replace last 3 characters of string On Mon, Apr 04,

Re: [R] replace last 3 characters of string

2011-04-03 Thread Dirk Eddelbuettel
On Mon, Apr 04, 2011 at 01:39:34AM +0200, Bert Jacobs wrote: > I would like to replace the last tree characters of the values of a certain > column in a dataframe. > > This replacement should only take place if the last three characters > correspond to the value "/:/" and they should be replaced w

[R] replace last 3 characters of string

2011-04-03 Thread Bert Jacobs
Hi, I would like to replace the last tree characters of the values of a certain column in a dataframe. This replacement should only take place if the last three characters correspond to the value "/:/" and they should be replaced with ""(blank) I cannot perform a simple gsub because the chara

[R] power of 2 way ANOVA with interaction

2011-04-03 Thread Timothy Spier
I've been searching for an answer to this for a while but no joy. I have a simple 2-way ANOVA with an interaction. I'd like to determine the power of this test for each factor (factor A, factor B, and the A*B interaction). How can I do this in R? I used to do this with "proc Glmpower" in SAS, b

Re: [R] Quick question about duplicating vectors

2011-04-03 Thread Xiaoxi Gao
Sorry for the imprecision. I haven't used r for a long time... > CC: r-help@r-project.org > From: dwinsem...@comcast.net > To: rhel...@hotmail.com > Subject: Re: [R] Quick question about duplicating vectors > Date: Sun, 3 Apr 2011 18:17:55 -0400 > > > On Apr 3, 2011, at 6:04 PM, Xiaoxi Gao wro

Re: [R] another question on shapefiles and geom_point in ggplot2

2011-04-03 Thread Felipe Carrillo
Manuel: As far as I know one needs gpclibPermit() in order to fortify see this: Note: polygon geometry computations in maptools     depend on the package gpclib, which has a     restricted licence. It is disabled by default;     to enable gpclib, type gpclibPermit()  I am going to guess

Re: [R] R-project: plot 2 zoo objects (price series) that have some date mis-matches

2011-04-03 Thread Gabor Grothendieck
On Sun, Apr 3, 2011 at 1:57 PM, algotr8der wrote: > I have 2 zoo objects - > > 1) Interest rate spread between 10-YR-US-Treasury and 2-YR-US-Treasury > (object name = sprd) > > 2) S&P 500 index (object name = spy) > > >> str(spy) > ‘zoo’ series from 1976-06-01 to 2011-03-31 >  Data: num [1:8791] 9

Re: [R] zoo:rollapply by multiple grouping factors

2011-04-03 Thread Gabor Grothendieck
On Sun, Apr 3, 2011 at 11:58 AM, Mark Novak wrote: > # Hi there, > # I am trying to apply a function over a moving-window for a large number of > multivariate time-series that are grouped in a nested set of factors.  I > have spent a few days searching for solutions with no luck, so any > suggesti

Re: [R] Quick question about duplicating vectors

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 6:04 PM, Xiaoxi Gao wrote: Hello R users, I have a quick question, if I have a data set a, say a <- c(1,2,3,4,5) You might as well learn to use more precise R terminology... that is a `vector`. and I want to get a new data set b, b <- c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5)

[R] Avoiding loops in creating a coinvestment matrix

2011-04-03 Thread Daniel Malter
Hi, I am working on a dataset in which a number of venture capitalists invest in a number of firms. What I am creating is an asymmetric matrix M in which m(ij) is the volume (sum) of coinvestments of VC i with VC j (i.e., how much has VC i invested in companies that VC j also has investments in). T

Re: [R] style question

2011-04-03 Thread David Winsemius
On May 1, 2008, at 1:46 PM, Sebastián Daza wrote: Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berke

Re: [R] another question on shapefiles and geom_point in ggplot2

2011-04-03 Thread Manuel Spínola
Thank you very much Felipe, Did you see the solution from ahmadou dicko? He doesn´t use gpclibPermit() I have another option but I cannot get the right fill for the id. See attached map. ai_biotica = readOGR(dsn="C:/ProyectosRespacial/ICE/SIG_Biotica_PHED", layer="AI_BIOTICA_010411_CRTM05")

[R] Quick question about duplicating vectors

2011-04-03 Thread Xiaoxi Gao
Hello R users, I have a quick question, if I have a data set a, say a <- c(1,2,3,4,5) and I want to get a new data set b, b <- c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5) Could anyone tell me how I can obtain this? Thank you! Xiaoxi [[alternative HTML ver

Re: [R] Function for finding NA's

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 3:46 PM, Tyler Rinker wrote: aThanks David, After seeing the simplicity of your function versus the convoluted mess I worked up I now understand why it's not necessary to have a package to find NA's (and from what you said is a part of other packages such as Hmisc alre

Re: [R] installing ncdf on Ubuntu 10.04

2011-04-03 Thread Thomas Lumley
On Mon, Apr 4, 2011 at 6:35 AM, Steve Friedman wrote: > Hello > > I'm trying to install the "ncdf" package on a Ubuntu 10.04 laptop. > > Can you offer suggestions to install this package? > checking for netcdf.h... no > configure: error: netcdf header netcdf.h not found > ERROR: configuration f

Re: [R] Plotting MDS (multidimensional scaling)

2011-04-03 Thread Daniel Malter
The header "metric mds" was actually a leftover because I initially used cmdscale and did not bother changing it for the example. Thanks, Daniel -- View this message in context: http://r.789695.n4.nabble.com/Plotting-MDS-multidimensional-scaling-tp3422670p3424135.html Sent from the R help mailin

[R] How do I modify uniroot function to return .0001 if error ?

2011-04-03 Thread eric
I am calling the uniroot function from inside another function using these lines (last two lines of the function) : d <- uniroot(k, c(.001, 250), tol=.05) return(d$root) The problem is that on occasion there's a problem with the values I'm passing to uniroot. In those instances uniroot stops and

Re: [R] Function for finding NA's

2011-04-03 Thread Tyler Rinker
aThanks David, After seeing the simplicity of your function versus the convoluted mess I worked up I now understand why it's not necessary to have a package to find NA's (and from what you said is a part of other packages such as Hmisc already). I am at the 2 1/2 month mark as an R user a

[R] R-project: plot 2 zoo objects (price series) that have some date mis-matches

2011-04-03 Thread algotr8der
I have 2 zoo objects - 1) Interest rate spread between 10-YR-US-Treasury and 2-YR-US-Treasury (object name = sprd) 2) S&P 500 index (object name = spy) > str(spy) ‘zoo’ series from 1976-06-01 to 2011-03-31 Data: num [1:8791] 99.8 100.2 100.1 99.2 98.6 ... Index: Class 'Date' num [1:8791] 2

Re: [R] Structural equation modeling in R(lavaan,sem)

2011-04-03 Thread jouba
Daer all, I have a question concerning longitudinal data: When we have a longitudinal data and we have to do sem analysis there is in the package lavaan some functions,options in this package that help to do this or we can treat these data like non longitudinal data Thanks you a lot Antr

[R] Standard Error for Cointegration Results

2011-04-03 Thread Kin Ming Wong
Dear Sir/Madam, I have used ca.jo in urca package to identify the cointegration and cajorls to estimate the vecm. Althought both return the coefficients for long run relationship (or ect1 in cajorls), I am unable to find the standard error and t statistics. I spend some weeks to search arou

[R] installing ncdf on Ubuntu 10.04

2011-04-03 Thread Steve Friedman
Hello I'm trying to install the "ncdf" package on a Ubuntu 10.04 laptop. Can you offer suggestions to install this package? > install.packages("ncdf", repositories = TRUE) Installing package(s) into ‘/home/steve/R/i486-pc-linux-gnu-library/2.12’ (as ‘lib’ is unspecified) --- Please select a CRA

[R] Homals package color function problem

2011-04-03 Thread dsg
Hello The Homals package and its plot options are excellent. However, I am unable to manipulate the colour in the plots. In a call such as: plot(mc_analysis, plot.dim = c(1,3), plot.type = "jointplot", col = 1) this should be straightforward - but I can't seem to effect the plotted colours I h

[R] style question

2011-04-03 Thread Sebastián Daza
Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Reje

[R] coefficients style

2011-04-03 Thread Sebastián Daza
Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Reje

[R] setCoefTemplate

2011-04-03 Thread Sebastián Daza
Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Reje

Re: [R] Discretizing data rows into regular intervals

2011-04-03 Thread wulei wong
i suggest you'd better translate the.data format to another form like in SAS, the data can be a number count from a fixed time like 2005-01-01, then you can translate the datas to a sequence. On 03-Apr-2011 5:40 PM, "Gabor Grothendieck" wrote: On Sat, Apr 2, 2011 at 9:31 PM, Linh Tran wrote: >

Re: [R] R gui on windows how to force to always show the last line of output

2011-04-03 Thread wulei wong
i wonder why you want to change R but not your program itself, you just print the last sentence is ok? wuleiwong statistic@CSU On 03-Apr-2011 9:06 PM, "lcn" wrote: CTRL+L helps. 2011/4/2 stan zimine > Hi. > Googled but did not found the answer for the following little issue. > > how to force

Re: [R] another question on shapefiles and geom_point in ggplot2

2011-04-03 Thread Felipe Carrillo
Manuel: I changed your variable names from x to 'long' and y to 'lat' on the riqueza_out.csv file. The code below should do what you want. Also, since the legend title is kind of long, I broke it down into three lines so you can see more plot area. I am cc'ing the other groups so more people use

[R] zoo:rollapply by multiple grouping factors

2011-04-03 Thread Mark Novak
# Hi there, # I am trying to apply a function over a moving-window for a large number of multivariate time-series that are grouped in a nested set of factors. I have spent a few days searching for solutions with no luck, so any suggestions are much appreciated. # The data I have are for the

[R] setCoefTemplate

2011-04-03 Thread Sebastián Daza
Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Reje

Re: [R] kernel density plot

2011-04-03 Thread Greg Snow
It would be easier to test and example if we knew what s22, s33, and s44 were, but you can try: density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t plot(t, xlim=c(-3,4), main="Net Income Distribution", axes=F, ylim=c(0,0.00035). xlab="Value in Rupees") density(s33

Re: [R] Help in splitting ists into sub-lists

2011-04-03 Thread andrija djurovic
Hi: here is one solution. Not so elegant but maybe it will give you some ideas: mylist1<-rep(mylist,c(2,2)) a<-matrix(c(index1,!index1,index2,!index2),ncol=4) mylist2<-list() for (i in 1:4) { mylist2[[i]]<-mylist1[[i]][,a[,i]] } mylist2 Andrija On Sun, Apr 3, 2011 at 6:56 PM, Axel Urbiz wrote

Re: [R] kernel density plot

2011-04-03 Thread Muzna Alvi
i am sorry greg, can you explain that with an example? On Mon, Apr 4, 2011 at 12:11 AM, Greg Snow wrote: > It is better to replace your later calls to plot with calls to lines > instead, then you don't need to use par(new=T) which as you see tends to > cause problems. > > -- > Gregory (Greg) L.

Re: [R] kernel density plot

2011-04-03 Thread Greg Snow
It is better to replace your later calls to plot with calls to lines instead, then you don't need to use par(new=T) which as you see tends to cause problems. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Messag

Re: [R] Function for finding NA's

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 1:44 PM, Tyler Rinker wrote: Quick question, I tried to find a function in available packages to find NA's for an entire data set (or single variables) and report the row of missing values (NA's for each column). I searched the typical routes through the blogs and th

Re: [R] converting "call" objects into character

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 1:22 PM, Douglas Bates wrote: On Sun, Apr 3, 2011 at 11:42 AM, David Winsemius > wrote: On Apr 3, 2011, at 12:14 PM, Samuel Le wrote: Dear all, I would like to log the calls to my functions. I am trying to do this using the function match.call(): fTest<-function(x)

Re: [R] :HELP

2011-04-03 Thread Dennis Murphy
Hi: Is this what you're after? x <- matrix(1:10, nrow = 2, byrow = TRUE) rowSums(x[, 1:3]) HTH, Dennis PS: This list expects submissions to be in text, not HTML. See the Posting Guide, please. On Sun, Apr 3, 2011 at 6:10 AM, mustafabinar wrote: > Hello, >   > I want to sum first three terms o

[R] Function for finding NA's

2011-04-03 Thread Tyler Rinker
Quick question, I tried to find a function in available packages to find NA's for an entire data set (or single variables) and report the row of missing values (NA's for each column). I searched the typical routes through the blogs and the help manuals for 15 minutes. Rather than spend any

Re: [R] converting "call" objects into character

2011-04-03 Thread Douglas Bates
On Sun, Apr 3, 2011 at 11:42 AM, David Winsemius wrote: > > On Apr 3, 2011, at 12:14 PM, Samuel Le wrote: > >> Dear all, >> >> >> >> I would like to log the calls to my functions. I am trying to do this >> using the function match.call(): > > fTest<-function(x) > > {  theCall<-match.call() >      

[R] Help in splitting ists into sub-lists

2011-04-03 Thread Axel Urbiz
Dear List, Let's say I have a list whose components are 2 matrices (as exemplified in the "mylist" object below). I'd like to create a list with components being 4 matrices based on an logical index vector. is there a way to simplify what I'm doing to obtain the results in "mylist2"? I'd like some

Re: [R] converting "call" objects into character

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 12:14 PM, Samuel Le wrote: Dear all, I would like to log the calls to my functions. I am trying to do this using the function match.call(): fTest<-function(x) { theCall<-match.call() print(theCall) return(list(x=x, logf = theCall)) } > fTest(x=2)$x [1

[R] converting "call" objects into character

2011-04-03 Thread Samuel Le
Dear all, I would like to log the calls to my functions. I am trying to do this using the function match.call(): fTest<-function(x) { theCall<-match.call() print(theCall) return(x) } > fTest(2) fTest(x = 2) [1] 2 I can see "theCall" printed into the console, bu

Re: [R] :HELP

2011-04-03 Thread John Kane
Your message is severely garbled. However a={1,2,3,4,5 6,7,8,9,10} is not going to work since {} is incorrect. you need c() Anyway try this: a=matrix(1:10,ncol=2) sum(a[1:3,1]) sum(a[1:3,2]) --- On Sun, 4/3/11, mustafabinar wrote: > From: mustafabinar > Subject: [R] :HELP > To: r-help@r-p

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread John Fox
Dear Krishna, > -Original Message- > From: Krishna Kirti Das [mailto:krishnaki...@gmail.com] > Sent: April-03-11 10:36 AM > To: John Fox > Cc: r-help@r-project.org > Subject: Re: [R] Unbalanced Anova: What is the best approach? > > Thank you, John. > > Yes, your answers do help. For me i

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread John Fox
Dear Spencer, > -Original Message- > From: Spencer Graves [mailto:spencer.gra...@prodsyse.com] > Sent: April-03-11 11:07 AM > To: Krishna Kirti Das > Cc: John Fox; r-help@r-project.org > Subject: Re: [R] Unbalanced Anova: What is the best approach? > > Hi, Krishna: > > > > > On 4/3/20

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread Spencer Graves
Hi, Krishna: On 4/3/2011 7:35 AM, Krishna Kirti Das wrote: Thank you, John. Yes, your answers do help. For me it's mainly about getting familiar with the "R" way of doing things. Thus your response also confirms what I suspected, that there is no explicit user-interface (at least one that i

[R] :HELP

2011-04-03 Thread mustafabinar
Hello,   I want to sum first three terms of each column of matrix. But I don't calculate with "apply" function.   skwkrt<-function(N=1,mu=0,sigma=1,n=100, nboot=1000,alpha=0.05){ x<-rnorm(N,mu,sigma)#population samplex<-matrix(sample(x,n*nboot,replace=T),nrow=nboot) #... }   is that: suppose a

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread Krishna Kirti Das
On Sun, Apr 3, 2011 at 3:10 AM, peter dalgaard wrote: > > On Apr 3, 2011, at 09:24 , Krishna Kirti Das wrote: > > > I have a three-way unbalanced ANOVA that I need to calculate (fixed > effects > > plus interactions, no random effects). But word has it that aov() is good > > only for balanced des

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread Krishna Kirti Das
Thank you, John. Yes, your answers do help. For me it's mainly about getting familiar with the "R" way of doing things. Thus your response also confirms what I suspected, that there is no explicit user-interface (at least one that is widely used) in terms of functions/packages that represents an

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread John Fox
Dear Krishna, Although it's difficult to explain briefly, I'd argue that balanced and unbalanced ANOVA are not fundamentally different, in that the focus should be on the hypotheses that are tested, and these are naturally expressed as functions of cell means and marginal means. For example, in a

Re: [R] Inverse noncentral Beta

2011-04-03 Thread Ted Harding
On 03-Apr-11 13:03:30, Eduardo M. A. M.Mendes wrote: > Hello > I could not find whether there is any R-function that implements > the inverse of a noncentral Beta. Could someone out there tell > me where I can find it? Or how to implement it? > > Many thanks > Ed Have a look at the 'ncp' paramat

Re: [R] recommendation on r scripting tutorial?

2011-04-03 Thread Duncan Murdoch
On 11-04-03 9:01 AM, lcn wrote: The documents accompanying the distribution can be a good start. And a suggestion for searching help on R over Google, use "r-help" as a basic keyword, coz a single letter of "r" hardly helps you find the desired topics. When I google for "R tutorial" or "r tuto

Re: [R] R gui on windows how to force to always show the last line of output

2011-04-03 Thread lcn
CTRL+L helps. 2011/4/2 stan zimine > Hi. > Googled but did not found the answer for the following little issue. > > how to force R gui on windows (maybe a specific setting) to always > show the last line of output in the window console. > > > My program in R makes measurements every 5 mins in

[R] Inverse noncentral Beta

2011-04-03 Thread Eduardo M. A. M.Mendes
Hello I could not find whether there is any R-function that implements the inverse of a noncentral Beta. Could someone out there tell me where I can find it? Or how to implement it? Many thanks Ed [[alternative HTML version deleted]] ___

Re: [R] recommendation on r scripting tutorial?

2011-04-03 Thread lcn
The documents accompanying the distribution can be a good start. And a suggestion for searching help on R over Google, use "r-help" as a basic keyword, coz a single letter of "r" hardly helps you find the desired topics. 2011/4/2 Wensui Liu > Good morning, dear listers > > I am wondering if you

Re: [R] Error in "color2D.matplot" : "Error in plot.new() : figure margins too large"

2011-04-03 Thread John Kane
Totally guessing but did you make an earlier call to par() to adjust margins? Otherwise you may want to supply the matrix here for other people to experiment with. See ?dput as probably the best way to do this. --- On Sun, 4/3/11, shahab wrote: > From: shahab > Subject: [R] Error in "color

Re: [R] I think I just broke R

2011-04-03 Thread Alexander Engelhardt
Am 03.04.2011 14:10, schrieb Duncan Murdoch: That's not weird: you just created an R variable named LANGUAGE, not an environment variable. Duncan Murdoch Silly me. It works now: alexx@derp:~$ cat ~/.Renviron LANGUAGE="EN" Thanks :) __ R-help@r-pro

Re: [R] I think I just broke R

2011-04-03 Thread Duncan Murdoch
On 11-04-03 7:50 AM, Alexander Engelhardt wrote: Am 03.04.2011 13:12, schrieb Prof Brian Ripley: On Sun, 3 Apr 2011, Uwe Ligges wrote: Use "LANGUAGE" rather than "LANG" as the environment variable. Also, set it outside your R session, e.g. in your .Renviron file. You are supposed to be able

[R] Error in "color2D.matplot" : "Error in plot.new() : figure margins too large"

2011-04-03 Thread shahab
Hi, I am using color2D.matplot (...) function of "plotrix" package. I used a matrix of size around 20*20 However, apparently it failed to visualize the matrix and gave the following exception, which I don't have any idea about possible source of this error. "Error in plot.new() : figure margins

Re: [R] help

2011-04-03 Thread Mike Marchywka
> Date: Sun, 3 Apr 2011 01:35:16 +0530 > From: nandan.a...@gmail.com > To: padmanabhan.vija...@gmail.com > CC: r-help@r-project.org > Subject: Re: [R] help > > One way that u might have thought of is to create plot in PDF in R and the > use pdftools. > Additionally one can also think of running

Re: [R] kernel density plot

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 7:55 AM, David Winsemius wrote: On Apr 3, 2011, at 6:56 AM, Muzna Alvi wrote: I am using the following commands for plotting kernel density for three kinds of crops density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t plot(t, xlim=c(-3,4)

Re: [R] kernel density plot

2011-04-03 Thread David Winsemius
On Apr 3, 2011, at 6:56 AM, Muzna Alvi wrote: I am using the following commands for plotting kernel density for three kinds of crops density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t plot(t, xlim=c(-3,4), main="Net Income Distribution", axes=F, ylim=c(0,0.0

Re: [R] I think I just broke R

2011-04-03 Thread Alexander Engelhardt
Am 03.04.2011 13:12, schrieb Prof Brian Ripley: On Sun, 3 Apr 2011, Uwe Ligges wrote: Use "LANGUAGE" rather than "LANG" as the environment variable. Also, set it outside your R session, e.g. in your .Renviron file. You are supposed to be able to change this during an R session, but if you rel

Re: [R] I think I just broke R

2011-04-03 Thread Prof Brian Ripley
On Sun, 3 Apr 2011, Uwe Ligges wrote: On 03.04.2011 09:30, Alexander Engelhardt wrote: Am 03.04.2011 03:51, schrieb Daniel Malter: Check whether x, y, or glm have been redefined. If not, restart R. I wouldn't call my function 'glm'. However, I did call one 'binomial'. That was my mistake.

[R] Plotting data on a US County Map

2011-04-03 Thread mathijsdevaan
Hi, I have a data frame listing US counties and a quantity ("number") per county and I have a shapefile of the US with county ID's. I would like to plot the "number" variable on a map (in the shapefile) using a color range per county (e.g. white = min(number) = 2, black = max(number) = 15). Can a

[R] kernel density plot

2011-04-03 Thread Muzna Alvi
I am using the following commands for plotting kernel density for three kinds of crops density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t plot(t, xlim=c(-3,4), main="Net Income Distribution", axes=F, ylim=c(0,0.00035). xlab="Value in Rupees") par(new=T) density(s3

Re: [R] I think I just broke R

2011-04-03 Thread Uwe Ligges
On 03.04.2011 09:30, Alexander Engelhardt wrote: Am 03.04.2011 03:51, schrieb Daniel Malter: Check whether x, y, or glm have been redefined. If not, restart R. I wouldn't call my function 'glm'. However, I did call one 'binomial'. That was my mistake. Thanks :) A few weeks ago I asked how t

Re: [R] Discretizing data rows into regular intervals

2011-04-03 Thread Gabor Grothendieck
On Sat, Apr 2, 2011 at 9:31 PM, Linh Tran wrote: > Hi guys, > > I'd like to thank you ahead of time for any help that you can offer me. > I'm kind of stuck trying to do this. > > I have a data frame with dates and values (note: only two columns shown): > > head(test) >        date     value      

Re: [R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread peter dalgaard
On Apr 3, 2011, at 09:24 , Krishna Kirti Das wrote: > I have a three-way unbalanced ANOVA that I need to calculate (fixed effects > plus interactions, no random effects). But word has it that aov() is good > only for balanced designs. I have seen a number of different recommendations > for workin

[R] Unbalanced Anova: What is the best approach?

2011-04-03 Thread Krishna Kirti Das
I have a three-way unbalanced ANOVA that I need to calculate (fixed effects plus interactions, no random effects). But word has it that aov() is good only for balanced designs. I have seen a number of different recommendations for working with unbalanced designs, but they seem to differ widely (car

Re: [R] I think I just broke R

2011-04-03 Thread Alexander Engelhardt
Am 03.04.2011 03:51, schrieb Daniel Malter: Check whether x, y, or glm have been redefined. If not, restart R. I wouldn't call my function 'glm'. However, I did call one 'binomial'. That was my mistake. Thanks :) A few weeks ago I asked how to set my error messages to english, and Richard H