[R] calculating p-values from t-values for a Bonferroni adjustment

2009-07-27 Thread NativeBuff2003
I am performing a sequential bonferroni adjustment on the results of an ANCOVA but the equation I have for calculating p-values from the t-values is not working. I can't seem to find it anywhere else. This is the code I have now: 2*(1-pt(t,df)) where t=t-value and df=degrees of freedom from the

Re: [R] Writing to a UDP server from R?

2009-07-27 Thread Warren Young
Saptarshi Guha wrote: This does not work u - socketConnection('localhost',9000) UDP and TCP are entirely different worlds. They don't share the same port space; TCP port 9000 and UDP port 9000 are different things. Plus, UDP is connectionless, so you can't connect to a UDP port. All

Re: [R] local regression using loess

2009-07-27 Thread Bert Gunter
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of cindy Guo Sent: Monday, July 27, 2009 4:06 PM To: r-help@r-project.org Subject: [R] local regression using loess Hi, All, I have a dataset with binary response ( 0 and 1) and some

Re: [R] local regression using loess

2009-07-27 Thread cindy Guo
Hi, Bert, Thanks for the response. But then in this case, can I use loess to fit the data? If yes, then how to interpret the results? Cindy On Mon, Jul 27, 2009 at 4:32 PM, Bert Gunter gunter.ber...@gene.com wrote: -Original Message- From: r-help-boun...@r-project.org

[R] Get bearing for cartesian coordinates

2009-07-27 Thread pecardoso
Is it possible to get bearing in degrees from Cartesian (not lat long) coordinates? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] skip plot/blank plot on purpose (multi-plot question)

2009-07-27 Thread Greg Snow
Also look at the frame and plot.new functions. -Original Message- From: Mark Knecht markkne...@gmail.com To: Bert Gunter gunter.ber...@gene.com Cc: r-help r-help@r-project.org Sent: 7/27/09 1:56 PM Subject: Re: [R] skip plot/blank plot on purpose (multi-plot question) On Mon, Jul 27,

[R] Multiple graphs

2009-07-27 Thread Data Analytics Corp.
Hi, I wrote a simple master function, run(), that has inside six qplot functions. The goal is to type run() and have all six graphs appear as separate windows so that I can copy them into PowerPoint for a client. When I type run(), only the last graph appears, the first five apparently

Re: [R] Splitting matrix into several small matrices

2009-07-27 Thread jim holtman
Here is one way to do it: x V1 V2 V3 V4 V5 [1,] 13 1 1 1 1 [2,] 12 0 0 0 0 [3,] 8 1 0 1 1 [4,] 9 0 1 0 0 [5,] 10 1 1 1 1 [6,] 3 0 1 0 0 [7,] 3 1 0 1 1 [8,] 6 1 1 1 1 # create the row indices based on columns 2-3 x.i - split(seq(nrow(x)), paste(x[,2],

Re: [R] Making a sub data.frame

2009-07-27 Thread jim holtman
Try this: x ID VAR1 1 11 blaaal 2 121 blalda 3 121 adada 4 234 baada 5 231 ddaaa 6 231 baada idm - c(121,234,231) subset(x, ID %in% idm) ID VAR1 2 121 blalda 3 121 adada 4 234 baada 5 231 ddaaa 6 231 baada On Mon, Jul 27, 2009 at 4:40 PM, desperides...@gmail.com wrote:

Re: [R] Multiple graphs

2009-07-27 Thread Cedrick Johnson
Try dev.new() after each graph is generated. Example code from my setup (adapted a bit to fit your situation): Alternate Method 1: Generate images and insert into PPT png(file = chart1.png, width=800, height=800) #do something pretty dev.off() Alternate Method 2: Generate images with new

Re: [R] Confirmatory factor analysis problems using sem package (works in Amos)

2009-07-27 Thread Solomon Messing
Dear John, Would it possible to use a different optimizer with the sem package? Perhaps optim(..., method = c(Nelder-Mead, BFGS, CG, L-BFGS-B, SANN),...) for example? Thank you very much, -Solomon -Original Message- From: John Fox [mailto:j...@mcmaster.ca] Sent: Friday, May 22,

