Re: [R] manipulating data in heatmap

2012-05-11 Thread fjucks
but, df1 is the row and column? As the call to dput defni them? eh my script this, but where I put these lines to order? library(gplots) arq <-read.table("table") matrix_l <-data.matrix(arq) pdf("heatmap.pdf", height = 10 , width=10) #paleta de 10 cores - sentido branco -> preto my.colors <- c

Re: [R] Multiple regression Categorical data

2012-05-11 Thread David Winsemius
On May 11, 2012, at 9:13 AM, Rosario Garcia Gil wrote: Hello It is possible to set up an lm() model where none of the categories of the categorical independent variable need to be used as references, I mean use the total mean instead. Yes. ?contrasts -- David Winsemius, MD Heritage Lab

[R] useR! 2012: Call for Late-Breaking Posters; REGULAR REGISTRATION ENDS 12May

2012-05-11 Thread Frank Harrell
*** Call for Late-breaking Posters *** Abstracts may be submitted for posters presenting recent developments and late-breaking applications of R, on topics as indicated in the earlier call for abstracts: http://biostat.mc.vanderbilt.edu/UseR-2012#Call_for_Abstracts_and_Tutorial Late-breaking

[R] Why can we combine design matrix and data-frame in R?

2012-05-11 Thread Michael
Hi all, Could you please help me? I am trying to understand why this line works: lm1x = lm(y~X-1, tmp) Here it seems that I was combining the design matrix and the data frame... And X below is not a single column, in fact, it's a bunch of columns in matrix form... I don't understand why this

[R] NLS sensitivity to start= values or poles in data range

2012-05-11 Thread Christopher Battles
Greetings R-help! I'm fairly new to R and am trying to expand my knowledge beyond using R for simple summary statistics and basic tests. To that end I am attempting to write an interactive R-script that will perform a general rational function fit to a given dataset based on the example given

[R] manipulating data in heatmap

2012-05-11 Thread fjucks
manipulating data in heatmap, have the following data: RF00505 RF00232 RF00493 RF00231 RF00505100.00 20.00 0.00 0.00 RF00232100.00 00.00100.00 100.00 RF00493100.00 60.00100.00 33.33 RF0023

Re: [R] summary for weekdays()

2012-05-11 Thread Rui Barradas
Hello Try table(mydata$day) Works with any type of character input. Hope this helps, Rui Barradas Kai Mx wrote > > Hi all, > probably really simple: I seem to be lacking some understanding for the > character class: I have a bunch of dates in a dataframe and I want to add > a > string variab

Re: [R] Matrix invert

2012-05-11 Thread R. Michael Weylandt
a[NROW(a):1, ] Michael On May 11, 2012, at 3:00 PM, Trying To learn again wrote: > I all, > > I have a matrix like this > > a= > > 1 4 > 2 7 > 3 6 > > > I want to create a new matrix > > b= > > 3 6 > 2 7 > 1 4 > > Anyone knows if there is a "reverse" function? > I can do it with loop

Re: [R] Result of clustering on plot

2012-05-11 Thread greyofthedawn
Yes, now it works! Thank you very much for your help! Jane -- View this message in context: http://r.789695.n4.nabble.com/Result-of-clustering-on-plot-tp4625043p4627431.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.

[R] Matrix invert

2012-05-11 Thread Trying To learn again
I all, I have a matrix like this a= 1 4 2 7 3 6 I want to create a new matrix b= 3 6 2 7 1 4 Anyone knows if there is a "reverse" function? I can do it with loops if no exits. [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] Error t value matrix

2012-05-11 Thread Trying To learn again
I have it, I have used a simply programmation using loops I have to run the multiple lm with at least 3 variables (because it needs at least 2 degrees of freedom). Many Thanks 2012/5/11 Trying To learn again > Hi all, > > I have I got it. I see the way on my mind. I will try this week en

[R] Calculating all possible ratios

2012-05-11 Thread genome1976
I have a data matrix with genes as columns and samples as rows. I want to create all possible gene ratios.Is there an elegant and fast way to do it in R and write it to a dataframe? Thanks for any help. Som. -- View this message in context: http://r.789695.n4.nabble.com/Calculating-all-possible

