[R] Digest, Vol90, Issue18 does not display properly

2010-08-18 Thread Jianyun Wu
Dear R team, I received Digest, Vol90, Issue18 today. However there are bunchs of threads not displaying properly. Could you please check it out?? Thanks and Regards Fred *After Message 15, :* Date/Time: 2010-08-16 14:14:57 +0300 OS Version: 10.6.4 (Build 10F569) Architecture:

Re: [R] ask for help about netcdf

2010-08-18 Thread Prof Brian Ripley
Please learn that capitalization is important both in R (sic) and in effective communication. There are CRAN packages RNetCDF, ncdf and ncdf4 on CRAN to read/write NetCDF (sic) format files from R. On Wed, 18 Aug 2010, mat chueng wrote: hi,r project I am a new user who has some questions

Re: [R] using color palette in persp

2010-08-18 Thread szisziszilvi
I remembered wrong, the problem rised int the next line to the mentioned one, z was wrong, and I didn't see way. But new day, new point of wiew. For anybody's interest here is the solution. The original try was this: x - seq(-1.95, 1.95, length = 30) y - seq(-1.95, 1.95, length = 35) z -

[R] How to read ANOVA output

2010-08-18 Thread Stephen Liu
Hi folks, Where can I find document re how to read anova output? Google found many of them. But seemingly non of them can explain to me following output:- tabA = c(5.67, 5.67, 5.55, 5.57) tabB = c(5.75, 5.47, 5.43, 5.45) tabC = c(4.74, 4.45, 4.65, 4.94) tabs = data.frame(tabA, tabB,

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 05:42 +0800, elaine kuo wrote: Thank you. Most of the answers solved the puzzles. Q2 I tried to display sub-model with only temp_ran using the code below but failed. Please kindly suggest the potential failure cause.

Re: [R] How to read ANOVA output

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 00:42 -0700, Stephen Liu wrote: Hi folks, Where can I find document re how to read anova output? Google found many of them. But seemingly non of them can explain to me following output:- tabA = c(5.67, 5.67, 5.55, 5.57) tabB = c(5.75, 5.47, 5.43, 5.45)

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 08:51 +0100, Gavin Simpson wrote: On Wed, 2010-08-18 at 05:42 +0800, elaine kuo wrote: Thank you. Most of the answers solved the puzzles. snip / = Please explain why fitted lm is better for GLM. Seriously? A GLM specified as glm(, family = gaussian) is the

Re: [R] how to paste a vector to expression?

2010-08-18 Thread Gavin Simpson
On Tue, 2010-08-17 at 13:49 -0700, Kay Cichini wrote: ..thank you! this generally works well - but after saving the result as jpeg i noticed that the printed = used with the expression looks different from the = signs used for the same graph for another simple text annotation (i.e. the

Re: [R] How to read ANOVA output

2010-08-18 Thread Ted Harding
On 18-Aug-10 07:42:23, Stephen Liu wrote: Hi folks, Where can I find document re how to read anova output? Google found many of them. But seemingly non of them can explain to me following output:- tabA = c(5.67, 5.67, 5.55, 5.57) tabB = c(5.75, 5.47, 5.43, 5.45) tabC = c(4.74, 4.45,

[R] reading lmer table

2010-08-18 Thread Nicola Spotorno
Dear all, I'm quite new in R and especially with linear mixed effects models and I'm not completely sure to read the lmer table in the right way. for example: head(march.f) fam subjID Cond Code reg total first second log.total log.second cat 3

[R] ifelse command

2010-08-18 Thread Philip Wong
hello people, I want to make a biased dice using the sample() function and print out the results after n number of runs, I've successfully generated the dice using the following command: mydie2-function(n=1000,y=NULL,...){ for(i in 1:n){ x-sample(1:6,n,replace=TRUE,prob=c(1,1,2,3,2,1)/10)

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
Hmmm. Thinking some more, I might not have answered your (unstated) question. What do your mean by GLM? = I meant generalized linear model as well. Thanks for the references. The first one was mentioned first in my life time after keeping asking the same question. I mean the Generalized

Re: [R] how to paste a vector to expression?

2010-08-18 Thread Kay Cichini
..of course that works, too. thanks a lot, kay Zitat von Gavin Simpson [via R] ml-node+2329504-916404301-40...@n4.nabble.com: On Tue, 2010-08-17 at 13:49 -0700, Kay Cichini wrote: ..thank you! this generally works well - but after saving the result as jpeg i noticed that the

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
Please suggest how to define subset in my case How would I know? I still haven't seen your data. You seem to be mistaken on what is and is not included in your model and you fitted it. What hope do we have...? However, given the model 'mig.stds' from above in this email: mig.stds

Re: [R] Weird differing results when using the Wilcoxon-test

2010-08-18 Thread Cedric Laczny
I was able to trace down the unexpected behavior to the following line SIGMA - sqrt((n.x * n.y/12) * ((n.x + n.y + 1) - sum(NTIES^3 - NTIES)/((n.x + n.y) * (n.x + n.y - 1 My calculations of the Z-score for the normal approximation where based on using the

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 17:37 +0800, elaine kuo wrote: snip / = I posted it for help, after following the manual with the command dredge but receiving an error message two days ago. command target-dredge(mig.stds, subset = temp_max) error in eval(expr,

Re: [R] ifelse command

2010-08-18 Thread Paul Hiemstra
Hi Philip, Why do you want to use ifelse? It is not an alternative to an if then else structure, but a way of iterating over a vector with a statement, returning one value it is true, another if false. An example: x = runif(100) y = ifelse(x 0.5, Larger, Smaller) x y cheers, Paul On

Re: [R] ifelse command

2010-08-18 Thread Joshua Wiley
Dear Philip, Trying to use 5 nested ifelse() statements will not be the same as a series of if statements, and will also be quite messy. Also, consider these lines of your code: x - sample(1:6,n,replace=TRUE,prob=c(1,1,2,3,2,1)/10) x = runif(n) First you assign n samples of 1:6 to 'x', then

[R] mtext coordinates

2010-08-18 Thread David Hajage
Hello, I am using mtext to write some text below a graph: plot(1) mtext(foo, side = 1, line = 2, at = seq(0.6, 1.6, 0.2)) I would like to draw something near the foo texts, for example a segment. For this, I need to know the coordinates of the text ploted by mtext function. I spent a few hours

[R] Using objectname in function

2010-08-18 Thread JesperHybel
Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey-c(0,0,0,1,1,1) Wax-c(1,0,1,0,1,0) f-function(x,y){ table(x,y) } f(Monkey,Wax) so that the printout is not y x 0 1 0 1 2 1 2 1 but Wax Monkey 0 1 0 1 2 1 2 1 -- View

Re: [R] mtext coordinates

2010-08-18 Thread Jim Lemon
On 08/18/2010 08:40 PM, David Hajage wrote: Hello, I am using mtext to write some text below a graph: plot(1) mtext(foo, side = 1, line = 2, at = seq(0.6, 1.6, 0.2)) I would like to draw something near the foo texts, for example a segment. For this, I need to know the coordinates of the text

Re: [R] Using objectname in function

2010-08-18 Thread Gabor Grothendieck
On Wed, Aug 18, 2010 at 7:10 AM, JesperHybel jesperhy...@hotmail.com wrote: Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey-c(0,0,0,1,1,1) Wax-c(1,0,1,0,1,0) f-function(x,y){ table(x,y) } f(Monkey,Wax) so that the printout is not   y x   0 1  

Re: [R] mtext coordinates

2010-08-18 Thread David Winsemius
On Aug 18, 2010, at 6:40 AM, David Hajage wrote: Hello, I am using mtext to write some text below a graph: plot(1) mtext(foo, side = 1, line = 2, at = seq(0.6, 1.6, 0.2)) I would like to draw something near the foo texts, for example a segment. For this, I need to know the coordinates of

Re: [R] Using objectname in function

2010-08-18 Thread Barry Rowlingson
On Wed, Aug 18, 2010 at 12:10 PM, JesperHybel jesperhy...@hotmail.com wrote: f-function(x,y){ table(x,y If you look at the code for 'plot', you'll see it does that with some deparse/substitute magic. Then use the dnn option to table to set the names: f-function(x,y){

Re: [R] ifelse command

2010-08-18 Thread David Winsemius
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've successfully generated the dice using the following command: mydie2-function(n=1000,y=NULL,...){ for(i in 1:n){

Re: [R] Using objectname in function

2010-08-18 Thread Duncan Murdoch
On 18/08/2010 7:10 AM, JesperHybel wrote: Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey-c(0,0,0,1,1,1) Wax-c(1,0,1,0,1,0) f-function(x,y){ table(x,y) } f(Monkey,Wax) so that the printout is not y x 0 1 0 1 2 1 2 1 but Wax Monkey

Re: [R] Using objectname in function

2010-08-18 Thread Gabor Grothendieck
On Wed, Aug 18, 2010 at 7:23 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Wed, Aug 18, 2010 at 7:10 AM, JesperHybel jesperhy...@hotmail.com wrote: Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey-c(0,0,0,1,1,1) Wax-c(1,0,1,0,1,0)

Re: [R] [solved] Weird differing results when using the Wilcoxon-test

2010-08-18 Thread Cedric Laczny
By using the help page for wilcox.test and looking in the there-mentioned book (Hollander and Wolfe), I could find the formula also in the book. People that might run into the same problem: The formula can also be found here, at least if you have access to it via your university or insitution,

Re: [R] ifelse command

2010-08-18 Thread JesperHybel
Or try: datbig-numeric(1000) for (i in 1:1000){ datbig[i] - mean(sample(1:6, 1000, prob = c(1, 1, 2, 3, 2, 1)/10, replace = TRUE)) hist(datbig,breaks=FD) -- View this message in context: http://r.789695.n4.nabble.com/ifelse-command-tp2329538p2329711.html Sent from the R help mailing list

Re: [R] ifelse command

2010-08-18 Thread JesperHybel
sorry missed an } datbig-numeric(1000) for (i in 1:1000){ datbig[i] - mean(sample(1:6, 1000, prob = c(1, 1, 2, 3, 2, 1)/10, replace = TRUE)) hist(datbig,breaks=FD) } -- View this message in context: http://r.789695.n4.nabble.com/ifelse-command-tp2329538p2329717.html Sent from the R help

Re: [R] Delete rpart/mvpart cross-validation output

2010-08-18 Thread Tal Galili
Hi Marie and Gavin, I do remember there is some command doing silent, so to suppress output (as if directing it to sink, but not really directing it anywhere). The problem is I don't remember the command at the moment - but some searching might yield you results. Cheers, Tal

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
A cause other than data based on standardized regression was identified. It is that the manual command added with target - at the left hand side. C1 did not work but C2 did. C1 target-dredge(mig.stds, subset = temp_max) C2 dredge(mig.stds, subset = temp_max) Elaine On Wed, Aug 18, 2010

Re: [R] Delete rpart/mvpart cross-validation output

2010-08-18 Thread Gabor Grothendieck
On Fri, Aug 13, 2010 at 1:52 PM, Marie-Hélène Ouellette mariehele...@gmail.com wrote: Dear all, I was wondering if there is a simple way to avoid printing the multiple cross-validation automatic output to the console of recursive partitionning functions like rpart or mvpart. For example...

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 21:11 +0800, elaine kuo wrote: A cause other than data based on standardized regression was identified. It is that the manual command added with target - at the left hand side. C1 did not work but C2 did. C1 target-dredge(mig.stds, subset = temp_max) C2

Re: [R] Sweave

2010-08-18 Thread Randall Wrong
Thanks Karen and Abhijit. I have read the section 4.1.2 of the Sweave user manual. Actually the manual lacks example code. I would like to change the option in includegraphics, not to change the true sizes of the pictures. This is what I usually do : \begin{figure} \centering

Re: [R] Delete rpart/mvpart cross-validation output

2010-08-18 Thread Matt Shotwell
Or, if using GNU Linux or other UNIX-like system: sink(/dev/null) # Issue commands sink() -Matt On Wed, 2010-08-18 at 09:14 -0400, Gabor Grothendieck wrote: On Fri, Aug 13, 2010 at 1:52 PM, Marie-Hélène Ouellette mariehele...@gmail.com wrote: Dear all, I was wondering if there is a

Re: [R] Sweave

2010-08-18 Thread Abhijit Dasgupta
No, the \setkeys statement should be in the main body of the Sweave file, not in the R code part. On Aug 18, 2010 9:20 AM, Randall Wrong randall.wr...@gmail.com wrote: Thanks Karen and Abhijit. I have read the section 4.1.2 of the Sweave user manual. Actually the manual lacks example code. I

Re: [R] Delete rpart/mvpart cross-validation output

2010-08-18 Thread Henrique Dallazuanna
Try this also: invisible(capture.output(x - mvpart(data.matrix(spider[,1:12])~herbs+reft+moss+sand+twigs+water,spider,xv=1se,xvmult=100))) x On Fri, Aug 13, 2010 at 2:52 PM, Marie-Hélène Ouellette mariehele...@gmail.com wrote: Dear all, I was wondering if there is a simple way to avoid

Re: [R] Delete rpart/mvpart cross-validation output

2010-08-18 Thread Gabor Grothendieck
On Wed, Aug 18, 2010 at 9:14 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Fri, Aug 13, 2010 at 1:52 PM, Marie-Hélène Ouellette mariehele...@gmail.com wrote: Dear all, I was wondering if there is a simple way to avoid printing the multiple cross-validation automatic output to the

[R] Controlling the dates in R;

2010-08-18 Thread Mangalani Peter Makananisa
Dear R gurus, I am currently using the chron package for date manipulation in R and I am able to control the date which take the format - 1998/07/03, and - 2006-09-15 (date without time) Library(chron) As.Date(2006-09-15, format = '%Y-%m-%d') However I have a

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

2010-08-18 Thread stompper33
Ok..So here is an example of some of my code Sites,P1,P2,P3,P4,P5,P6,P7,S1,S2,S3,S4,S4,S5 trap01,0,1,0,1,1,0,1,1,0,0,0,0 trap02,0,1,0,0,0,0,0,0,1,1,0,0 trap03,0,1,0,1,1,1,1,1,0,0,0,1 trap04,0,1,0,0,0,0,0,0,1,1,0,0 trap05,0,1,0,0,0,1,1,1,1,0,0,1 So here is a SMALL example of what the code looks

[R] Yet another memory limit problem

2010-08-18 Thread Tim Clark
Dear List, I have read and read and still don't get why I am getting a memory issue.  I am using a Samsung PC running Windows 7.  I have set memory in the target field:  C:\Program Files (x86)\R\R-2.11.1\bin\Rgui.exe --max-mem-size=3G But when I try a simple plot I get:   

[R] Rolling window linear regression

2010-08-18 Thread siddharth . garg85
Hi Does there exists an efficient way of performing linear regression on rolling windows in R. The exact problem is: We have a dataset of length l. The window size is w. Now, I perform linear regression on window i to (i+w) . Using this model can I perform linear regression over window (i+1)

[R] recurrent events

2010-08-18 Thread avsha38
Hello everyone, I have a question about using frailtypack to run recurrent events data. I have a file with 1,580 events. The file structure fits the andersen gill model. When I run the file with all records, the R stops and shuts down. Once I run only part of the file, about 400 events, R runs

Re: [R] creation package

2010-08-18 Thread anderson nuel
Dear r-help, No, I find errors in the file C:/Rp/namepackage.Rcheck/00check.txt : * using log directory 'C:/Rpackage/namepackage.Rcheck' * using R version 2.10.1 (2009-12-14) * using session charset: ISO8859-1 * checking for file 'namepackage/DESCRIPTION' ... OK * checking extension type ...

[R] Array Data

2010-08-18 Thread navishkumarb
Hello Can any one let me know how to delete a value from an array and then push back rest of the remaining elements up into an array. Thanks -- View this message in context: http://r.789695.n4.nabble.com/Array-Data-tp2329643p2329643.html Sent from the R help mailing list archive at

Re: [R] Weird differing results when using the Wilcoxon-test

2010-08-18 Thread peter dalgaard
On Aug 18, 2010, at 11:55 AM, Cedric Laczny wrote: I was able to trace down the unexpected behavior to the following line SIGMA - sqrt((n.x * n.y/12) * ((n.x + n.y + 1) - sum(NTIES^3 - NTIES)/((n.x + n.y) * (n.x + n.y - 1 My calculations of the Z-score

Re: [R] Yet another memory limit problem

2010-08-18 Thread ONKELINX, Thierry
Dear Tim, Are mydat$Date and mydat$Time what you think they are? What output do you get from str(mydat)? HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie Kwaliteitszorg

Re: [R] Fwd: R SOFTWARE

2010-08-18 Thread Erik Iverson
Im busy doing an exercise on Foundation of Statistical Description and Analysis,Please check the attached exercise on page 3.In this exercise the first step is to check the content of Limestone.data so I entered the following command on R console

[R] combinations

2010-08-18 Thread Maas James Dr (MED)
I would appreciate any suggestions on which function to use to write subsequent functions analysing combinations of treatments. This refers to experimental trials of medical treatments. I want to write routines to analyse various comparisons (combinations) So if 5 treatments are

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

2010-08-18 Thread Ista Zahn
If you read your email back to yourself, you will notice that you never actually describe the problems you are having. There are people on this list who can help you, but you need to tell them what you need help with. We don't need your whole data set and script (in fact sending the whole shebang

Re: [R] Array Data

2010-08-18 Thread Erik Iverson
On 08/18/2010 06:16 AM, navishkumarb wrote: Hello Can any one let me know how to delete a value from an array and then push back rest of the remaining elements up into an array. Thanks Reproducible examples help us understand what you mean. x - 1:10 x[-1] Is that what you mean?

[R] Tcl/Tk interface - autolaunch point-and-click problem

2010-08-18 Thread Joanna Zhuang
Dear all, I'm trying to run function qvalue that requires Tcl/Tk interface. But after loaded the tcltk package, it would automatically launch point-and-click that forces me to run analysis in a new window whenever I want to use the function. I would be grateful if anyone knows how to remove

[R] Linear regression equation and coefficient matrix

2010-08-18 Thread ashz
Hi, I have 20*60 data matrix (with some NAs) and I wish to perfom a Pearson correlation coefficient matrix as well as simple linear regression equation and coefficient of determination (R2) for every possible combination. Any tip/idea/library/script how do to so. Thanks, As hz -- View

Re: [R] combinations

2010-08-18 Thread Jorge Ivan Velez
Hi Jim, How about expand.grid()? expand.grid(c('t01','t02','t03','t04','t05'), 't04') Var1 Var2 1 t01 t04 2 t02 t04 3 t03 t04 4 t04 t04 5 t05 t04 expand.grid(c('t01','t02','t03','t04','t05'), 't01') Var1 Var2 1 t01 t01 2 t02 t01 3 t03 t01 4 t04 t01 5 t05 t01 HTH, Jorge

Re: [R] Controlling the dates in R;

2010-08-18 Thread Gabor Grothendieck
On Wed, Aug 18, 2010 at 3:33 AM, Mangalani Peter Makananisa pmakanan...@sars.gov.za wrote: Dear R gurus, I am currently using the chron package for date manipulation in R and I am able to control  the date which take the format -          1998/07/03, and -          2006-09-15 (date

Re: [R] Fwd: R SOFTWARE

2010-08-18 Thread david.jessop
HI The problem is the file is called limestone.dat.txt. In Explorer in Windows you should go to Tools, Folder Options and under View you should find an option Hide Extensions for known file types - turn this off and then press Apply to all folders. Then you'll see the real file name. It's

Re: [R] Fwd: R SOFTWARE

2010-08-18 Thread Wu Gong
Hi Nokuzola, I guess the capital T. 'c:/temp/M5/limestone.dat' C:\Temp\M5\limestone.dat Regards, Wu - A R learner. -- View this message in context: http://r.789695.n4.nabble.com/Fwd-R-SOFTWARE-tp2329861p2329887.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] combinations

2010-08-18 Thread Dimitris Rizopoulos
one possibility is: treats - t(combn(c(t01, t02, t03, t04, t05), 2)) treats # extract only t04 treats[apply(treats == t04, 1, any), ] I hope it helps. Best, Dimitris On 8/18/2010 1:44 PM, Maas James Dr (MED) wrote: I would appreciate any suggestions on which function to use to write

Re: [R] Linear regression equation and coefficient matrix

2010-08-18 Thread Joshua Wiley
On Wed, Aug 18, 2010 at 6:09 AM, ashz a...@walla.co.il wrote: Hi, I have 20*60 data matrix (with some NAs) and I wish to perfom a  Pearson correlation coefficient matrix as well as simple linear regression equation The correlation matrix can be readily obtained by calling cor() on the entire

Re: [R] Sweave

2010-08-18 Thread Randall Wrong
I am sorry, I don't get it :-( What's the main body of the Sweave file ? Randall 2010/8/18 Abhijit Dasgupta adasgu...@araastat.com No, the \setkeys statement should be in the main body of the Sweave file, not in the R code part. On Aug 18, 2010 9:20 AM, Randall Wrong

Re: [R] Sweave

2010-08-18 Thread Duncan Murdoch
On 18/08/2010 10:24 AM, Randall Wrong wrote: I am sorry, I don't get it :-( What's the main body of the Sweave file ? Abhijit must have misread what you posted. You did put it outside of the R code part. I'd normally put it earlier in the file, so I got consistent sizing for all of the

Re: [R] creation package

2010-08-18 Thread Michael Dewey
At 10:27 18/08/2010, anderson nuel wrote: Dear r-help, No, I find errors in the file C:/Rp/namepackage.Rcheck/00check.txt : * using log directory 'C:/Rpackage/namepackage.Rcheck' * using R version 2.10.1 (2009-12-14) * using session charset: ISO8859-1 * checking for file

Re: [R] Linear regression equation and coefficient matrix

2010-08-18 Thread Joshua Wiley
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

Re: [R] Missing values in function calls

2010-08-18 Thread Greg Snow
?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 AM

Re: [R] functions and multiple levels

2010-08-18 Thread Wu Gong
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

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

2010-08-18 Thread Simon Urbanek
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), closing

Re: [R] mtext coordinates

2010-08-18 Thread David Hajage
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

[R] functions and multiple levels

2010-08-18 Thread chris20
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] How to Perform CCA in ??!! Help please

2010-08-18 Thread stompper33
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]

Re: [R] Linear regression equation and coefficient matrix

2010-08-18 Thread ashz
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,

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

2010-08-18 Thread Gil Gamesh
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 again!

Re: [R] ifelse command

2010-08-18 Thread Dennis Murphy
Hi: Here's the problem I had with the OP's function: On Wed, Aug 18, 2010 at 4:26 AM, David Winsemius dwinsem...@comcast.netwrote: 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

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

2010-08-18 Thread Thomas Lumley
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

Re: [R] functions and multiple levels

2010-08-18 Thread Charles C. Berry
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

[R] dimnames

2010-08-18 Thread Jimmy Söderly
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] Coinertia randtest

2010-08-18 Thread Petar Milin
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

[R] question about unwanted addresses in contact list

2010-08-18 Thread Scott Compton
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

Re: [R] dimnames

2010-08-18 Thread Joshua Wiley
On Wed, Aug 18, 2010 at 8:58 AM, Jimmy Söderly jimmy.sode...@gmail.com 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) -

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

2010-08-18 Thread Erik Iverson
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] question about unwanted addresses in contact list

2010-08-18 Thread Scott Compton
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 added to my

Re: [R] dimnames

2010-08-18 Thread David Winsemius
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) Is it

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

2010-08-18 Thread David Winsemius
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: David

[R] Problem using 'get'

2010-08-18 Thread Johnny Tkach
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

Re: [R] Problem using 'get'

2010-08-18 Thread Henrique Dallazuanna
Try this: get('datasplit')$'3'$Area On Wed, Aug 18, 2010 at 1:53 PM, Johnny Tkach johnny.tk...@utoronto.cawrote: 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

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

2010-08-18 Thread Johan Jackson
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 +

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

2010-08-18 Thread David Winsemius
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

[R] error with bitmap

2010-08-18 Thread kayj
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.

[R] Help with ANOVA

2010-08-18 Thread Sneeketeeke
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

Re: [R] dimnames

2010-08-18 Thread David Winsemius
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

[R] syntax for batching rbind process

2010-08-18 Thread Toby Gass
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] syntax for batching rbind process

2010-08-18 Thread Erik Iverson
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(test))

[R] Again: wget parameters

2010-08-18 Thread christiaan pauw
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

Re: [R] dimnames

2010-08-18 Thread Peter Dalgaard
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

Re: [R] Again: wget parameters

2010-08-18 Thread David Winsemius
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

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

2010-08-18 Thread Sue Turner
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

Re: [R] Yet another memory limit problem

2010-08-18 Thread Peter Dalgaard
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

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

2010-08-18 Thread Paul Murrell
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

Re: [R] syntax for batching rbind process

2010-08-18 Thread Toby Gass
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

Re: [R] syntax for batching rbind process

2010-08-18 Thread Erik Iverson
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

  1   2   >