Re: [R] Converting factor data into Date-time format

2014-09-30 Thread PIKAL Petr
Hi Everything seems to be OK. Please check your if your objects are as required by. str(objects) (and post result if you are in doubt :) from your attempts I presume that datetime was not created properly Maybe this can do it, without actual data it is hard to tell exactly. datetime=as.POSIXct(

Re: [R] Inverse Student t-value

2014-09-30 Thread Andre
Hi JLucke, Maybe the old excel function. TINV and T.INV.2T is same function for Two-Tailed Inverse of the student`s t-distribution but T.INV use for Left-Tailed inverse of the Student's t-distribution and can be use for Inverse of the student`s t-distribution. I know automatic or functions any

[R] uninstalled and reinstalled R, now cannot install packages

2014-09-30 Thread Sarah Hardy
Hello, I have a student (running Windows 8) who ran into a problem installing the RcmdrPlugin.IPSUR package and, despite instructions to the contrary, uninstalled R and reinstalled. He uninstalled it using the Windows 8 utility. Now he cannot install any packages. Here are the error messages he

[R] Using R for modelling Australian Senate Election in NSW?

2014-09-30 Thread Philip Rhoades
People, I am setting up an Australian Science and Technology party with Human Health and Longevity as it's defining, first platform plank: http://lestp.org and intend to run Senate candidates at the next Federal Election (due late 2016). Senators are elected from each State and Territory

Re: [R] Best Distribution

2014-09-30 Thread eliza botto
Dear Turner and Murphy, Thankyou so very much for replying. It fixed the issue. :) Eliza > Date: Wed, 1 Oct 2014 11:04:32 +1300 > From: r.tur...@auckland.ac.nz > To: eliza_bo...@hotmail.com; r-help@r-project.org > Subject: Re: [R] Best Distribution > > On 01/10/14 08:50, eliza botto wrote: > > D

Re: [R] Reading text file with fortran format

2014-09-30 Thread Steven Yen
Thanks to all. Steven Yen At 06:18 PM 9/30/2014, Nordlund, Dan (DSHS/RDA) wrote: > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Steven Yen > Sent: Tuesday, September 30, 2014 2:04 PM > To: r-help > Subject: [R] Reading text

Re: [R] Reading text file with fortran format

2014-09-30 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Steven Yen > Sent: Tuesday, September 30, 2014 2:04 PM > To: r-help > Subject: [R] Reading text file with fortran format > > Hello > > I read data with fortran format: > mydata<

Re: [R] Best Distribution

