Re: [R] Kolmogorov-Smirnov Test

2013-08-02 Thread peter dalgaard
On Aug 2, 2013, at 03:24 , Roslina Zakaria wrote: Dear r-users, I am using KS test to test the goodness of fit for my data and the got the following output. However, I don't understand about the warning messages. What does it mean by horizontals is not a graphical parameter It's

Re: [R] Kolmogorov-Smirnov Test

2013-08-02 Thread Roslina Zakaria
Hi Peter,   Thank you so much for your explaination.  I draw histogram before that, so maybe that warning messages are meant for that. From: peter dalgaard pda...@gmail.com Cc: r-help@r-project.org r-help@r-project.org Sent: Friday, August 2, 2013 3:11 PM

Re: [R] readBin into a data frame

2013-08-02 Thread Uwe Ligges
On 01.08.2013 10:36, Zhang Weiwu wrote: Hello. readBin is designed to read a batch of data with the same spec, e.g. read 1 floats into a vector. In practise I read into data frame, not vector. For each data frame, I need to read a integer and a float. for (i in 1:1000) {

[R] optim with(out) box constraints!

2013-08-02 Thread Anera Salucci
Hello  R users, Does optimizing a function using optim with  method= L-BFGS-B and without box constraints lead to L-BFGS optimization in R?  [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Correlation Loops in time series

2013-08-02 Thread TMiller
Thanks David, that went perfectly fine. Best Tom -- View this message in context: http://r.789695.n4.nabble.com/Correlation-Loops-in-time-series-tp4672732p4672879.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] optim with(out) box constraints!

2013-08-02 Thread Uwe Ligges
On 02.08.2013 10:10, Anera Salucci wrote: Hello R users, Does optimizing a function using optim with method= L-BFGS-B and without box constraints lead to L-BFGS optimization in R? Sort of, but the question is why this would be beneficial with today's computers ... Best, Uwe Ligges

Re: [R] skipping columns in read.table; was: (no subject)

2013-08-02 Thread Uwe Ligges
On 02.08.2013 05:29, Babu Guha wrote: I have a comma delimited file with 62 fields of which some are comments. There are about 1.5 million records/lines. Sme of the fields which has comments and which i do not need have 40 characters. Of the 62 fields, I will need at most 12 fields. What's

Re: [R] X11 and printing/copying graphics from mac

2013-08-02 Thread Prof Brian Ripley
On 02/08/2013 06:54, Marlin Keith Cox wrote: Hello all, I frequently use R, but recently have switched to a mac. When creating a new graphic, I have used X11() for a new graphics window and have always (when I used my PC) have been able to copy or print from this With which OS? graphic

Re: [R] readBin into a data frame

2013-08-02 Thread Jeff Newmiller
I find the hexView package quite helpful for this. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

[R] Write-table and dynamic access path

2013-08-02 Thread Robert U
Dear R-users, I am quite sure this is a beginner question, but i cannot manager to find the answer on the Internet... I am using write.table function inside some kind of loop. I'd like to write my tables in different folders without having to change the path inside the function every time,

Re: [R] Write-table and dynamic access path

2013-08-02 Thread Jim Holtman
?paste Sent from my iPad On Aug 2, 2013, at 6:20, Robert U tacsun...@yahoo.fr wrote: Dear R-users, I am quite sure this is a beginner question, but i cannot manager to find the answer on the Internet... I am using write.table function inside some kind of loop. I'd like to write my

Re: [R] Creating DAGS with plate notation in R

2013-08-02 Thread Raghu Naik
Thanks Sacha. The code works well although, as you alluded, the plates are being drawn after some positional fine tuning. I was hoping that the plates could come out automatically as they do in graphviz. Possibly, whenever you get to it, your package can alleviate that. Cheers. On Thu, Aug

Re: [R] Forest plot with sub-group analyses

2013-08-02 Thread Michael Dewey
At 18:19 01/08/2013, wanjiku gichohi wrote: Greetings. I am conducting a meta analyses that has Beta and SE values. I am wondering how to place the command to split the analyses. Could you assist me to develop the command within the syntax? I would also like to add the N and Year values in

Re: [R] optim with(out) box constraints!

2013-08-02 Thread Prof J C Nash (U30A)
It's possibly the L in L-BFGS-B that is more interesting for some problems, namely for the Limited Memory, so running without bounds can make sense. Unfortunately, the version of L-BFGS-B in R is from the 1990s, and Nocedal et al. released an update in 2011. Maybe someone will want to work on

Re: [R] Editing code in a function? related to allowing zero weights in lmer()

2013-08-02 Thread Mike Rennie
To close the thread, I wanted to post my solution to this issue. Searching some other help threads, I found this: https://stat.ethz.ch/pipermail/r-help/2008-August/171321.html So, here's what I did: I got at the script for lmerFrames by doing this: edit(lme4:::lmerFrames) then edited the

[R] verbose=1 setting in lmer gives no iteration results

2013-08-02 Thread jroyrobertson
I am using R version 3.0.1 (2013-05-16) on Platform: x86_64-w64-mingw32/x64 (64-bit) and lme4_0.9911-3. I run the following code from Douglas Bates lme4: Mixed-effects modeling with R (June 25, 2010) pg 14 and do not get any iteration information. I get the same results if I use

Re: [R] GO TO function in R

2013-08-02 Thread Barry Rowlingson
No. We figured this out in 1968. http://dl.acm.org/citation.cfm?doid=362929.362947 On Fri, Aug 2, 2013 at 3:20 AM, David Winsemius dwinsem...@comcast.net wrote: On Aug 1, 2013, at 6:58 PM, Peter Maclean wrote: Does R has a GO TO function? For example, I want to run this script at a given

[R] Regression Column names instead of numbers

2013-08-02 Thread TMiller
Hi guys I am new to R and I am currently trying to do a regression: I have two matrices with 200 time series each. In order to achieve a loop, I used the following command: sapply(1:200, function(x) summary(lm(formula=matrix1[,x]~matrix2[,x]))) Each column/time series has a unique name, in case

Re: [R] R runs on Windows Server 2012?

2013-08-02 Thread Jeff Newmiller
The FAQ may seem vague to you, but your question seems vague to me. I don't run R on any Windows Server, so I cannot answer your question, but there is more than one way to run R on any server and you have not mentioned how you want it to run. For example, you might run it via Remote Terminal

Re: [R] Write-table and dynamic access path

2013-08-02 Thread Jeff Newmiller
Or ?file.path for a more OS- independent approach. --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] Forest plot with sub-group analyses

2013-08-02 Thread David Winsemius
On Aug 2, 2013, at 5:20 AM, Michael Dewey wrote: At 18:19 01/08/2013, wanjiku gichohi wrote: Greetings. I am conducting a meta analyses that has Beta and SE values. I am wondering how to place the command to split the analyses. Could you assist me to develop the command within the

Re: [R] Regression Column names instead of numbers

2013-08-02 Thread David Carlson
It depends on how fancy you want to get. A quick fix would be pairs - paste0(colnames(matrix1), ., colnames(matrix2)) # lapply will be faster since you are returning a list results - lapply(1:200, function(x) summary(lm(formula=matrix1[,x]~matrix2[,x]))) names(results) - pairs results The

Re: [R] rms plot.Predict when type=p

2013-08-02 Thread Frank Harrell
Mike, 2 out of 3 of your requests are handled simply. I've tidied up your code and show how to include extra parameters into the plot( ) call below. The request to change plotting symbol color for each factor level is not consistent with the lattice philosophy for non-superposed variables

[R] lme: object is not a matrix

2013-08-02 Thread Puschner, Bernd
Dear all, when running in R a simple lme call (3 variables in data set) which works perfectly fine in S-PLUS, I keep getting the error message Error in model.frame.default(formula = ~t + can2p + code, data = list( : object is not a matrix Any help would be very much appreciated.

[R] cannot base64decode string which is base64encode in R

2013-08-02 Thread Qiang Wang
Hi, I'm struggling with encode/decode strings in R. Don't know why the second example below would fail. Thanks in advance for your help. succeed: s - saf x - base64encode(s) y - base64decode(x, character) fail: s - safs x - base64encode(s) y - base64decode(x, character) [[alternative

[R] Forest plot with sub group analyses

2013-08-02 Thread wanjiku gichohi
Greetings. I am conducting a meta analyses that has Beta and SE values. I am wondering how to place the command to split the analyses. Could you assist me to develop the command within the syntax? I would also like to add the N and Year values in the plot. I would like my plot to look like the

[R] HiPLARb installation failed when magma-lib included

2013-08-02 Thread Clara Anton Fernandez
Hi all, I am trying to install HiPLARb package but I am running into an error. Any help would be greatly appreciated. I used the installer available from their website (http://www.hiplar.org/hiplar-b-installation.html) and followed their instructions. The installer finished without showing

[R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread Klot Lee
Hi, The function locator() allows me to add points by mouse on function plot() I have draft. But the points can not be edit once I have made the click. How can I remove or just move the points I made on the draft? I use the locator function like this: plot(1,1) locator(type=o) Thanks.

[R] Creating a model with fixed and random variables

2013-08-02 Thread cbarl3y
I am trying to recreate a model that would mimic results a peer obtained using SAS, but in R. The goal of this model would be to determine if there is any variation among parental sets and lines (essentially is there a significant p-value for the variables “set” and “line”). The problems I am

Re: [R] Forest plot with sub group analyses

2013-08-02 Thread John Kane
No attachments came through. Either include the code and the data in the email---see ?dupt for a way to include the data or include the files as .txt files. John Kane Kingston ON Canada -Original Message- From: wngich...@gmail.com Sent: Fri, 2 Aug 2013 10:45:33 +0200 To:

Re: [R] Regression Column names instead of numbers

2013-08-02 Thread TMiller
Hello David Thanks for your answer. It works with the number in the double bracket above each of the regression results. However, the x's still remain in the call formula. I really appreciate any help Best Tom -- View this message in context:

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread David Winsemius
On Aug 2, 2013, at 2:26 AM, Klot Lee wrote: Hi, The function locator() allows me to add points by mouse on function plot() I have draft. But the points can not be edit once I have made the click. How can I remove or just move the points I made on the draft? You cannot. The plotting model

Re: [R] Regression Column names instead of numbers

2013-08-02 Thread arun
You could try: set.seed(25) mt1- matrix(sample(c(NA,1:40),20*200,replace=TRUE),ncol=200) colnames(mt1)- paste0(X,1:200) set.seed(487) mt2- matrix(sample(c(NA,1:80),20*200,replace=TRUE),ncol=200) colnames(mt2)- colnames(mt1) res-lapply(colnames(mt1),function(x) {x1-data.frame(mt1[,x],mt2[,x]);

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread Daniel Guetta
Depending on the complexity of your graphic, you could also just re-draw the entire thing every time you want to make a change. Daniel On Fri, Aug 2, 2013 at 10:17 AM, David Winsemius dwinsem...@comcast.netwrote: On Aug 2, 2013, at 2:26 AM, Klot Lee wrote: Hi, The function locator()

Re: [R] Removing null vector that is not in .csv file data

2013-08-02 Thread arun
Hi, You could remove it by: vec1- c( ,1003PI_pilot0101, 1003PI_pilot0102,1003PI_pilot0103, 1003PI_pilot0104, 1003PI_pilot0105, 1003PI_pilot0107) vec1[vec1!=] #[1] 1003PI_pilot0101 1003PI_pilot0102 1003PI_pilot0103 1003PI_pilot0104 #[5] 1003PI_pilot0105 1003PI_pilot0107 #or  vec1[nchar(vec1)0]

Re: [R] lme: object is not a matrix

2013-08-02 Thread Bert Gunter
R is not splus. Read ?lme. Bert Sent from my iPhone -- please excuse typos. On Aug 2, 2013, at 5:45 AM, Puschner, Bernd bernd.pusch...@bkh-guenzburg.de wrote: Dear all, when running in R a simple lme call (3 variables in data set) which works perfectly fine in S-PLUS, I keep getting

Re: [R] R runs on Windows Server 2012?

2013-08-02 Thread Chris Solomon
Hi Jeff- Thanks for your input. I will try to provide a little more detail. First, it is useful to know that I know very little about high-performance computing with R or in general. We occasionally run models that take a lot of processor time. In the past I have used quad-core desktops and

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread David Carlson
You cannot move or edit them interactively, but you can save the locations of those points by using: a - locator(type=o) # e.g. click on three locations a $x [1] 0.8559376 1.1778126 1.2340626 $y [1] 1.2256976 1.2501162 0.8402324 You could now delete or modify the points and then use lines()

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread Greg Snow
The function put.points.demo in the TeachingDemos package allows you to add points, delete points, and move points on a scatterplot (and also computes some summaries based on the points). It actually uses locator behind the scenes and redraws the plot each time, but it gives the effect. If that

[R] CSPADE error: system invocation error

2013-08-02 Thread Cen Guo
Hi Sue Michael, I'm running into the same situation, i first tried everything with my own data and it didn't work, and then i tried with zaki.txt and i got the same error Here is my sessioninfo: sessionInfo() R version 3.0.1 (2013-05-16) Platform: i386-w64-mingw32/i386 (32-bit) locale: [1]

[R] R Job taking longer in Linux than in Windows

2013-08-02 Thread alina andrei
Hello,   I am running an R job on a Windows 7 machine, having 4 cores and 16GB RAM , R 3.0.1, and it takes 1.5 hours to complete. I am running the same job in R on a Linux enviroment (Platform: x86_64-redhat-linux-gnu (64-bit)) with huge amounts of memory: 40 cores and .5 TB RAM., and the job

Re: [R] problem about mean function in ffbase package

2013-08-02 Thread Milan Bouchet-Valat
Le jeudi 01 août 2013 à 00:10 +0800, Chaos Chen a écrit : Hi all, I experienced some unmatched result using mean function in ffbase package and cannot figure out what's wrong. I have a simulated ff vector with 10 numbers inside and want to calculate its mean. But the results are

[R] How to index each loop results

2013-08-02 Thread Dadong Zhang
*Hi R Users, * *I was trying to write a simple function for data imputation. I need to index the results from each step and combine it with some data of the original dataset for the next step, until the end. I attached a sample data and the work-well R code function I wrote, but the function is

[R] Error executing R CMD check --as-cran caused by missing Maintainer field

2013-08-02 Thread Stefano Costa
Hello, while checking a small package I wrote (http://cran.r-project.org/web/packages/tipom/) before submitting a new version, I was blocked by an error in the checking CRAN incoming feasibility section of R CMD check --as-cran, as pasted below: steko@ganymede:~/code$ LC_ALL=C R CMD check

Re: [R] problem about mean function in ffbase package

2013-08-02 Thread Milan Bouchet-Valat
Le jeudi 01 août 2013 à 00:10 +0800, Chaos Chen a écrit : Hi all, I experienced some unmatched result using mean function in ffbase package and cannot figure out what's wrong. I have a simulated ff vector with 10 numbers inside and want to calculate its mean. But the results are

Re: [R] lme: object is not a matrix

2013-08-02 Thread zelfortin
Hi, make sure your data is in a data frame: Data - data.frame(t,can2p,code) just specify in which object you find those specific variables by using object$ before each one. I hope this solve the problem. If not, post your script, it's really hard to work with only an error message! Cheers

Re: [R] Removing null vector that is not in .csv file data

2013-08-02 Thread zelfortin
Hi, you could try going in your spreadsheet and delete columns/rows that are seemingly empty (delete a lot a of them), then try again. Sometimes in your spreadsheet you enter values then you delete them. This leaves the cells active and when you save the csv file, it writes some commas with

Re: [R] R Job taking longer in Linux than in Windows

2013-08-02 Thread Uwe Ligges
On 02.08.2013 22:14, alina andrei wrote: Hello, I am running an R job on a Windows 7 machine, having 4 cores and 16GB RAM , R 3.0.1, and it takes 1.5 hours to complete. I am running the same job in R on a Linux enviroment (Platform: x86_64-redhat-linux-gnu (64-bit)) with huge amounts of

[R] How to download this data?

2013-08-02 Thread Ron Michael
Hi all,   I need to download the data from this web page:   https://www.theice.com/productguide/ProductSpec.shtml?specId=219#expiry   I used the function readHTMLTable() from package XML, however could not download that.   Can somebody help me how to get the data onto my R window?   Thank you.

Re: [R] using rollapply to calculate a moving sum or running sum?

2013-08-02 Thread Anika Masters
This is not critical, but I am curious to learn. Are there any suggestions for speeding up the process to calculate a moving row sum? (Ideally from within R, as opposed to suing C, etc.) Using rollapply on a matrix of 45,000 rows and 400 columns takes 83 minutes. date() mymatrix -

Re: [R] using rollapply to calculate a moving sum or running sum?

2013-08-02 Thread William Dunlap
Try the filter() function. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Anika Masters Sent: Friday, August 02, 2013 2:48 PM To: arun Cc: R help Subject: Re: [R]

Re: [R] How to download this data?

2013-08-02 Thread Duncan Temple Lang
That URL is an HTTPS (secure HTTP), not an HTTP. The XML parser cannot retrieve the file. Instead, use the RCurl package to get the file. However, it is more complicated than that. If you look at source of the HTML page in a browser, you'll see a jsessionid and that is a session identifier. The

Re: [R] (no subject)

2013-08-02 Thread Jim Lemon
On 08/02/2013 01:29 PM, Babu Guha wrote: I have a comma delimited file with 62 fields of which some are comments. There are about 1.5 million records/lines. Sme of the fields which has comments and which i do not need have 40 characters. Of the 62 fields, I will need at most 12 fields. What's

[R] access the tcl history stream in tcltk

2013-08-02 Thread Christopher Marcum
Hello, Does anyone have a clue about how dump tcl's history stream out to a file from R pointers to tcl? This should be accessible via tcl(history,...) but the stream seems to be always empty. library(tcltk) #Load tcltk tcl(history,nextid) #should be 1, and it is a1-tclVar(1) tclvalue(a1)-2

[R] Internalization of help pages

2013-08-02 Thread Tomáš Greif
Is it possible to translate help files? I see there are some localization options for GUI, but not for help. I think this would be really helpful for users who don't have English as their primary language. Because all the online help is created from text files, it should be not that difficult to

Re: [R] Internalization of help pages

2013-08-02 Thread Bert Gunter
Tomáš 1. Your concern is valid, especially now that R has become so widely used throughout the world. But language translation is hard. 2. Warning and error messages are short, and relatively trivial to translate. 3. There are tens of thousands of Help pages in thousands of packages, most

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread Klot Lee
Maybe I have to white a function for the purpose base on loacator. Thanks for the reply. 2013/8/3 David Carlson dcarl...@tamu.edu You cannot move or edit them interactively, but you can save the locations of those points by using: a - locator(type=o) # e.g. click on three locations a

[R] Singular fit error with R rlm function

2013-08-02 Thread Krishnanand Khambadkone
Has anyone used the rlm regression function. I am using a factor variable as part of this regression. If the number of occurrences of any factor is 0 or 1 it complains with the error that 'x' is singular: singular fits are not implemented in 'rlm' The droplevels function removes 0 occurrence

Re: [R] could I remove or move points in locator() based on plot()?

2013-08-02 Thread Klot Lee
Thanks. I have tried put.points.demo, it helps, maybe some more modify then it would be suitable for me. Thank you again. 2013/8/3 Greg Snow 538...@gmail.com The function put.points.demo in the TeachingDemos package allows you to add points, delete points, and move points on a scatterplot