[R] Random forests prediction

2012-05-11 Thread matt
Hi all, I have a strange problem when applying RF in R. I have a set of variables with which I obtain an AUC of 0.67. I do have a second set of variables that have an AUC of 0.57. When I merge the first and second set of variables, the AUC becomes 0.64. I would expect the prediction to becom

Re: [R] How to re-order clusters of hclust output?

2012-05-11 Thread Emmanuel Levy
Hi Michael, Thanks for the info! I think the key function I was missing is: order.dendrogram and not REorder.dendrogram. It returns me the new order, so I think I should get going with that :) Emmanuel On 11/05/2012, R. Michael Weylandt wrote: > I don't have a general answer to your questio

[R] Difference of AIC computation between R (>2.12) and Splus (7.0.6) during stepwise GAM analysis

2012-05-11 Thread Sebastien Bihorel
Dear R Users, I was wondering if some members of the list could shed some light on the difference in AIC computation existing between R (>2.12; gam package) and Splus (7.0.6). Because I am not a statistician by training, I would like to apologize in advance if I use wrong terms or dot not desc

Re: [R] How to re-order clusters of hclust output?

2012-05-11 Thread R. Michael Weylandt
I don't have a general answer to your question, but 1L and 2L are just the integers 1 and 2 (the L makes them integers instead of doubles which is useful for some things) Michael On May 11, 2012, at 2:15 PM, Emmanuel Levy wrote: > Hello, > > The heatmap function conveniently has a "reorder.d

Re: [R] Matching - finding and listing data that is one object but not in another

2012-05-11 Thread R. Michael Weylandt
? setdiff Michael On May 11, 2012, at 8:50 AM, James Holland wrote: > What is the best way to find out what elements/numbers that are in one > object are not in another. > > > I came up with this method, but I'm wondering if there is a more efficient > way (and one that doesn't seem so "clunk

[R] How to re-order clusters of hclust output?

2012-05-11 Thread Emmanuel Levy
Hello, The heatmap function conveniently has a "reorder.dendrogram" function so that clusters follow a certain logic. It seems that the hclust function doesn't have such feature. I can use the "reorder" function on the dendrogram obtained from hclust, but this does not modify the hclust object it

[R] contrasts with an imbalance in a factor

2012-05-11 Thread Connolly, Colm
Hi everybody, I have an experiment examining risky choice behavior where two groups of subjects were unevenly divided across two different MRI scanners while they performed a task. Each subject's data was recorded once and only once on a particular scanner. The table describing the distribution

Re: [R] Strange "Error: subscript out of bounds"

2012-05-11 Thread Ethan Brown
Hi, this looks like a typo to me. The name of the argument to your function is 'pre.mat', but you're trying to print an object called 'pred.mat' (with an extra 'd') that never appears before. It's easier to help when you give a reproducible example that we can execute on our own computers, as reco

Re: [R] Replicate weights in R?

2012-05-11 Thread Thomas Lumley
Yes. You want the survey package, Use the svrepdesign() function to create a survey object containing the replicate weights, and the rest of the functions in the package to do analyses -thomas On Fri, May 11, 2012 at 9:24 AM, Samantha Tornello wrote: > Hello everyone! I am analyzing a publi

[R] Interpreting summary of lme

2012-05-11 Thread a.lesp
Dear mixed-modelers, I have built a mixed model and I'm having serious trouble with interpreting the output. I want to test differences in the coefficient of variation (CV) of light across 3 tree crown exposures (Depth). I have measured direct and diffuse radiation (RF) and I want to test for dif

Re: [R] Error t value matrix

2012-05-11 Thread Trying To learn again
Hi all, I have I got it. I see the way on my mind. I will try this week end and if I get it I will send. Thanks. 2012/5/10 Trying To learn again > Hi all, > > I want to make the following: > > I want to run a linear regression on each column of a matrix "estima" on > the correspondent column

Re: [R] Multiple regression Categorical data

2012-05-11 Thread Greg Snow
?contr.sum On Fri, May 11, 2012 at 8:13 AM, Rosario Garcia Gil wrote: > Hello > > It is possible to set up an lm() model where none of the categories of the > categorical independent variable need to be used as references, I mean use > the total mean instead. > > Thanks > /R > _

