[R] R reports

I don't see much in the way of an ability to write reports in R the way you can with SAS. You basically have to write a program with R in a 3G way unlike SAS with it's 4G proc print and proc report. Are there similar R functions and packages? -- View this message in context: http://r.789695.n

Re: [R] modifying factors in a data frame

Kiaora Koe Peter It worked. I hope I am getting my Maori right. -Abhi On Thu, Aug 19, 2010 at 1:50 AM, Peter Alspach wrote: > > Tena koe Abhi > > The easiest way is probably > > data1$sample <- as.numeric(data1$sample) > > If you want it to remain as a factor then > > data1$sample <- as.factor(

[R] (no subject)

  To   R group Help Desk     I am a user of R software.  I am facing a problem while using "frailtyPenal" command in R.2.11.1. When I use these command, R closes completely without any prior alert message. Can I know what would be the reason? My data size is 7050 records with atleast 25 variable

Re: [R] modifying factors in a data frame

Tena koe Abhi The easiest way is probably data1$sample <- as.numeric(data1$sample) If you want it to remain as a factor then data1$sample <- as.factor(as.numeric(data1$sample)) HTH Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-

[R] modifying factors in a data frame

Hi All I have a data frame in which on the columns is a factor with following levels Levels: M006 M0114 M0379 M0432 M0823 M1012 M1096 M1107 SW393 SW708 I want to change all the M006 to "1", M0114 to "2" and so on. What I am trying is not working. data1$sample[data1$sample =='M006'] <- "1" Warn

Re: [R] 'panel.smooth' error

panel.smooth is a base graphics, not a lattice, function. They are incompatible. ?panel.smooth -- Bert On Wed, Aug 18, 2010 at 2:05 PM, Sebastian P. Luque wrote: > On Wed, 18 Aug 2010 13:56:23 -0700 (PDT), > Phil Spector wrote: > >> Sebastian - panel.loess will be more suitable for use with xy

Re: [R] Displaying Results in Two Columns

You could find several good examples on the first several results here: http://www.google.com/search?sourceid=chrome&ie=UTF-8&q=plot.ci+R Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talga

Re: [R] Rolling window linear regression

Thanks, I will try it. Regards Sid Sent on my BlackBerry® from Vodafone -Original Message- From: Dennis Murphy Date: Wed, 18 Aug 2010 08:46:49 To: Subject: Re: [R] Rolling window linear regression This is called kernel-based regression; the most popular version is loess. Try library(

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

Hi all, I figured out why this was happening. It is because my actual code was: lmer(Y~X + (1|as.factor(labs)),data=DATA) In this case, the as.factor function looks for object 'labs' not object 'DATA$labs.' Scope is something you hear about don't worry about until it bites you on your ass I gue

Re: [R] Different way of aggregating

Hi: On Wed, Aug 18, 2010 at 4:43 PM, skan wrote: > > Hi > > > Usually "aggregate" is used to calculate things such as the sum of all data > on the first day, the sum next day, and so on. > But how can I calculate the mean of the first hour of all days, the mean of > the second hour of all days,

Re: [R] Does R always insist on sending plot output to a file?

Donald, At least for the PDF device (I know you asked about png, but I believe they are similar), the answer no. Ultimately, this device calls the standard C function fopen, and writes its data to the resulting file stream. If you're using GNU Linux, you might trick R into writing to a fifo (a na

Re: [R] ifelse command

Hi Dennis, I see your point about using a different name for the runif() functions, other then the one I used to stimulate the bias coins. I start to get what you and David meant after thinking it through for a while regarding with comparison the biased dice and the uniform distribution, but the

Re: [R] ifelse command

On Wed, Aug 18, 2010 at 7:46 PM, Philip Wong wrote: > > well to be honest, it is a assignment for the Bayesian statistic paper I wish > to take later in the academic year.  But I'm a slow learner, so I'm going to > try out some of the assignments posted in the university forum hoping to get > some

Re: [R] Help with ANOVA

From: Tal Galili To: Stephen Liu Cc: r-help@r-project.org Sent: Thu, August 19, 2010 11:05:38 AM Subject: Re: [R] Help with ANOVA >#first: >install.packages("reshape") >#Then try: >library(reshape) >#and then >?melt >?cast Hi Tal, I got it done as advised. Thanks B.R. Stephen L --

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

On 2010-08-18 18:41, Johan Jackson wrote: Hi all, I figured out why this was happening. It is because my actual code was: lmer(Y~X + (1|as.factor(labs)),data=DATA) In this case, the as.factor function looks for object 'labs' not object 'DATA$labs.' Scope is something you hear about don't worr

Re: [R] Help with ANOVA

On 08/18/2010 10:05 PM, Tal Galili wrote: Hi Stephen, #first: install.packages("reshape") #Then try: library(reshape) #and then ?melt ?cast This is also useful: http://www.jstatsoft.org/v21/i12/paper __ R-help@r-project.org mailing list https://sta

Re: [R] Help with ANOVA

Hi Stephen, #first: install.packages("reshape") #Then try: library(reshape) #and then ?melt ?cast Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.

Re: [R] Help with ANOVA

- Original Message From: Tal Galili To: Sneeketeeke Cc: r-help@r-project.org Sent: Thu, August 19, 2010 10:23:04 AM Subject: Re: [R] Help with ANOVA - snip - >Also, you might need to move your data from wide to long. For which you >might wish to look at >?reshape >(or better yet the

Re: [R] ifelse command

Hi Josh, I think I know where does the 50+ warning is coming from, because I used n in runif(n), similarly I get 50+ warnings if I use runif(1000). Yet if I use runif(1) the warnings() doesn't show. -- View this message in context: http://r.789695.n4.nabble.com/ifelse-command-tp2329538p2330597.

Re: [R] reading a text file, one line at a time

Hi Jim, I was trying to use your template without success. With the toy data below, could you explain how to use this template to change all "b"s to "z"s -- just as an exercise, reading in 3 lines at a time. I need to use this strategy for a larger problem, but I haven't been able to get the basic

Re: [R] ifelse command

well to be honest, it is a assignment for the Bayesian statistic paper I wish to take later in the academic year. But I'm a slow learner, so I'm going to try out some of the assignments posted in the university forum hoping to get some practice in advance. Could you please elaborate more on the

Re: [R] 'panel.smooth' error

Hi: This also works: xyplot(incidence ~ year, melanoma, type = c('p', 'smooth')) See p. 75 of the Lattice book for the allowable types. HTH, Dennis On Wed, Aug 18, 2010 at 1:34 PM, Sebastian P. Luque wrote: > Hi, > > The following call: > > xyplot(incidence ~ year, melanoma, panel=panel.smoot

[R] Displaying of Results

I apologize if this posting shows up again, for some reason I wasn't able to post from a different account. So, here I am (reborn). Could I have some suggestions as to how I can display my results and their respective CIs in an aesthetically pleasing manner? Below is the example code. rm(list

Re: [R] Plotting K-means clustering results on an MDS

Hi trzeszutek, Another approach you might want to look at is the clustergram: http://www.r-statistics.com/2010/06/clustergram-visualization-and-diagnostics-for-cluster-analysis-r-code/ Chee

Re: [R] Help with ANOVA

Hi Sneeketeeke, Did you happen to have a look at the example on ?aov and ?TukeyHSD Also, you might need to move your data from wide to long. For which you might wish to look at ?reshape (or better yet the ?melt ?cast commands from the reshape package) Tal Contact Details:---

[R] Hot deck package

Hi.. Anyone knows about package that i can use to make analysis about hot deck imputation? Thank You __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guid

Re: [R] question about unwanted addresses in contact list

The suggestion that "a virus floating around on this distribution list" is baseless. All source code distributed on the list is embedded in messages in plaintext -- i.e. there is no distribution of binary files as part of the mailing list. How would a virus propagate under these circumstances?

Re: [R] Does R always insist on sending plot output to a file?

Take a look on help page for ?recordPlot On Wed, Aug 18, 2010 at 10:49 PM, Donald Paul Winston < satchwins...@yahoo.com> wrote: > > I need to write the output of a R plot to a Java OutputStream. It looks > like > R insists on sending it's output to a file. Is there anyway to get bytes > directly

[R] GLMM random effects

Hello, I have a couple questions regarding generalized linear mixed models specifically around fitting the random effects terms correctly to account for any pseudo-replication. I am reading through and trying to follow examples from Zuur et al. Mixed Effects Models and Extensions in Ecology with

[R] Does R always insist on sending plot output to a file?

I need to write the output of a R plot to a Java OutputStream. It looks like R insists on sending it's output to a file. Is there anyway to get bytes directly from the output of a plot so I can write it with Java? Writing it to a file is too slow. Is there a parameter in the graphics device funct

[R] Adding column to dataframe

Two questions: 1) Are there any good R guides/sites with information/techniques for dealing with large datasets in R? (Large being ~2 mil rows and ~200 columns) 2) My specific problem with this dataset. I am essentially trying to convert a date and add it to a data frame. I imagine any 'data mani

Re: [R] sample() question

Yes it is! I just came up with a method, which needs 5 lines. Thank you very much! -Original Message- From: jim holtman [mailto:jholt...@gmail.com] Sent: Wednesday, August 18, 2010 8:10 PM To: Chen,Shaofei Cc: R Newbie; r-help@r-project.org Subject: Re: [R] sample() question I think thi

Re: [R] sample() question

I think this is what was previously proposed on the list; does this do it for you: > sample(c(1:3, sample(1:3,2,TRUE)), 5) [1] 1 2 3 1 1 > sample(c(1:3, sample(1:3,2,TRUE)), 5) [1] 3 2 1 2 3 > sample(c(1:3, sample(1:3,2,TRUE)), 5) [1] 1 2 3 2 1 > sample(c(1:3, sample(1:3,2,TRUE)), 5) [1] 1 2 3 2 2

Re: [R] sample() question

Yes, sample size is 5, but the size of set is 3. There is no probability to control, or we can say all elements in the set have equal probability to be picked, but I want everyone has to be sampled at least once. Thanks! -Original Message- From: R Newbie [mailto:help0938...@gmail.com] Sen

Re: [R] sample() question

I believe you mean c(sample(1:3,3), sample(1:3, 2, replace = TRUE) is that right? However, we cannot see the combination like 33312. This is a problem. -Original Message- From: Erik Iverson [mailto:er...@ccbr.umn.edu] Sent: Wednesday, August 18, 2010 7:56 PM To: Chen,Shaofei Cc: r-help@

Re: [R] sample() question

On 08/18/2010 07:51 PM, Chen,Shaofei wrote: Consider an example: There are 3 stores (1,2,3), and 5 customers. Each store must have at least one customer. Thus a possible combination is 11213. On the other hand, a combination 22333 is not what I want I have considered your solution earlier, but i

Re: [R] sample() question

Consider an example: There are 3 stores (1,2,3), and 5 customers. Each store must have at least one customer. Thus a possible combination is 11213. On the other hand, a combination 22333 is not what I want I have considered your solution earlier, but in this case, first three customers have to go

Re: [R] sample() question

Just to clarify some of your language here before some others rip you apart, you mean to say you'd like to take a random sample of size 5 and not 5 random samples. Now, I believe you can control the probability with which each element of your original data set is sampled (using weights), but in

[R] memory problem

Hi, when i run the following code i get this massege: "The instruction at 0x reference memory at 0x###, the memory cannot be "read". and then i have to close R. what is the problem and how can i solve it? thanks in advance Avi my code # frailtypack library(frailtypack) cgd.ag <

[R] Displaying Results in Two Columns

Could I have some suggestions as to how (various ways) I can display my confidence interval results? rm(list = ls()) set.seed(1) func <- function(d,t,beta,lambda,alpha,p.gamma,delta,B){ d <- c(5,1,5,14,3,19,1,1,4,22) t <- c(94.32,15.72,62.88,125.76,5.24,31.44,1.048,1.048,2.096,10.

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

On Aug 18, 2010, at 6:45 PM, Peter Ehlers wrote: On 2010-08-18 11:49, Johan Jackson wrote: No, apologies (good catch David!), I merely copied the script incorrectly. It was lmer(Y~X + (1|labs),data=DATA) in my original script. So my question still stands: is it expected behavior for lmer

Re: [R] sample() question

On 08/18/2010 06:58 PM, Chen,Shaofei wrote: Hello all, I have a question regarding sample() in R. For example, I have a set: set<- c(2,3,5) and I want to draw 5 samples from this set, so replacement is true: m<- sample(set, 5, replace=TRUE) However here comes a problem, for example, I will

Re: [R] How to measure correlations in terms of distance and draw them on a 2-dimmentional plot?

Weijian, Look at multidimensional scaling functions such as cmdscale or isoMDS (in MASS) or factor analysis (e.g., factanal in core R or fa in the psych package) or item cluster analysis (ICLUST, also in psych). What you will need to think about is how many dimensions best represent the 10

[R] sample() question

Hello all, I have a question regarding sample() in R. For example, I have a set: set <- c(2,3,5) and I want to draw 5 samples from this set, so replacement is true: m <- sample(set, 5, replace=TRUE) However here comes a problem, for example, I will have (2,3,3,2,5), but I will also get (3

Re: [R] Different way of aggregating

Skan - I'll answer this question: "how can I calculate the mean of the first hour of all days, the mean of the second hour of all days, and so on." Since you didn't include a reproducible example, I have to make some assumptions. Since you are interested in hours, I'm guessing your dates are

[R] Different way of aggregating

Hi Usually "aggregate" is used to calculate things such as the sum of all data on the first day, the sum next day, and so on. But how can I calculate the mean of the first hour of all days, the mean of the second hour of all days, and so on. ??? That's Most examples: today at 1am + today at 2

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

On 2010-08-18 11:49, Johan Jackson wrote: No, apologies (good catch David!), I merely copied the script incorrectly. It was lmer(Y~X + (1|labs),data=DATA) in my original script. So my question still stands: is it expected behavior for lmer to access the object 'labs' rather than the object 'DAT

Re: [R] How to measure correlations in terms of distance and draw them on a 2-dimmentional plot?

Do a search for "multidimensional scaling" and that should give some ideas (use 1-cor or something like that as the distance measure). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message

Re: [R] How to measure correlations in terms of distance and draw them on a 2-dimmentional plot?

Probably I put it too complicated. I just try to draw 10 points on a plot. Each point stand for one variable. If two variables are highly correlated, the corresponding dots are to be drawn closer. Say, variable X1, X2, X3, are highly correlated among them, so they will be clustered; variables X4,

Re: [R] Linear regression equation and coefficient matrix

Well you already have the correlation matrix, so you can just work from there. One of the profound equations in statistics is that r^2=r^2, meaning that to get r2 (as you call it, assuming you mean the coefficient of determination) just square the elements of the correlation matrix. To get reg

Re: [R] How to measure correlations in terms of distance and draw them on a 2-dimmentional plot?

On Aug 18, 2010, at 5:11 PM, weijian21cn wrote: Hi, Suppose I have n variables, X1, X2, ..., Xn and want to visualize the pairwise correlations on a 2-dimmensional plot. In particular, I plan to mark down n dots on a piece of paper and dot i represents Xi, in a way such that the highly c

Re: [R] Looking for Rsafd Macintosh package

On Aug 18, 2010, at 1:39 PM, hankblalo...@yahoo.com wrote: Hi all, Does anybody know where I can get a mac copy of the Rsafd package? I have tried http://orfe.princeton.edu/~rcarmona/SVbook/svbook.html but the link is down there. Thanks! http://R-Forge.R-project.org/ -- David Winsemiu

Re: [R] replacing values in a vector

What do you plan to do with those vectors? Creating a set of vectors and 'asign'ing them as global variables will likely make future tasks harder than they need to be. In the long run it will probably work better for you to either keep them as a matrix (and just access the row/column you need)

[R] How to measure correlations in terms of distance and draw them on a 2-dimmentional plot?

Hi, Suppose I have n variables, X1, X2, ..., Xn and want to visualize the pairwise correlations on a 2-dimmensional plot. In particular, I plan to mark down n dots on a piece of paper and dot i represents Xi, in a way such that the highly correlated variables are placed closer; that is, shorter

Re: [R] 'panel.smooth' error

On Wed, 18 Aug 2010 13:56:23 -0700 (PDT), Phil Spector wrote: > Sebastian - panel.loess will be more suitable for use with xyplot. To > get both points and the smoothed line, you can use > xyplot(incidence~year,melanoma,panel=function(x,y,...){ > panel.xyplot(x,y,...) panel.loess(x,y,...)}) T

Re: [R] 'panel.smooth' error

Sebastian - panel.loess will be more suitable for use with xyplot. To get both points and the smoothed line, you can use xyplot(incidence~year,melanoma,panel=function(x,y,...){ panel.xyplot(x,y,...) panel.loess(x,y,...)})

Re: [R] Looking for Rsafd Macintosh package

Take a look on R-Forge: http://r-forge.r-project.org/R/?group_id=592 On Wed, Aug 18, 2010 at 2:39 PM, hankblalo...@yahoo.com < hankblalo...@yahoo.com> wrote: > Hi all, > Does anybody know where I can get a mac copy of the Rsafd package? I have > tried > http://orfe.princeton.edu/~rcarmona/SV

Re: [R] Looking for Rsafd Macintosh package

On 8/18/2010 10:39 AM, hankblalo...@yahoo.com wrote: Hi all, Does anybody know where I can get a mac copy of the Rsafd package? I have tried http://orfe.princeton.edu/~rcarmona/SVbook/svbook.html but the link is down there. Thanks! A quick search on RSeek.org turns up http://r-forge.r-project

Re: [R] AIC in MuMIn

Yes, I tried the example in the ?dredge and agreed that something else caused the mistake. Aside from the cause which takes time to clarify (16 explanatory variables in the model), I would like to ask another question. Please kindly advise if it is possible to show the singular model with only on

[R] 'panel.smooth' error

Hi, The following call: xyplot(incidence ~ year, melanoma, panel=panel.smooth) produces a blank plot region with an error message: Error using packet 1 plot.new has not been called yet > sessionInfo() R version 2.11.1 (2010-05-31) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

No, apologies (good catch David!), I merely copied the script incorrectly. It was lmer(Y~X + (1|labs),data=DATA) in my original script. So my question still stands: is it expected behavior for lmer to access the object 'labs' rather than the object 'DATA$labs' when using the data= argument? JJ

[R] Looking for Rsafd Macintosh package

Hi all, Does anybody know where I can get a mac copy of the Rsafd package?  I have tried http://orfe.princeton.edu/~rcarmona/SVbook/svbook.html but the link is down there. Thanks! [[alternative HTML version deleted]] __ R-help@r-projec

Re: [R] syntax for batching rbind process

It works perfectly now. Thank you all. Toby On 18 Aug 2010 at 15:04, Erik Iverson wrote: > > > Toby Gass wrote: > > Thank you for the suggestions for the more efficient code. The > > problem remains, however, that the final dataframe does not contain > > the correct values. So, in the cas

Re: [R] syntax for batching rbind process

Toby Gass wrote: Thank you for the suggestions for the more efficient code. The problem remains, however, that the final dataframe does not contain the correct values. So, in the case of the code you suggested, imp <- lapply(test, read.csv) do.call(rbind, imp) imp does contain all the da

Re: [R] syntax for batching rbind process

Thank you for the suggestions for the more efficient code. The problem remains, however, that the final dataframe does not contain the correct values. So, in the case of the code you suggested, imp <- lapply(test, read.csv) do.call(rbind, imp) imp does contain all the data from each datafram

Re: [R] Fwd: \ell symbol (log-likelihood)

Hi On 19/08/2010 4:25 a.m., David Winsemius wrote: I sent this privately to ivo welch yesterday, and he thinks it might be useful to someone else as well. Since I'm on a Mac the screen device is quartz(): > quartz() > plot( c(0,1), c(0,1) ); > text( 0.5, 0.5, "\u2113" ) # and then Fil

Re: [R] Yet another memory limit problem

On 08/18/2010 11:19 AM, Tim Clark wrote: > Dear List, ... > I would appreciate some help. I think I set the target field correctly, both > memory.size and memory.limit indicate I have over 2G of memory, yet I can't > allocate 831.3Mb? It just doesn't make sense to me. That shouldn't by itself

[R] R Courses***September - October 2010 schedule by XLSolutions Corp

We've scheduled R courses for September-October 2010 in Washington DC, San Francisco, Seattle, New York City and other USA cities http://www.xlsolutions-corp.com/Rcourses September-October 2010 *** R/S: Programming Essentials *** R Fundamentals and Programming Techniques *** R/S-PLUS Functions b

Re: [R] Again: wget parameters

On Aug 18, 2010, at 2:58 PM, christiaan pauw wrote: Hi Ricardo and everybody In this old post to rhelp you say that the problem was solved but not what the sollution was. I have the same problem now. I want to read a Excel file from a google site that has restricted access so I need to giv

Re: [R] dimnames

On 08/18/2010 06:20 PM, Joshua Wiley wrote: > matrix(1:4, nrow = 4, ncol = 1, dimnames = list(NULL, "true value")) > or even cbind("true value"=1:4) -- Peter Dalgaard Center for Statistics, Copenhagen Business School Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com _

[R] Again: wget parameters

Hi Ricardo and everybody In this old post to rhelp you say that the problem was solved but not what the sollution was. I have the same problem now. I want to read a Excel file from a google site that has restricted access so I need to give a username and password at some point. I could not find an

Re: [R] syntax for batching rbind process

Toby Gass wrote: Dear helpeRs, I am attempting to read in a series of csv files so I can bind them into one large dataframe. I have written the following script: test <- list.files(".", pattern = "csv") #lline 1 imp <- list()#line 2 for (i in 1:length(t

[R] syntax for batching rbind process

Dear helpeRs, I am attempting to read in a series of csv files so I can bind them into one large dataframe. I have written the following script: test <- list.files(".", pattern = "csv") #lline 1 imp <- list()#line 2 for (i in 1:length(test)) {

Re: [R] dimnames

On Aug 18, 2010, at 12:38 PM, David Winsemius wrote: On Aug 18, 2010, at 11:58 AM, Jimmy Söderly wrote: Hi everybody, I wanted to name a column vector (a variable). This is what I did : try<-matrix(1:4,nrow=4,ncol=1) ("try" is a bad name for an object because it is a useful function n

[R] Help with ANOVA

Hello everyone, I am doing an analysis on some data for a project, and do not have much background in either R or statistics, but have been pointed in the direction of doing ANOVA and then multiple comparisons on the data to give me the information I need. I was wondering if some of you could hel

[R] error with bitmap

Hi All I still have a problem with “bitmap” function in R, I have downloaded ghostscript 8.71 and added the path to the executable to the path environment variable, by going to the control panel , system, advanced system settings, added C:\Program Files(x86)\gs\gs8.71\bin to the path variable. I

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

On Aug 18, 2010, at 1:19 PM, Johan Jackson wrote: Hi all, Thanks for the replies (including off list). I have since resolved the discrepant results. I believe it has to do with R's scoping rules - I had an object called 'labs' and a variable in the dataset (DATA) called 'labs', and app

Re: [R] what does it mean when my main effect 'disappears' when using lme4?

Hi all, Thanks for the replies (including off list). I have since resolved the discrepant results. I believe it has to do with R's scoping rules - I had an object called 'labs' and a variable in the dataset (DATA) called 'labs', and apparently (to my surprise), when I called this: lmer(Y~X + (1|

Re: [R] Problem using 'get'

Try this: get('datasplit')$'3'$Area On Wed, Aug 18, 2010 at 1:53 PM, Johnny Tkach wrote: > Hi all, > > I have loaded a table called 'data' and used the split command to create a > new table called 'datasplit' based on the variable 'a.ImageNumber' in > 'data'. a.ImageNumber just specifies the im

[R] Problem using 'get'

Hi all, I have loaded a table called 'data' and used the split command to create a new table called 'datasplit' based on the variable 'a.ImageNumber' in 'data'. a.ImageNumber just specifies the image number that certain measurements are associated with (so there could be 330 measurements from

[R] Fwd: \ell symbol (log-likelihood)

I sent this privately to ivo welch yesterday, and he thinks it might be useful to someone else as well. Since I'm on a Mac the screen device is quartz(): > quartz() > plot( c(0,1), c(0,1) ); > text( 0.5, 0.5, "\u2113" ) # and then File/Save As/ -- David. Begin forwarded message: From: Da

Re: [R] dimnames

On Aug 18, 2010, at 11:58 AM, Jimmy Söderly wrote: Hi everybody, I wanted to name a column vector (a variable). This is what I did : try<-matrix(1:4,nrow=4,ncol=1) ("try" is a bad name for an object because it is a useful function name.) attr(try,"dimnames")<-list(NULL,"true value")

Re: [R] question about unwanted addresses in contact list

It doesn't happen to other lists that I belong to, just the R lists. But I will check. Thanks for the tip. Best, Scott On Aug 18, 2010, at 12:23 PM, Erik Iverson wrote: > > > Scott Compton wrote: >> Hi, >> It looks like about 200+ addresses from people who post on this list >> have been adde

Re: [R] question about unwanted addresses in contact list

Scott Compton wrote: Hi, It looks like about 200+ addresses from people who post on this list have been added to my contact list without my permission. Has this happened to other people? I use Mac OS. This is quite disconcerting and suggests a virus floating around on this distribution list.

Re: [R] dimnames

On Wed, Aug 18, 2010 at 8:58 AM, Jimmy Söderly wrote: > Hi everybody, > > I wanted to name a column vector (a variable). > > This is what I did : > > try<-matrix(1:4,nrow=4,ncol=1) > attr(try,"dimnames")<-list(NULL,"true value") > > Is it OK ? Sure. You can also do: dimnames(try) <- list(NULL,"

[R] question about unwanted addresses in contact list

Hi, It looks like about 200+ addresses from people who post on this list have been added to my contact list without my permission. Has this happened to other people? I use Mac OS. This is quite disconcerting and suggests a virus floating around on this distribution list. Thanks, Scott __

[R] Coinertia randtest

Hello! I dunno why, but I cannot make randtes.coinertia() from ade4 package working. I have two nice distance matrices (Euclidean): > dist1 1 2 3 4 5 6 7 2 2.5776799 3 1.7892825 1.0637487 4 1.0557991 2.4270728 2.0626604 5 1.6745483 4

[R] dimnames

Hi everybody, I wanted to name a column vector (a variable). This is what I did : try<-matrix(1:4,nrow=4,ncol=1) attr(try,"dimnames")<-list(NULL,"true value") Is it OK ? Thanks for your help, Jimmy [[alternative HTML version deleted]] __ R-

Re: [R] functions and multiple levels

On Wed, 18 Aug 2010, chris20 wrote: Hi, I am trying to write a function; I want to subtract the mean of each class in level 2 from the mean of each class in level 1 and square the answer, eg. level.1 level.2 observation 1 1 0.5 1

Re: [R] svyquantile w/ svyby is returning an error

On Tue, 17 Aug 2010, Marc Zodet wrote: svymean w/ svyby is working for me... svyby(~visitcnt, ~agegrp3.f, svymean, design=svydes) agegrp3.f visitcnt se.visitcnt 18-44 18-44 8.72 0.4953235 45-64 45-64 10.131555 0.5347806 65+ 65+ 9.588802 0.4323629 svyquantile

Re: [R] ifelse command

Hi: Here's the problem I had with the OP's function: On Wed, Aug 18, 2010 at 4:26 AM, David Winsemius wrote: > > On Aug 18, 2010, at 5:18 AM, Philip Wong wrote: > > >> hello people, >> I want to make a biased dice using the sample() function and print out the >> results after n number of runs, I

Re: [R] Combining values from a data.frame row with matrix rownames.

Terrific! Many thanks, Josh, you put me on the right path with the tip re dput and the remark about factors. Using "stringsAsFactors=FALSE" in my read.table and then "IndexList <- c( unlist(seed_panel[ i, ]), rownames( Control[ j, ] ) )" did the trick. Your help is hugely appreciated. Thanks aga

Re: [R] Linear regression equation and coefficient matrix

Hi, Thanks, the cor() works. Regarding the simple linear regression equation (mainly, the slope parameter) and r2. I think I was not writing it well. I need to do it just for the columns. If I have a, b, c, d columns I wish to compute the relation of there data, e.g., between a-b, a-c, a-d, b-a

Re: [R] How to Perform CCA in ??!! Help please

Hello Ista, I recently posted a new message on my own post that should clarify what you are talking about. Let me know if you can't see it. I been having problems post things on the forum. thanks On Wed, Aug 18, 2010 at 6:57 AM, Ista Zahn-2 [via R] < ml-node+2329880-957103028-137...@n4.nabble.co

[R] functions and multiple levels

Hi, I am trying to write a function; I want to subtract the mean of each class in level 2 from the mean of each class in level 1 and square the answer, eg. level.1 level.2 observation 1 1 0.5 1 1 0.2 1

Re: [R] mtext coordinates

Thank, it is quite what I was looking for. But with your code, when I tried to change the size of the plot window, the segments are moving... I don't know if you see what I want, just try to expand the height of the window. Is there any solution for this ? 2010/8/18 Jim Lemon > On 08/18/2010 0

Re: [R] rJava: System.exit(1) in Java causes R crash

On Aug 18, 2010, at 12:18 AM, Sigal Blay wrote: > I'm writing a package that successfully calls Java from R using rJava. > When R sends my Java function bad data (through .jcall()), > I get the error message details (when I run it on a console) > and then: > > Java requested System.exit(1), clo

Re: [R] functions and multiple levels

Hi Chris, Try ?ave will help you. Anyway, I guess you are computing a statistic. strs <- " level.1 level.2 observation 1 1 0.5 1 1 0.2 1 2 0.6 1 2 0

Re: [R] Missing values in function calls

?missing -- 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 vo...@cbox.cz > Sent: Tuesday, August 17, 2010 6:22

Re: [R] Linear regression equation and coefficient matrix

Hmm, after reading one of your other posts, I am thinking you may *just* want all pairwise combinations. This worked for me: # Create a sample data frame with 60 named columns x <- data.frame(matrix(rnorm(1200), ncol = 60, dimnames = list(NULL, paste("Col", 1:60, sep='' # calculate the corre

  1   2   >