Re: [R] local regression using loess

2009-07-27 Thread Ryan
Hi, All, I have a dataset with binary response ( 0 and 1) and some numerical covariates. I know I can use logistic regression to fit the data. But I want to consider more locally. So I am wondering how can I fit the data with 'loess' function in R? And what will be the response:

Re: [R] Forumla format?

2009-07-27 Thread Bernd Kreuss
Steve Lianoglou wrote: Is it possible to build up your formula as a string, and then convert to formula w/ as.formula? what about simply using it this way instead: svm(label ~ ., data=mydata[,-c(22,23,25,31)]) __ R-help@r-project.org mailing list

Re: [R] Assigning rank based on total count

2009-07-27 Thread ws
I am using ACS micro data (PUMS) with one of the columns as a factor for the place of birth (POBPF). I would like to create a column (POBR) containing a rank corresponding to the place of the observation in the POBPF rankings. I wrote a blog entry on my solution:

Re: [R] Confirmatory factor analysis problems using sem package (works in Amos)

2009-07-27 Thread John Fox
Dear Solomon, When I originally programmed the sem() function, I used optim() and experimented with the different methods provided, settling on BFGS as the default. Shortly after that, I compared optim() to nlm() for a range of problems and found that the latter performed better. It would not be

Re: [R] How should i change the SAS Codes into R Codes?

2009-07-27 Thread zhijie zhang
Thanks a lot. Matt's code works very well. Could i use several 'for' arguments to get the results as what i think first? 2009/7/28 John Kane jrkrid...@yahoo.ca I see Matt Aldridge has given you the answers to your specific questions. If you are used to using SAS you might find Bob

Re: [R] local regression using loess

2009-07-27 Thread cindy Guo
Hi, Ryan, Thank you for the information. I tried it. But there are some error messages. When I use fit - locfit(Y~X1*X2,family='binomial'), the error message is error lfproc(x, y, weights = weights, cens = cens, base = base, geth = geth, : compparcomp: parameters out of bounds And when I use

Re: [R] calculating p-values from t-values for a Bonferroni adjustment

2009-07-27 Thread Ben Bolker
NativeBuff2003 wrote: I am performing a sequential bonferroni adjustment on the results of an ANCOVA but the equation I have for calculating p-values from the t-values is not working. I can't seem to find it anywhere else. This is the code I have now: 2*(1-pt(t,df)) where t=t-value and

Re: [R] calculating p-values from t-values for a Bonferroni adjustment

2009-07-27 Thread NativeBuff2003
Ben Bolker wrote: NativeBuff2003 wrote: I am performing a sequential bonferroni adjustment on the results of an ANCOVA but the equation I have for calculating p-values from the t-values is not working. I can't seem to find it anywhere else. This is the code I have now:

Re: [R] dataframe to list conversion

2009-07-27 Thread Bill.Venables
All you need is csum - function(m) sapply(m, sum) (in which case making csum a function does not achieve very much). If for some reason you want to hang on to your code and just modify the last line (though I cannot think why, but still...) you could do csum - function(m) { a -

Re: [R] dataframe to list conversion

2009-07-27 Thread Jorge Ivan Velez
Dear voidobscura, Try also Csum - function(x) apply(x, 2, sum) Csum(m) HTH, Jorge On Mon, Jul 27, 2009 at 11:03 AM, voidobscura wrote: Hi all, I have been experimenting with writing my own matrix column sum function. I want it to return a list. csum-function(m) { a =

[R] storing output in html or pdf table format.

2009-07-27 Thread Albert EINstEIN
Hi every one, Thanks for every one who are all supporting to us. we want some clarification on output in R. I have generated summary statistics output for dataset (E.g. sales) in output window. Now i want to store that output in a html or pdf in a table format. if possible can any one provide

Re: [R] storing output in html or pdf table format.

2009-07-27 Thread cls59
Albert EINstEIN wrote: Hi every one, Thanks for every one who are all supporting to us. we want some clarification on output in R. I have generated summary statistics output for dataset (E.g. sales) in output window. Now i want to store that output in a html or pdf in a table format. if

<    1   2