2014-09-30 Thread Rolf Turner
On 01/10/14 08:50, eliza botto wrote: Dear useRs, I have this following data dput(Prec) c(42.2, 45.2, 46, 48, 54, 54.1, 59.4, 61, 62.2, 63.5, 65.024, 71.9, 73.4, 76.6, 76.708, 77.5, 77.724, 78, 81.3, 84.7, 84.836, 85.09, 88.2, 91.4, 94, 95.8, 96, 97.3, 101, 101, 101.5, 102.3, 102.87, 108.7, 1

Re: [R] Regarding R package spatstat

2014-09-30 Thread Rolf Turner
Why are you trying to install version 1.26-0? That is very ancient; it dates from April 2012. The version currently on CRAN is 1.38-1. I would suggest that you try that. I have no idea why your installation is hanging at the byte-compile stage, nor any idea how to diagnose the problem. I

[R] Reading text file with fortran format

2014-09-30 Thread Steven Yen
Hello I read data with fortran format: mydata<-read.fortran('foo.txt', c("4F10.4","F8.3","3F3.0","20F2.0")) colnames(mydata)<-c("q1","q2","q3","q4","income","hhsize", "weekend","dietk","quart1","quart2","quart3","male","age35", "age50","age65","midwest","south","west","no

Re: [R] Converting factor data into Date-time format

2014-09-30 Thread David L Carlson
First, use stringsAsFactors=FALSE with the read.csv() function. That will prevent the conversion to factors. Then try to convert date and time to datetime objects. - David L Carlson Department of Anthropology Texas A&M University College Station, TX 77840-435

[R] Best Distribution

2014-09-30 Thread eliza botto
Dear useRs, I have this following data > dput(Prec) c(42.2, 45.2, 46, 48, 54, 54.1, 59.4, 61, 62.2, 63.5, 65.024, 71.9, 73.4, 76.6, 76.708, 77.5, 77.724, 78, 81.3, 84.7, 84.836, 85.09, 88.2, 91.4, 94, 95.8, 96, 97.3, 101, 101, 101.5, 102.3, 102.87, 108.7, 109.5, 110.5, 110.7, 112, 114.3, 118.11,

Re: [R] Inverse Student t-value

2014-09-30 Thread JLucke
The website has your answer. The t-distribution is a regularized incomplete beta function. The incomplete beta function is given by R's pbeta function. You regularize it with R's beta function. Then you use R's uniroot function to find the inverse. Good homework problem. Andre Sent by:

[R] Converting factor data into Date-time format

2014-09-30 Thread tandi perkins
Hello R help: I am new to this forum so I apologize in advance for any protocol missteps. I have a data set that is comprised of eight birds with GPS; each of which transmit everyday at 8:00 am, 4:00 pm, and midnight for 1 year (although I have some missing relocation's). I am trying to forma

Re: [R] Inverse Student t-value

2014-09-30 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of jlu...@ria.buffalo.edu > Sent: Tuesday, September 30, 2014 11:49 AM > To: Andre > Cc: r-help@r-project.org > Subject: Re: [R] Inverse Student t-value > > My Excel (2013) returns

Re: [R] Inverse Student t-value

2014-09-30 Thread Andre
Hi Duncan, Actually, I am trying trace the formula for the "Critical value of Z" and manual formula is =(I7-1)/SQRT(I7)*SQRT((TINV(0.05/I7,I7-2))^2/(I7-2+TINV(0.05/I7,I7-2))) So, I got new problem for TINV formula. I just need a manual equation for TINV. Hope solve this problem. Cheers! On We

Re: [R] Inverse Student t-value

2014-09-30 Thread Andre
Hi Duncan, Let me explain again, I just need a manual expression for inverse student t value. You could go to web page http://www.danielsoper.com/statcalc3/calc.aspx?id=10 That's inverse student t value calculator. Do you know a manual expression use it. Cheers! On Wednesday, October 1, 2014,

Re: [R] Inverse Student t-value

2014-09-30 Thread Andre
Hi Duncan, No, that's correct. Actually, I have data set below; N= 1223 alpha= 0.05 Then probability= 0.05/1223=0.408831 degree of freedom= 1223-2= 1221 So, TINV(0.408831,1221) returns 4.0891672 Could you show me more detail a manual equation. I really appreciate it if you may give mo

Re: [R] Inverse Student t-value

2014-09-30 Thread peter dalgaard
On 30 Sep 2014, at 20:36 , Duncan Murdoch wrote: >> TINV(0.408831,1221) Just for the record: The above _does_ give 4.117456652 in Excel for Mac OS X. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (

Re: [R] Inverse Student t-value

2014-09-30 Thread Marc Schwartz
FWIW, I get: 4.117456652 in Excel 2011 on OS X and: 4.117457 in R 3.1.1 on OS X There is a KB article on the TINV function here, suggesting that the threshold for the iterative algorithm in Excel has been tightened in recent versions: http://support.microsoft.com/kb/828340 Regards, Ma

Re: [R] Inverse Student t-value

2014-09-30 Thread JLucke
My Excel (2013) returns exactly what R does. I used both T.INV and T.INV.T2There is no TINV. Has Excel been updated? Duncan Murdoch Sent by: r-help-boun...@r-project.org 09/30/2014 02:36 PM To Andre , cc r-help@r-project.org Subject Re: [R] Inverse Student t-value On 30/09/20

Re: [R] Inverse Student t-value

2014-09-30 Thread Duncan Murdoch
On 30/09/2014 2:26 PM, Andre wrote: Hi Duncan, Actually, I am trying trace the formula for the "Critical value of Z" and manual formula is =(I7-1)/SQRT(I7)*SQRT((TINV(0.05/I7,I7-2))^2/(I7-2+TINV(0.05/I7,I7-2))) So, I got new problem for TINV formula. I just need a manual equation for TINV.

Re: [R] Inverse Student t-value

2014-09-30 Thread Ted Harding
And, with 1221 degrees of freedom, one cannot be far off a Normal distribution. So: abs(qnorm(0.408831/2)) [1] 4.102431 which is nearer to Duncan's answer (4.117, and a bit smaller, as it should be) than to yours (4.0891672, which, if accurate, should be greater than the Normal value 4.10

Re: [R] Inverse Student t-value

2014-09-30 Thread Duncan Murdoch
On 30/09/2014 2:11 PM, Andre wrote: Hi Duncan, No, that's correct. Actually, I have data set below; Then it seems Excel is worse than I would have expected. I confirmed R's value in two other pieces of software, OpenOffice and some software I wrote a long time ago based on an algorithm publ

Re: [R] Inverse Student t-value

2014-09-30 Thread Duncan Murdoch
On 30/09/2014 1:31 PM, Andre wrote: Dear Sir/Madam, I am trying to use calculation for two-tailed inverse of the student`s t-distribution function presented by Excel functions like =TINV(probability, deg_freedom). For instance: The Excel function =TINV(0.408831,1221) = returns 4.0891672.

[R] Inverse Student t-value

2014-09-30 Thread Andre
Dear Sir/Madam, I am trying to use calculation for two-tailed inverse of the student`s t-distribution function presented by Excel functions like =TINV(probability, deg_freedom). For instance: The Excel function =TINV(0.408831,1221) = returns 4.0891672. Would you like to show me a manual ca

Re: [R] Issue with a function

2014-09-30 Thread Yara Abu Awad
Yes, thank you. I was able to resolve the issue using this code: as.matrix(newX%*%Diagonal(length(beta), x=beta)) On Tuesday, September 30, 2014 3:38 AM, peter dalgaard wrote: On 30 Sep 2014, at 03:47 , Jeff Newmiller wrote: > library(Matrix) > result <- newX %*% Diagonal( unlist(beta) ) T

[R] Course Halifax: Introduction to Linear mixed effects models, GLMM and MCMC with R

2014-09-30 Thread Highland Statistics Ltd
Apologies for cross-posting We would like to announce the following statistics course: Course: Introduction to Linear mixed effects models, GLMM and MCMC with R Location: Halifax, Canada Date: 19 - 23 January 2015 Remaining seats: 10 Course website: http://www.highstat.com/statscourse.

Re: [R] modify function in a package

2014-09-30 Thread C Lin
David, thank you for pointing me in the right direction! For the record, I was able to override the NanoStringNorm function while still calling other functions in the package using the following: source('my.nanostringnorm.R") # inside my.nanostringnorm.R, my modified function is called NanoStrin

[R] Using "survey" package with ACS PUMS

2014-09-30 Thread Michael . Laviolette
I'm trying to reproduce some results from the American Community Survey PUMS data using the "survey" package. I'm using the one-year 2012 estimates for New Hampshire (http://www2.census.gov/acs2012_1yr/pums/csv_pnh.zip) and comparing to the estimates for user verification from http://www.census.go

Re: [R] Loop does not work: Error in else statement (II)

2014-09-30 Thread PIKAL Petr
Hi Slightly better but still html scrambled. see in line > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Frank S. > Sent: Tuesday, September 30, 2014 2:55 PM > To: r-help@r-project.org > Subject: [R] Loop does not work: Erro

[R] Regarding R package spatstat

2014-09-30 Thread shivali gangwar
Hi , I am installing a R Package as following - R CMD INSTALL spatstat_1.26-0.tar.gz It seems my installation hangs at "byte-compile and prepare package for lazy loading" I kept on waiting for 1/2anhr but it hanged at below stage. Please help. ~~~ Target "all" i

Re: [R] Using "survey" package with ACS PUMS

2014-09-30 Thread Anthony Damico
hi michael, you probably need options( "survey.replicates.mse" = TRUE ) i also think you don't want type = "Fay" and you do want scale = 4/80 and rscales = rep( 1 , 80 ) as well, but what you have might be equivalent (not sure) regardless, this blog post details how to precisely replicate the

Re: [R] Using "survey" package with ACS PUMS

2014-09-30 Thread Michael . Laviolette
Setting the options did the trick. The code using "Fay" came from another post, but this works: options( "survey.replicates.mse" = TRUE) pums_p.rep <- svrepdesign(repweights = pums_p[207:286], weights = ~PWGTP, combined.weights = TRUE,

Re: [R] [rJava] RJavaClassLoader and system classloader

2014-09-30 Thread Benoît Thiébault
Thank you very much for your answer. I will look further into R Service Bus. In the meantime, I have posted an issue on GitHub https://github.com/s-u/rJava/issues/32 where I describe a way to reproduce the unexpected behaviour. Kind regards, Ben Le 30 sept. 2014 � 12:23, Mohan Radhakrishnan

Re: [R] Loop does not work: Error in else statement

2014-09-30 Thread Frank S.
Dear Berend and Petr, I do apologise for the disorderly code I posted. I have tried to solve it in a new mail. Frank S. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list htt

[R] Loop does not work: Error in else statement (II)

2014-09-30 Thread Frank S.
Hi to all members of R list, I�m working with data.table package, and with 6 variables: ID: Subject identifier born: Birthdate start: Starting date register: date of measurement value: Value of measurement end: date of expiration of the measurements. So, the natural order of d

Re: [R] ..nlme: start values necessary even though using self-start function?..

2014-09-30 Thread Prof Brian Ripley
This seems to be a bug in the package. Maybe the code was written for S, or maybe R has changed in the many years since the former maintainer decided to be unhelpful. I'll put this on the list for the next release. On 30/09/2014 09:13, Matthias Kuhn wrote: Dear list. I am stuck with an "su

Re: [R] Obtain 00Index.html

2014-09-30 Thread Duncan Murdoch
On 30/09/2014, 5:41 AM, Sven E. Templer wrote: > Hello, > > how can I open (by an R command) the index page in html mode, as obtained by: > > options(browser="firefox") # or any other > options(help_type="html") > ?help > # and then following the html reference on the page bottom named "Index" >

[R] coplot: extract regression output

2014-09-30 Thread Marco Helbich
Dear list, I would like to extract the model details (coefficients, p-values etc.) for the regression lines within a coplot. How are these results accessible? A code example is as follows: dat <- data.frame(a=rnorm(111), b=rnorm(111), c=rnorm(111)) coplot(dat$a ~ dat$b | dat$c, data=dat, cex=

Re: [R] [rJava] RJavaClassLoader and system classloader

2014-09-30 Thread Mohan Radhakrishnan
Hi, You could ask the author or post an issue here .( https://github.com/s-u/rJava) I had faced a problem and I got a response from him. In my case I am streaming JVM data and this( http://www.openanalytics.eu/r-service-bus) may be a better idea. I haven't still tried this. Thanks, Mohan

[R] Obtain 00Index.html

2014-09-30 Thread Sven E. Templer
Hello, how can I open (by an R command) the index page in html mode, as obtained by: options(browser="firefox") # or any other options(help_type="html") ?help # and then following the html reference on the page bottom named "Index" In text mode I know library(help='utils') to open the utils pack

[R] ..nlme: start values necessary even though using self-start function?..

2014-09-30 Thread Matthias Kuhn
Dear list. I am stuck with an "subscript out of bounds" error when I play with an nlme-example from the yellow Pinheiro/Bates book. It is about using nlme() using the formula interface for non-linear mixed models. The following code snippet is inspired by the book (in section 8.2), using the Oran

Re: [R] Median of streaming data

2014-09-30 Thread Mohan Radhakrishnan
Hi, I came across this project(http://jwijffels.github.io/RMOA/) which seems to be directly addressing streaming prediction. >From the manual : A data stream environment has different requirements from the traditional setting. The most significant are the following: Requirement 1 Process a

Re: [R] Issue with a function

2014-09-30 Thread peter dalgaard
On 30 Sep 2014, at 03:47 , Jeff Newmiller wrote: > library(Matrix) > result <- newX %*% Diagonal( unlist(beta) ) The code in Matrix does this efficiently, I hope? Otherwise, sweep() is the traditional ticket: result <- sweep(newX, 2, unlist(beta), "*") Some of my students do t(t(newX) * unl