Re: [R] How to specify multiple regular expressions for pattern argument

2012-05-11 Thread Berend Hasselman
On 11-05-2012, at 20:18, Jun Shen wrote: > Dear list, > > I want to show all the objects starting with "d" and ending with a digit. > How do I specify these conditions in the pattern argument > > I can do one condition but not two > > ls(pattern='^d') > ls(pattern='[[:digit:]]$') > > are work

Re: [R] How to specify multiple regular expressions for pattern argument

2012-05-11 Thread Tyler Rinker
I'm curious about the 'real' answer myself but this would work: ls(pattern='^d')[ls(pattern='^d') %in% ls(pattern='[[:digit:]]$')] Cheers,Tyler Rinker > Date: Fri, 11 May 2012 13:18:59 -0500 > From: jun.shen...@gmail.com > To: r-h...@stat.math.ethz.c

Re: [R] identify() doesn't return "true" numbers

2012-05-11 Thread Jean V Adams
Daniel, There are a few ways to deal with this. You could sort your data by y before you apply these functions. Then the point labelled 100 will be the 100th row in the data frame. df <- data.frame(x=1:100, y=(1:100)+rnorm(100, sd=100) df2 <- df[order(df$y), ] windows() qqInteractive(df$x, df$y)

Re: [R] How to specify multiple regular expressions for pattern argument

2012-05-11 Thread Jeff Newmiller
ls(pattern='^d.*[[:digit:]]$') --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. P

Re: [R] text(): combine expression and line break

2012-05-11 Thread Bert Gunter
Try this: plot(1,type="n", xaxt='n', yaxt='n', ann=FALSE) text(1,1,labels=expression(atop(sigma,"(log scale, m)")),cex = 2) You can add extra blank characters to align the text like this: plot(1,type="n", xaxt='n', yaxt='n', ann=FALSE) text(1,1,labels=expression(atop(sigma,phantom("W")~"(l

Re: [R] text(): combine expression and line break

2012-05-11 Thread jim holtman
try this: plot(1,type="n", xaxt='n', yaxt='n', ann=FALSE) text(1,1,labels=expression(atop(sigma, "log scale,m")),cex = 2) On Fri, May 11, 2012 at 1:40 PM, Johannes Radinger wrote: > Hi, > > I would like to plot some extra text in my plot. > This should be a two line text including a special c

[R] How to specify multiple regular expressions for pattern argument

2012-05-11 Thread Jun Shen
Dear list, I want to show all the objects starting with "d" and ending with a digit. How do I specify these conditions in the pattern argument I can do one condition but not two ls(pattern='^d') ls(pattern='[[:digit:]]$') are working. But, ls(pattern='^d'&'[[:digit:]]$') is not working. Appr

Re: [R] Result of clustering on plot

2012-05-11 Thread Jean V Adams
You don't give much information about x, but assuming that the "names" you want to add to the plot are the row names. Try this text(x, col=cl$cluster, rownames(x)) Jean greyofthedawn wrote on 05/10/2012 06:46:05 PM: > Hello! I 'm new to R and need your help in one question... > I did cluste

[R] Replicate weights in R?

2012-05-11 Thread Samantha Tornello
Hello everyone! I am analyzing a public use dataset which has both basic weights and replicated weights. I have found that only Stata seems to work with replicate weights but I would rather use R. Does anyone know if R can use replicate weights? If this is possible, could someone share some links

Re: [R] as.function parameters

2012-05-11 Thread Rui Barradas
Hello, I'm not completely sure if I understand, but maybe using environments. If each function copy in the list has a different environment, you can set the parameters values in that environment. Something like f <- as.function(alist(y=,x=tmp,y+x)) env <- list( e1=new.env(), e2=

Re: [R] Interweaving of two datasets

2012-05-11 Thread lunarossa
Ok, merging in this way the total number of rows of the new dataset is the number of rows of the dataset "request". Is it possible to obtain a new dataset with the number of rows of the dataset "phone"? -- View this message in context: http://r.789695.n4.nabble.com/Interweaving-of-two-datasets-t

Re: [R] Fisher Test in R

2012-05-11 Thread arun
Hi Aayush, You are getting different results for fisher.test with > mat2<-matrix(c(5,10,60,30), > nrow=2,dimnames=list(c("Dieting","Non-Dieting"),c("Men","Women"))) is because the first test used one-tailed (alternative="greater") while the default without the alternative option is two-tail

[R] Multiple regression Categorical data

2012-05-11 Thread Rosario Garcia Gil
Hello It is possible to set up an lm() model where none of the categories of the categorical independent variable need to be used as references, I mean use the total mean instead. Thanks /R __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

[R] Matching - finding and listing data that is one object but not in another

2012-05-11 Thread James Holland
What is the best way to find out what elements/numbers that are in one object are not in another. I came up with this method, but I'm wondering if there is a more efficient way (and one that doesn't seem so "clunky"). #Example id <- c(1,2,3,4,5,6,7,9,10) example.1 <- data.frame(id) #Second ob

[R] Strange "Error: subscript out of bounds"

2012-05-11 Thread Petri Lankoski
Dear all, I am trying to write a function for visualizing ordinal model results. The function works fine with some values, but then I get "Error: subscript out of bounds" even there the index should be pointing a legal item. Code is below as well as the example of failure: plotProb <- func

Re: [R] Remove a number from a vector

2012-05-11 Thread Zhou Fang
Better yet, remove the which altogether, and it'll run a slight bit faster and maybe look a little neater. x <- x[x!="bobo"] -- View this message in context: http://r.789695.n4.nabble.com/Remove-a-number-from-a-vector-tp851865p4626413.html Sent from the R help mailing list archive at Nabble.com.

[R] identify() doesn't return "true" numbers

2012-05-11 Thread greatest.possible.newbie
Dear R community. I am using the identify() function to identify outliers in my dataset. This is the code I am using: # Function to allow identifying points in the QQ plot (by mouseclicking) qqInteractive <- function(..., IDENT

Re: [R] as.function parameters

2012-05-11 Thread jackl
Thanks for the fast answer.. The problem with your "free" definition of f is, that it does not really fit my task. I have to build a tree where each node contains the same function just with adjusted parameters. Setting the parameters in the function call is not really an option. Is there anothe

[R] summary for weekdays()

2012-05-11 Thread Kai Mx
Hi all, probably really simple: I seem to be lacking some understanding for the character class: I have a bunch of dates in a dataframe and I want to add a string variable with the weekday for each date. If I use something like mydata$day <- weekdays (mydata$date), I can create subsets for each wee

[R] text(): combine expression and line break

2012-05-11 Thread Johannes Radinger
Hi, I would like to plot some extra text in my plot. This should be a two line text including a special character (sigma). I tried so far a to use expression in combination with paste and "\n"... but I can't get the line break... Here what I've done so far: plot(1,type="n", xaxt='n', yaxt='n', a

[R] bigmemory

2012-05-11 Thread ya
Hi Jay, Thanks for the reply:) Could you show me a link of the c++ examples? Thank you very much. ya ya From: Jay Emerson Date: 2012-05-11 20:12 To: xinxi813 CC: r-help Subject: Re: bigmemory R internally uses 32-bit integers for indexing (though this may change). For this and other reaso

[R] Using xpathapply or getnodeset to get text between two distinct tags

2012-05-11 Thread Simon Kiss
Hello: The following code extracts the links to the daily transcripts of Canada's House Of Commons. 'links' is a matrix of URLs (ncol=1), each of which points to one day's transcripts. If you inspect the code for scrape(links[1]), you will find that periodically there appears an italicitze t

Re: [R] bigmemory

2012-05-11 Thread Jay Emerson
R internally uses 32-bit integers for indexing (though this may change). For this and other reasons these external objects with specialized purposes (larger-than-RAM, shared memory) simply can't behave exactly as R objects. Best case, some R functions will work. Others would simply break. Others

Re: [R] mgcv: inclusion of random intercept in model - based on p-value of smooth or anova?

2012-05-11 Thread Martijn Wieling
Dear Simon, Thanks for your concise reply, this is very helpful. With respect to my second question, however, I was not entirely clear - or perhaps I'm misunderstanding your answer. What I meant is: suppose I have a model with a random effect s(X, bs="re"). Now I want to test if a certain (fixed-

[R] bigmemory

2012-05-11 Thread ya
Hi Jay, I have a question about your reply. You mentioned that "the more serious problem is that you can't expect to run just any R function on a big.matrix (or on an ff object, if you check out ff for some nice features). " I am confused why the packages could not communicate with each ot

[R] Fisher Test in R

2012-05-11 Thread Aayush Raman
Suppose we have the following data set: MenWomen Dieting 10 30 Non-dieting 5 60 If I run the Fisher exact test in R then what does alternative = greater (or less) imply? For example: mat = matrix(c(10,5,30,60), 2,2) fisher.test(mat,alternative ="greate

Re: [R] mgcv: inclusion of random intercept in model - based on p-value of smooth or anova?

2012-05-11 Thread Simon Wood
Dear Martijn, Thanks for the off line code and data: very helpful. The answer to this is something of a 'can of worms'. Starting with the p-value inconsistency. The problem here really is that neither test is well justified in the case of s(...,"re") terms (and not having realised the extent

Re: [R] domain/number line/range reduction problem

2012-05-11 Thread Ben quant
Clarification/correction: Ex 5 isn't consistent with the other examples To be consistent with the other examples the resulting ranges would be something like: xa_rng = c(-100,-100) xb_rng = c(100,100) or just... xa_rng = -100 xb_rng = 100 However, my original Ex 5 would be a good solution if

[R] domain/number line/range reduction problem

2012-05-11 Thread Ben quant
Hello, Currently I'm only coming up with brute force solutions to this issue. Wondering if anyone knows of a better way to do this. The problem: I have endpoints of one x range (x_rng) and an unknown number of s ranges (s[#]_rng) also defined by endpoints. What I want are the parts of the x range

Re: [R] Remove a number from a vector

2012-05-11 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf > Of Berend Hasselman > Sent: Friday, May 11, 2012 4:21 AM > To: dkkhireche > Cc: r-help@r-project.org > Subject: Re: [R] Remove a number from a vector > > > On 11-05-2012, at 12:4

Re: [R] ANOVA question

2012-05-11 Thread Robert Latest
Hello Thierry, thanks for your answer! There is one thing, however, that I don't understand. The values labeled "B" in my data are generated with 1/20th the variance of the others, yet the standard error and confidence intervals are the same for all levels of the factor. How come? > summary(mod_

Re: [R] Output of system() not shown in if

2012-05-11 Thread syrvn
Ah I managed to do it by simple adding a print() statement around the system call... easy :) -- View this message in context: http://r.789695.n4.nabble.com/Output-of-system-not-shown-in-if-tp4626383p4626461.html Sent from the R help mailing list archive at Nabble.com. ___

[R] Output of system() not shown in if

2012-05-11 Thread syrvn
Hello, I am using system(...) to invoke some unix commands. Everything is working fine and the output of the programme is returned to the R console. However, as soon as I put the system(...) call into an if-clause the output is no longer shown in the R console. Is there any workaround for that pro

Re: [R] plot: many bars with wide bars

2012-05-11 Thread Jim Lemon
On 05/11/2012 09:33 PM, Marion Wenty wrote: Hello, I have got another problem. I wrote a function which is able to create plots with different numbers of bars. In order to give them a similar look I determined the width of the bars using the parameters width=0.59 AND ylim=c(0,10)) (see fun

Re: [R] Possible artifacts in cross-correlation function ("ccf")?

2012-05-11 Thread Duncan Murdoch
On 11/05/2012 7:06 AM, Tim Dorscheidt wrote: Dear R-users, I have been using R and its core-packages with great satisfaction now for many years, and have recently started using the "ccf" function (part of the "stats" package version 2.16.0), about which I have a question. The "ccf"-algorithm

[R] set specific contrasts using lapply

2012-05-11 Thread Frank Paetzold
I have the following data set > data A B X1 X2 Y 1 A1 B1 1.1 2.9 1.2 2 A1 B2 1.0 3.2 2.3 3 A2 B1 1.0 3.3 1.6 4 A2 B2 0.5 2.6 3.1 > sapply(data, class) A BX1X2 Y "factor" "factor" "numeric" "numeric" "numeric" I'd like to set a specific type of

[R] Odp: arguments must have same length

2012-05-11 Thread Petr PIKAL
Hi > > > score<-read.csv("http://users.stat.umn.edu/~chen2285/hw/ACT.csv";) > score<-read.csv("http://users.stat.umn.edu/~chen2285/hw/ACT.csv";) Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : unable to connect to 'users.stat.umn.edu'

[R] Possible artifacts in cross-correlation function ("ccf")?

2012-05-11 Thread Tim Dorscheidt
Dear R-users, I have been using R and its core-packages with great satisfaction now for many years, and have recently started using the "ccf" function (part of the "stats" package version 2.16.0), about which I have a question. The "ccf"-algorithm for calculating the cross-correlation between t

Re: [R] bigmemory

2012-05-11 Thread Jay Emerson
To answer your first question about read.big.matrix(), we don't know what your acc3.dat file is, but it doesn't appear to have been detected as a standard file (like a CSV file) or -- perhaps -- doesn't even exist (or doesn't exist in your current directory)? Next: > In addition, I am planning to

[R] plot: many bars with wide bars

2012-05-11 Thread Marion Wenty
Hello, I have got another problem. I wrote a function which is able to create plots with different numbers of bars. In order to give them a similar look I determined the width of the bars using the parameters width=0.59 AND ylim=c(0,10)) (see function below) Now, this only works if I have j

[R] Cook's distance for lme?

2012-05-11 Thread Beatriz De Francisco
Hi Is there any function that I can use to calculate Cook's distance for an lme? Beatriz de Francisco Mora PhD Student The Scottish Association for Marine Science Scottish Marine Institute Oban PA37 1QA Tel: 06131 559000 (switchboard) Fax: 01631559001 E. beatriz.defranci...@sams.ac.uk

Re: [R] Remove a number from a vector

2012-05-11 Thread Berend Hasselman
On 11-05-2012, at 12:45, dkkhireche wrote: > Example : > *x=c("abba","bobo","cocoa")* > In order to remove "bobo" from *x* you need to do the following: > *x=x[which(x=="bobo")]* > No. That does the opposite of what you describe. You do x <- x[-which(x=="bobo")] or better x <- x[which(x!="bo

Re: [R] converting raster image

2012-05-11 Thread Mintewab Bezabih
Thanks Mike and Micheal for suggestions. It now worked! Many thanks Mintewab Från: Michael Sumner [mdsum...@gmail.com] Skickat: den 11 maj 2012 01:20 Till: Mintewab Bezabih Kopia: r-help@r-project.org Ämne: Re: [R] converting raster image You should provid

Re: [R] Outcome~predictor model evaluation, repeated measurements

2012-05-11 Thread ONKELINX, Thierry
Dear nameless, A mixed model seems reasonable for your kind of data. lme() from nlme or lmer() from lme4 are good starting points. Please note that there is R-sig-mixed-models: a R mailing list dedicated to mixed models. Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoe

Re: [R] barplot: legend: two rows

2012-05-11 Thread Marion Wenty
Hello, thank you very much for your answers. Using the command ncol worked in my special case, if I have got a legend with 4 elements, but with 5 it doesn't work. I am using two functions which I had created and checked what might be the difference, so that the second one doesn't work, but I could

Re: [R] Remove a number from a vector

2012-05-11 Thread dkkhireche
Example : *x=c("abba","bobo","cocoa")* In order to remove "bobo" from *x* you need to do the following: *x=x[which(x=="bobo")]* In case you already know the index of the element you want to get rid of, it is even easier *x=x[-elem.index]* Good luck -- View this message in context: http://r.7

[R] Outcome~predictor model evaluation, repeated measurements

2012-05-11 Thread rad mac
Dear all, I have simple question regarding how to fit a model (i.e. linear) to the data. Say I have 10 subjects with different phenotypes (dependent var Y, identical for a particular subject) and one predictor variable measured 3 times for each subject (X). By other words: Y Subj X 1 1 1.2 1 1 1.

[R] Contour lines within a spplot

2012-05-11 Thread Nicola Cortesi
Hi all, As you know, spplot function of sp package is very powerful and even allows to superimpose contour lines with the option 'contour=TRUE'. The color of the contour line can be changed with the 'col' option. I'd like to know if it possible to improve on this: 1) changing the color of the

[R] Overlapping area script

2012-05-11 Thread achughes
Dear All I would really appreciate some help with a script which a colleague wrote for me (attached), but I am having problems running (and have not been able to contact my colleague). The script is designed to compare the area of suitable habitat in binary projections of a large number of specie

[R] arguments must have same length

2012-05-11 Thread suntime
> score<-read.csv("http://users.stat.umn.edu/~chen2285/hw/ACT.csv";) > interaction.plot(sex,rep(1,861),score,fun=mean,legend=F,main="profile of > sex") Error: tapply(response, list(x.factor, trace.factor), fun) : arguments must have same length > length(sex) [1] 861 > length(type) [1] 861 > le

Re: [R] fda problems - version correct?

2012-05-11 Thread Spencer Graves
Hi, Jim: A user has reported an error in "fdarm-ch09.R": The way I read the following, the exponentiation function for an fd function currently works only with a spline basis, but the "harmonics" variable in the example in "fdarm-ch09.R" has a fourier basis. I'm pretty confident that

[R] ff question

2012-05-11 Thread ya
Hi all, Thank you very much for reading this email. After trying the bigmemory package, I could not figure it out. So I turn to try ff, hopefully someone here have some experience about it? I have a data set (data.frame) with 16459 rows and 457 columns. I am trying to do a multiple imputation

Re: [R] ANOVA question

2012-05-11 Thread ONKELINX, Thierry
Dear Robert, It is easier to use lm instead of aov if you want coefficients for each group. Note that you can use rnorm vectorised. set.seed(0) N <- 100 # sample size MEAN <- c(10, 20, 30, 40, 50) VAR <- c(20,20,1, 20, 20) LABELS <- factor(c("A", "B", "C", "D", "E")) # create a data frame with

[R] ANOVA question

2012-05-11 Thread Robert Latest
Hello all, I'm very satisfied to say that my grip on both R and statistics is showing the first hints of firmness, on a very greenhorn level. I'm faced with a problem that I intend to analyze using ANOVA, and to test my understanding of a primitive, one-way ANOVA I've written the self-contained p

[R] Error on easy way for JoSAE Package

2012-05-11 Thread Ana-Maria Dobre
Dear Mr/Mrs, I have got the the data sets in the attachement (txt). The fact is that on applying the easy method for JoSAE package on these data sets I got the following error: > result <- eblup.mse.f.wrap(domain.data = d.data, lme.obj = fit.lme1) Error in solve.default(asympt.var.covar) : Lapa

Re: [R] randomForest Tutorial

2012-05-11 Thread Nagagautami Kodali
Hi Jim I am unable to access the webinar pdf on RandomForest Implementation. Can you please send it over? Regards -- Gautami This email message may contain proprietary, private and confidential information. The information transmitted is intended only for the

[R] Additional info: help with SMATR: help with pairwise comparisons using MA regression?

2012-05-11 Thread Ioulia Bespalova
Also, this works (taking out multcomp=TRUE, multcompmethod="adjusted"): com.test=ma(Head.W1~Leg.3.1+Site, type="elevation", data=queens) print(com.test) so for some reason it will do an MA regression on all my data point together, but shows an error when I try to do pairwise comparisons betwe

[R] help with SMATR: help with pairwise comparisons using MA regression?

2012-05-11 Thread Ioulia Bespalova
Hi there, I've been using the SMATR package to do standardized major axis (SMA) regression on allometric data, and I've been able to to pairwise comparisons of slope & elevation between multiple study sites. Now I'm trying to do the same thing using major axis (MA) regression, but I'm getting err

[R] how to make http proxy setting as "automatically detect"?

2012-05-11 Thread sagarnikam123
i am connecting through my mobile connection,which don't have proxy setting i can access internet in browser using "automatically detect setting",i get error when trying below > source("http://bioconductor.org/biocLite.R";) Error in file(file, "r", encoding = encoding) : cannot open the connec