[R] poLCA problem

2021-03-05 Thread Scott Colwell
e variables but shortened for this purpose df2 <- df[items] i <- cbind(x1, x2)~1 poLCA (i, df2, nclass=2, maxiter=100, nrep=10, verbose =TRUE) It is after the poLCA that I get the error. Any thoughts on what is causing this? Thanks, Scott

[R] cooks distance for repeated measures anova

2020-12-23 Thread Walker Scott Pedersen
Hi all, Is there a way to get cook's distance for a repeated measures anova? Neither cooks.distance or CookD from the predictmeans package seem to allow for this. For example, if I have the model data(iris) mod<-aov(Sepal.Length ~ Petal.Length + Petal.Width + Error(Species), data=iris) both

[R] Quantitative Methods Workshops in May 2020

2020-04-06 Thread Scott Colwell
in both R (using R-Studio) and Mplus. Syntax and output for both programs will be provided for all examples covered in the workshop. Thank you, Scott [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To

[R] choose(n, k) as n approaches k

2020-01-14 Thread Wright, Erik Scott
This struck me as incorrect: > choose(3.99, 4) [1] 0.979 > choose(3.999, 4) [1] 0 > choose(4, 4) [1] 1 > choose(4.001, 4) [1] 4 > choose(4.01, 4) [1] 1.02 Should base::choose(n, k) check whether n is within machine precision of k and return 1? Thanks, Erik *** sessionIn

Re: [R] Creating hanging bar plot in r from dplyr

2019-04-20 Thread Scott Colwell
That is perfect. Thanks! -- Scott R. Colwell, PhD On 2019-04-20, 1:23 PM, "Jeff Newmiller" wrote: Not really sure I understand what you want. Here is some code to consider: library(ggplot2) library(dplyr) library(tidyr)

[R] Creating hanging bar plot in r from dplyr

2019-04-20 Thread Scott Colwell
ly appreciated. Thank you! Scott -- Scott R. Colwell, PhD [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guid

[R] 2019 Spring Workshops using R

2019-04-13 Thread Scott Colwell
without covariates, growth modeling with time invariant and varying covariates, centering points, piecewise growth modeling, autoregressive latent trajectory modeling (ALT -- Scott R. Colwell, PhD, CStat, PStat [[alternative HTML version deleted

[R] plot one levelplot over another

2018-11-19 Thread Waichler, Scott R
surface show through. How can I do this? Regards, Scott Waichler Pacific Northwest National Laboratory Richland, WA USA __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] GLM Model Summary

2018-10-16 Thread Neslin, Scott A.
an Excel table easily rather than having to count the number of characters. Is this possible for us to do by ourselves? Or could you modify the package in some way? We appreciate your attention. Thank you! Scott Neslin Prasad Vana Dartmouth College [[alternative HTML version deleted

[R] seq() problem with chron

2018-09-05 Thread Waichler, Scott R
chron) dt1 <- chron("02/20/13", "00:00:00") dt2 <- chron("07/03/18", "15:30:00") dt <- seq(from=dt1, to=dt2, by=1/(24*4)) dt[length(dt)] #[1] (07/03/18 15:15:00) Thanks, Scott Waichler Pacific Northwest National Laboratory scott.waich...@pnnl.g

[R] download.file() problems with binary files containing EOF byte in Windows

2018-08-20 Thread Scott Sherrill-Mix
;) download.file('file://gunzip.c.gz','dl.gz',mode='wb',method='internal') file.size('dl.gz') It looks like the 'auto' and 'internal' methods both resolve to the 'wininet' method on Windows and mode is automatically

Re: [R] [Rd] source(echo = TRUE) with a iso-8859-1 encoded file gives an error

2018-05-05 Thread Scott Kostyshak
On Fri, May 04, 2018 at 10:58:26PM +, Ista Zahn wrote: > On Fri, May 4, 2018 at 4:47 PM, Scott Kostyshak wrote: > > I have very little knowledge about file encodings and would like to > > learn more. > > > > I've read the following pages to

[R] [Rd] source(echo = TRUE) with a iso-8859-1 encoded file gives an error

2018-05-04 Thread Scott Kostyshak
_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] compiler_3.4.3 > Sys.getlocale() [1] "LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en

[R] Newbie - Scrape Data From PDFs?

2018-01-23 Thread Scott Clausen
Hello, I’m new to R and am using it with RStudio to learn the language. I’m doing so as I have quite a lot of traffic data I would like to explore. My problem is that all the data is located on a number of PDFs. Can someone point me to info on gathering data from other sources? I’ve been to the

[R] convertTime package.

2017-10-31 Thread Scott Anderwald via R-help
, Scott Anderwald __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

[R] ggplot's aes_ doesn't work as expected for x=factor()

2017-08-13 Thread Waichler, Scott R
the x variable a factor ahead of time df2 <- df df2$Var1 <- as.factor(df2$Var1) ggplot(df2, aes_(x=as.name(varnames[2]), y=as.name(varnames[3]), fill=as.name(varnames[1]))) + geom_bar(stat="identity") # works Regards, Scott Waichler Pacific Northwest National Laboratory Richland, WA,

[R] plot3D color ramp not working as expected

2017-06-29 Thread Waichler, Scott R
colvar, theta = 60, col = c("blue", "green", "red"), breaks = c(-0.5, -0.1, 0.1, 0.5)) Please help. Thanks, Scott Scott Waichler, PhD Pacific Northwest National Laboratory scott.waich...@pnnl.gov Richland, Washington, USA __

Re: [R] customizing color key with plot3D

2017-06-21 Thread Waichler, Scott R
Previously I had to use a completely different software, VisIt, to do these kinds of composite plots. R is my preferred tool and it is great to finally be able to do these at "home". Thanks, Scott Waichler scott.waich...@pnnl.gov Pacific Northwest National Laboratory Richland, Washin

[R] customizing color key with plot3D

2017-06-20 Thread Waichler, Scott R
cex.axis = 0.8, cex.clab = 0.85, addlines=TRUE, tick=FALSE, at = 1:3, labels=paste("Class", 1:3)), clab = c("","Elev Classes"), plot = TRUE) Thanks, Scott Scott Waichler Pacific Northwest National Laboratory Richland, Washington, USA

[R] how apply.monthly() in package xts works

2017-03-09 Thread Waichler, Scott R
x, df2$dates) #[,1] # 1992-06-182 # 1993-06-224 # 1994-06-075 # 1995-06-086 Two questions: 1) Is there a more elegant way to do this? 2) Shouldn't the xts documentation discuss the problem of sparse data? Regards, Scott Waichler Pacific Northwest National Laborator

[R] odd behavior of numeric()

2016-04-27 Thread Scott Rifkin
Why does: > numeric(0.2*25) return [1] 0 0 0 0 0 but > numeric((1-0.8)*25) returns [1] 0 0 0 0 [running version 3.2.0] [Apologies if this has been asked before - it's a hard question to find specific search terms for] Thanks, Scott _

[R] R packages for Mac Users

2016-02-23 Thread Scott Colwell
Hello, Does anyone know if all the R packages that are available for Windows users are also available for Mac users? Thank you, Scott -- Scott R. Colwell, Ph.D. Associate Professor, Dept. of Mkt/Cons Studies Adjunct Professor, Dept. of Psychology University of Guelph Guelph, Ontario

Re: [R] Output In R

2015-08-18 Thread Scott Tetrick
?write.csv and look at with the editor of choice. On Tue, Aug 18, 2015 at 6:41 AM, Shivi82 wrote: > Hello All, > > As i am a newbie in R so most of you would have seen this question zillion > times. I searched for the answer on this forum as well on other various > forums however could not fi

[R] How to 2d cross sections from a 3d finite element mesh

2015-08-10 Thread Waichler, Scott R
interpolation using Delaunay triangulation, but I wonder if there is a package out there to simplify things. I don't know at all how to go about doing it the first way. Can anyone suggest or point me to existing methods? Thanks, Scott Waichler Pacific Northwest National Laboratory Richland, WA

Re: [R] Can't install rgl: installed package can't be loaded; 'memory not mapped'

2015-08-08 Thread Waichler, Scott R
d_square GL_NV_texture_rectan! gle GL_NV_texgen_reflection GL_SGIS_generate_mipmap GL_SGIS_texture_lod GL_SGIX_depth_texture GL_SGIX_shadow 24839 frames in 6.0 seconds = 4153.389 FPS 7152 frames in 6.0 seconds = 1192.594 FPS . . . Scott > > ** testing if installed package can be loaded > > sh: li

[R] Can't install rgl: installed package can't be loaded; 'memory not mapped'

2015-08-08 Thread Waichler, Scott R
installing to /files3/R/R-3.2.1_install/rgl.Rcheck/rgl/libs ** R ** demo ** inst ** preparing package for lazy loading ** help *** installing help indices ** building package indices ** installing vignettes ** testing if installed package can be loaded sh: line 1: 11949 Segmentation fault '

[R] write.table with append=T after using cat on same file

2015-07-27 Thread Waichler, Scott R
uot;junk.txt", sep="", "# An introductory note.\n") write.table(df, file="junk.txt", sep=",", append=T, quote=F, row.names=F, col.names=F) Error in file(file, ifelse(append, "a", "w")) : invalid 'open&#

[R] Looping and break

2015-03-02 Thread Scott Colwell
of a which would be a = 2, and so on. for (a in 1:3){ if (a == 1) { b <- c(1:5) } if (a == 2) { b <- c(1:5) } if (a == 3) { b <- c(1:5) } for (n in 1:50){ if (n > 15) next total <- 2*b if (total > 8) next } } Any help would be grea

Re: [R] Saving Mean Relative Difference from all.equal()

2015-02-26 Thread Scott Colwell
I think I have one solution. Not very pretty though. Relies on the text not changing at all. as.numeric(gsub("Mean relative difference: ", "", all.equal(cov2cor(ITEMCOV),cor(item.data))[2])) Is there a better way? -- View this message in context: http://r.789695.n4.nabble.com/Saving-Mean-Rela

[R] Saving Mean Relative Difference from all.equal()

2015-02-26 Thread Scott Colwell
Hello, Does anyone know how to save the numeric value of the "mean relative difference" when using the all.equal() command? For example this: all.equal(cov2cor(ITEMCOV),cor(item.data)) Gives: [1] "Attributes: < Length mismatch: comparison on first 1 components >" [2] "Mean relative difference

Re: [R] Extracting Factor Pattern Matrix Similar to Proc Factor

2015-02-26 Thread Scott Colwell
Thanks everyone -- View this message in context: http://r.789695.n4.nabble.com/Extracting-Factor-Pattern-Matrix-Similar-to-Proc-Factor-tp4703704p4703904.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing lis

Re: [R] Extracting Factor Pattern Matrix Similar to Proc Factor

2015-02-23 Thread Scott Colwell
Thanks David. What do you do when the input is a covariance matrix rather than a dataset? -- View this message in context: http://r.789695.n4.nabble.com/Extracting-Factor-Pattern-Matrix-Similar-to-Proc-Factor-tp4703704p4703719.html Sent from the R help mailing list archive at Nabble.com. _

[R] Extracting Factor Pattern Matrix Similar to Proc Factor

2015-02-23 Thread Scott Colwell
ERN'; DROP _TYPE_ _NAME_; RUN; Would any SAS IML to R converts be able to help me with this? Thanks, Scott Colwell, PhD -- View this message in context: http://r.789695.n4.nabble.com/Extracting-Factor-Pattern-Matrix-Similar-to-Proc-Factor-tp4703704.html Sent from the R help mailing

Re: [R] Rewriting the Biplot Function

2015-01-20 Thread Scott Robinson
Boris, Thanks very much, this looks just like what I need! All the best, Scott From: Boris Steipe [boris.ste...@utoronto.ca] Sent: 20 January 2015 20:21 To: Scott Robinson Cc: r-help@r-project.org Subject: Re: [R] Rewriting the Biplot Function Scott

[R] Rewriting the Biplot Function

2015-01-20 Thread Scott Robinson
s outside [0, 1]") if (scale != 0) lam <- lam^scale else lam <- 1 if (pc.biplot) lam <- lam/sqrt(n) colouredBiplot.internal(t(t(scores[, choices])/lam), t(t(x$rotation[, choices]) * lam), yCol, ...) invisible() } I have looked into a few alte

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 1:15 PM, Martin Morgan wrote: > On 11/25/2014 04:11 AM, Scott Kostyshak wrote: >> >> On Mon, Nov 24, 2014 at 12:34 PM, Sarah Goslee >> wrote: >>> >>> I took a look at apparent gender among list participants a few years ago: >>

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
On Tue, Nov 25, 2014 at 8:24 AM, Maarten Blaauw wrote: > Nice graph, Scott, thanks! > > Based on your code I plotted not the absolute numbers but the ratios, which > show slowly increasing relative participation of female Rhelpers over time > (red = women, blue=men, black=unknown

Re: [R] Gender balance in R

2014-11-25 Thread Scott Kostyshak
ry low percentage of regular female participants is one of the > things that keeps me active on this list: to demonstrate that it's not > only men who use R and participate in the community. Thank you for that! Scott -- Scott Kostyshak Economics PhD Candidate Princeton University >

Re: [R] Symbolic equations to R code?

2014-11-20 Thread Scott Rifkin
bypasses the whole process. Thanks, Scott On 11/19/14 1:52 PM, Ista Zahn wrote: Hi Scott, Can you give a couple of examples of the equations you have in mind along with how those should be translated to R? Thanks, Ista __ R-help@r-project.or

Re: [R] Symbolic equations to R code?

2014-11-19 Thread Scott Rifkin
ically translate the equations of 150+ students into R code. Thanks, Scott __ 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-guide.html and provide commented, min

[R] Symbolic equations to R code?

2014-11-19 Thread Scott Rifkin
eces of it would be very much appreciated. Thanks, Scott __ 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-guide.html and provide commented, minimal, self

Re: [R] Testing general hypotheses on regression coefficients

2014-09-06 Thread Scott Kostyshak
conometrics class it should not be difficult and can provide some insights. If not, at least consider that the delta method can lead to misleading estimates (biased standard errors) in many cases for finite samples. You might want to run some simulations to get a feel for it. Best, Scott

[R] Reading chunks of data from a file more efficiently

2014-08-09 Thread Waichler, Scott R
ot;) tmp[[i]] <- scan(plt.file, skip=ind, nmax=num.nodes, quiet=T) } return(tmp) } # end READ.PLOT.OUTPUT6() Regards, Scott Waichler Pacific Northwest National Laboratory Richland, WA, USA scott.waich...@pnnl.gov __ R-help@r-project.org maili

[R] Using axis limits with plot3D

2014-08-06 Thread Waichler, Scott R
imits, scale=F, ticktype="detailed") Thanks, Scott Waichler Pacific Northwest National Laboratory Richland, WA, USA __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-

Re: [R] how to subset based on other row values and multiplicity

2014-07-16 Thread Williams Scott
Thanks guys - amazingly prompt solutions from the R community as always. Yes, the c-y value reverts to just the first date event - the spirit of this is that I am trying to identify and confirm a list of diagnoses that a patient has coded in government administrative data. Once a diagnosis is made

[R] how to subset based on other row values and multiplicity

2014-07-16 Thread Williams Scott
his via some crude loops and subsetting, but I am looking for as much efficiency as possible as the dataset has around 50 million rows to assess. Any suggestions welcomed. Thanks in advance Scott Williams MD Melbourne, Australia This email (including any attachments or links) may cont

[R] Error evaluating partitioning around medoids clustering method R clValid package

2014-07-09 Thread Scott Davis
Single No 6 13286 @gmail.com Please let me know if I can provide more information. -- Scott Davis Cell: (408)826-9561 Skype ID: Scdavis61 San Jose, CA. [[alternative HTML version deleted]] __

[R] Error creating daisy matrix in R cluster package - Cannot allocate vector of size 66.0 Gb

2014-06-21 Thread Scott Davis
type = list(ordratio = c(1:35))) #Error: cannot allocate vector of size 66.0 Gb How can I fix the error? -- Scott Davis Cell: (408)826-9561 Skype ID: Scdavis61 San Jose, CA. [[alternative HTML version deleted]] __ R-help

Re: [R] Using apply with more than one matrix

2014-05-01 Thread Waichler, Scott R
ow my computing will be more efficient. Thanks also to David and the others who responded to my question. It all helped. --Scott Waichler > -Original Message- > From: arun [mailto:smartpink...@yahoo.com] > Sent: Thursday, May 01, 2014 6:24 AM > To: R. Help > Cc: Waich

Re: [R] Using apply with more than one matrix

2014-05-01 Thread Waichler, Scott R
r that could be any length. It is a selection on the vector a1[i,j,]. I want to get the first element of that selection. The key relationship is that the dimensions of the matrices and the first two dimensions of the 3d arrays are the same. That is, i and j have the same range for all of them.

Re: [R] Using apply with more than one matrix

2014-04-30 Thread Waichler, Scott R
if(length(ind.not.na) > 0) { return(condition1*(a1[idx,ind.not.na[1]] + m2[idx])) # doesn't work; improper indexing for a1 } } vec <- sapply(seq(length(m2)), evaluate) Scott Waichler > -Original Message- > From: David Winsemius [mailto:dwinsem...@comcast.net] >

Re: [R] Using apply with more than one matrix

2014-04-30 Thread Waichler, Scott R
I tried to demonstrate in my example above. Thanks, Scott __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Using apply with more than one matrix

2014-04-30 Thread Waichler, Scott R
n1[i,j] && length(ind.not.na) > 0) { ans[i,j] <- a1[i,j,ind.not.na[1]] + m2[i,j] } } } Scott Waichler > -Original Message- > From: Bert Gunter [mailto:gunter.ber...@gene.com] > Sent: Wednesday, April 30, 2014 12:18 PM > To: Waichler, Scott R > Cc: r-he

[R] Using apply with more than one matrix

2014-04-30 Thread Waichler, Scott R
Hi, I want to apply a custom function to all the elements of one matrix. The function uses the value in the same i,j in a second matrix. How can I use apply() or similar function to avoid nested loops in i and j? Thanks, Scott Waichler Pacific Northwest National Laboratory Richland, WA, USA

Re: [R] problem with pip2d() from ptinpoly

2014-04-19 Thread Waichler, Scott R
Thank you for catching that, Boris. I'm surprised, given that there is no mention of sense of direction in the package documentation. Scott Waichler > -Original Message- > From: Boris Steipe [mailto:boris.ste...@utoronto.ca] > Sent: Friday, April 18, 2014 3:49 PM > To

[R] problem with pip2d() from ptinpoly

2014-04-18 Thread Waichler, Scott R
;- matrix(c(594885.0,115435.0), ncol=2, byrow=T) pip2d(Vertices = verts, Queries = query) # result = -1 # contrary to -1 output of pip2d, plot shows point lies within triangle plot(c(594400, 595000), c(115000, 117500), type="n") polygon(verts, border="red") points(x=query[,1], y=query[,

[R] Unexpected behaviour with seq() and %in%

2013-08-23 Thread Finlay Scott
Hi, During the course of putting together a function I came across some unexpected behaviour when using seq() and %in%. I am creating two numeric vectors using seq(), and then using %in% to find the values in one vector that are in the other. Sometimes all the values are found, but sometimes a val

Re: [R] BH correction with p.adjust

2013-07-21 Thread Scott Robinson
)/1 [1] 0.991 Any help with this would be very much appreciated. It seems like it ought to be such a simple and commonly used method and yet I am struggling and not sure what to do about it. Thanks, Scott From: David Winsemius [dwinsem...@comcast.net]

[R] BH correction with p.adjust

2013-07-20 Thread Scott Robinson
ted there. I was wondering if there is some subtle difference about how p.adjust operates? Is there something more complicated about how to calculate 'n' or 'i' - perhaps due to identical p values being assigned the same rank or something? Does anyone have an idea what might be

Re: [R] Survey imputation

2013-06-13 Thread Scott Raynaud
I paln on using R to do the imputation once I figure out how to do it.   - Original Message - From: Bert Gunter To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Thursday, June 13, 2013 12:45 PM Subject: Re: [R] Survey imputation Is this an R question? Seems like it be

[R] Survey imputation

2013-06-13 Thread Scott Raynaud
I'm working with NHIS survye data.  I'd like the to use muliple imputation to cover the missing data for the variables in which I'm interested.  My question concerns the use of certain variables in the imputation model.  For example, race would be an important predictor in the imputation model,

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
for helping to diagnose and fix the problem. - Original Message - From: Ista Zahn To: Scott Raynaud Cc: William Dunlap ; "r-help@r-project.org" Sent: Thursday, June 6, 2013 9:15 AM Subject: Re: [R] SPlus script Presumably something like r <- sshc(50) print(r) But if you wer

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
how to call the function appreciated. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2:17 PM Subject: RE: [R] SPlus script Both the R and S+ versions (which seem to differ only in the use of _ for assignm

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
Ok. I tried copying the original program, changing all _ to <- and running in 3.0.1.  Still no error message or output.  It's a mystery to me. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2

Re: [R] SPlus script

2013-06-05 Thread Scott Raynaud
the historical group.  Maybe I'm forgetting something here...  It's been a while since I ran this.  Thanks for your help. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2:17 PM Subject: RE: [R] SPlus

Re: [R] SPlus script

2013-06-05 Thread Scott Raynaud
See my comments below. From: Pascal Oettli To: Scott Raynaud Sent: Wednesday, June 5, 2013 10:02 AM Subject: Re: [R] SPlus script Hello, 1) It is always nice to say something as "Hello", Tried, but could make it come out in anything other than fo

[R] SPlus script

2013-06-05 Thread Scott Raynaud
This originally was an SPlus script that I modifeid about a year-and-a-half ago.  It worked perfectly then.  Now I can't get any output despite not receiving an error message.  I'm providing the SPLUS script as a reference.  I'm running R15.2.2.  Any help appreciated.    ***

[R] vector field from a 3D scalar field

2013-03-21 Thread Waichler, Scott R
I have a 3D field of a scalar variable (x, y, z, value). Is there a way to generate a vector field from this data--gradient at defined points? I found the rasterVis package for 2D data, but as yet nothing for 3D data. Thanks, Scott Waichler Pacific Northwest National Laboratory Richland, WA

[R] how to get xmlToList() to retry if http fails

2013-03-06 Thread Waichler, Scott R
he webservice failed to respond in this instance. How can I trap the error and have it retry after waiting a second or two, instead of exiting? Thanks. --Scott Waichler Pacific Northwest National Laboratory Richland, WA, USA scott.waich...@pnnl.gov __

[R] Iterating through slots of an S4 object

2013-02-15 Thread Scott Robinson
successful. I presume it is possible though? Any help would be much appreciated, Scott [[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 ht

[R] WriteXLS: 'object not found' error within function

2013-02-13 Thread Scott Robinson
alling the WriteXLS function: print(class(tables)) flush.console() [1] "list" At least I would have thought this would have ruled out any scope-related issues. Has anyone else had this problem or have any ideas why it might be happening? Thanks, Scott PS here is the function in full:

Re: [R] sweave and lessR

2012-12-16 Thread Scott Cushman
Thank you for the reply. With many calls to Density() I miss Sweave taking care of the naming of the files and the generating of the \includegraphics statements, but I guess this will work. On Dec 16, 2012, at 2:54 PM, Duncan Murdoch wrote: > On 12-12-16 12:42 PM, Scott Cushman wr

[R] sweave and lessR

2012-12-16 Thread Scott Cushman
Hello I am using sweave with the Density() function from lessR package. However, unlike the older color.density() function, Density() does not work with the standard graphic output functions in R, such as pdf. Is there away to include figures fgenerated rom lessR in sweave documents? thank yo

[R] analysis of bitmaps

2012-09-05 Thread Scott Swigart
I'm interested in using R to perform some statistical analysis of bitmaps. When I search, I see a lot of information about outputting bitmaps, but I'm not finding much about loading a bitmap into a data frame so that it can be analyzed. Would someone have a quick pointer to help me out?

Re: [R] generating random samples of IG distribution

2012-08-29 Thread David Scott
I don't read R-help these days so have just seen this. Both generalized inverse Gaussian and normal inverse Gaussian are in GeneralizedHyperbolic. HyperbolicDist is no longer being maintained. David Scott On 12/06/2012 5:41 a.m., David L Carlson wrote: Should have been For the n

Re: [R] Error message

2012-08-29 Thread Scott Raynaud
rson so someone will have to tell me the command to obtain the Linux information requested.  Then maybe I can figure whether to go to the Debian board or not. - Original Message - From: Duncan Murdoch To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Tuesday, August 28, 20

Re: [R] Error message

2012-08-28 Thread Scott Raynaud
  - Original Message - From: Duncan Murdoch To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Tuesday, August 28, 2012 10:40 AM Subject: Re: [R] Error message On 28/08/2012 10:31 AM, Scott Raynaud wrote: > I suddenly started getting the error message below. > Not sure why.

[R] Error message

2012-08-28 Thread Scott Raynaud
I suddenly started getting the error message below. Not sure why.  If I type intalled.packages() it shows Matrix and lme4 installed.  Can someone tell what's going on and what I need to do to remedy the problem?  I'm running on a Linux box.             Loading required package: Matrix Loading requ

Re: [R] RGDAL OGRwrite question

2012-08-17 Thread Scott Duke-Sylvester
ermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > > > > On 8/17/12 9:28 AM, "Scott Duke-Sylvester" > wrote: > >>I have a quick question: It appears that in rgdal v0.7-12 (R version >>2.15.1, OSX

[R] RGDAL OGRwrite question

2012-08-17 Thread Scott Duke-Sylvester
1] "ID" "ANID" "F_AREA" "Avg_z" "Manning" [6] "IniWL" "IniSal" "DispCoeff" "HydRad" "veg1_DW" [11] "veg2_IW""veg3_SW""veg4_DC&q

Re: [R] Accents and special character using hwriter (on Windows)

2012-08-07 Thread David Scott
character codes, see for example http://www.ascii.cl/htmlcodes.htm David Scott On 31/07/2012 9:21 p.m., ramonovelar wrote: > Thanks Arun, > Yes, I have Windows 7. I have tried 2 versions of R, 2.14.1 and 2.15.x, but > it did not change anything. > Right now I can't try a di

[R] Can't read a binary file

2012-04-13 Thread Waichler, Scott R
using double): > v<-readBin("plotb.251", numeric(), size=4, n=1) > v [1] 1.614296e-39 > v<-readBin("plotb.251", numeric(), size=4, n=1, endian="swap") > v [1] 1.359775e-38 Platform is Intel Linux. How can I read the file described above? Thanks, S

Re: [R] Can't get R to recognize Java for rJava installation

2012-04-13 Thread Waichler, Scott R
w are the cpp flags: cpp flags : '' Could this be the problem now? How can I set those, and what value should I give? Scott Waichler > So I guess you need to find out what package provides this file on your > distribution (which you did not mention). First check the file is &g

[R] Can't get R to recognize Java for rJava installation

2012-04-12 Thread Waichler, Scott R
variables are not set. Make sure R is configured with full Java support (including JDK). Run R CMD javareconf as root to add Java support to R. If you don't have root privileges, run R CMD javareconf -e to set all Java-related variables and then install rJava. ERROR: configuration faile

Re: [R] Apt-get

2012-03-16 Thread Scott Raynaud
Ok.  When I type in sudo add-apt-repository ppa:marutter/rrutter I get a prompt for a password.  I enter my domain password and get some print about adding the ppa and requeting that I press enter.  I press enter and this is followed by print stating that the key B04C661B is being requested from h

Re: [R] Apt-get

2012-03-15 Thread Scott Raynaud
iller To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, March 7, 2012 12:34 PM Subject: Re: [R] Apt-get Google is really useful for questions like this. http://cran.r-project.org/bin/linux/ubuntu/ -

[R] Apt-get

2012-03-07 Thread Scott Raynaud
I have a box set up with Kubuntu as the OS.  I didn't perform the R install but was told the version of R available via the apt-get command was 2.13.1.  Is there any way to get 2.14.0 in that same manner? __ R-help@r-project.org mailing list https://sta

[R] Access to OpenBLAS

2012-02-14 Thread Scott Raynaud
My IT people have set up a Kubuntu box with an RWkard front end.  I have OpenBLAS set up as a shared BLAS but I'm not sure how to get R to see it.  A.3.1 of the installation docs talks about it but I'm not clear if I need a option on my startup line or if I need to find a config file.  The BLAS is

[R] Simulation confidence interval

2012-02-01 Thread Scott Raynaud
The follwing is a code snippet from a power simulation program that I'm using:   estbeta<-fixef(fitmodel)  sdebeta<-sqrt(diag(vcov(fitmodel)))   for(l in 1:betasize)   {     cibeta<-estbeta[l]-sgnbeta[l]*z1score*sdebeta[l]     if(beta[l]*cibeta>0)  powaprox[[l]]<-powaprox[[l]]+1  

[R] Recompile

2012-01-23 Thread Scott Raynaud
Looked at several posts and the installtions docs and still not clear.  If I compile source codes and then somewhere down the line add a new package, then I have to recompile my entire installation, correct?  Seems like this is the sentiment of the emails I read. ___

[R] GUI

2012-01-18 Thread Scott Raynaud
I'm setting up a Linux box to run R.  I ususally run in a Windows envrionment but after reading the docs I'm not sure what to expect in terms of the front end appearance in Linux.  Does it resemble Windows or will I need Rkward or R Commander? __ R-h

[R] BLAS

2012-01-17 Thread Scott Raynaud
I'm setting up an Ubuntu virtual machine that will use 4-Intel Xeon CPU x5650.  I'd like to compile R with a BLAS but the question is whcih one.  Seems like the only free ones are GotoBLAS which I'm not sure is being maintained for newer CPUs and OpenBLAS for Loongson CPUs.  I saw a favorable re

Re: [R] grplasso

2012-01-13 Thread Scott Raynaud
So does anyone use this package?   - Original Message - From: Scott Raynaud To: "r-help@r-project.org" Cc: Sent: Tuesday, January 10, 2012 1:40 PM Subject: grplasso I want to use the grplasso package on a data set where I want to fit a linear model.  My interest is in i

[R] grplasso

2012-01-10 Thread Scott Raynaud
I want to use the grplasso package on a data set where I want to fit a linear model.  My interest is in identifying significant beta coefficients.  The documentation is a bit cryptic so I'd appreciate some help.   I know this is a strategy for large numbers of variables but consider a simple cas

Re: [R] R report generator (for Word)?

2012-01-02 Thread David Scott
the package. David Scott From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of Joshua Wiley [jwiley.ps...@gmail.com] Sent: Monday, January 02, 2012 9:31 AM To: Michael Cc: r-help Subject: Re: [R] R report generator (for Word

Re: [R] rbinom

2011-12-30 Thread Scott Raynaud
but I'd need to square my x[,3] values before multiplying them by beta.  Can I say: x[,3]<-(rnorm(length,meanpred[3],sqrt(varpred[3])))^2 in lieu of x[,3]<-rnorm(length,meanpred[3],sqrt(varpred[3]))? - Original Message - From: peter dalgaard To: Scott Raynaud Cc: "r-hel

[R] rbinom

2011-12-27 Thread Scott Raynaud
I have the following code (which I did not write) that generates data based on a logistic model.  I'm only getting a single record with y=1.  It seems implausible that in 50k cases that have a single y=1.  Does that ring alarm bells for anyone else?   beta<-c(-1.585600,-0.246900) betasize<-lengt

[R] simulation

2011-12-16 Thread Scott Raynaud
I'm using an R program (which I did not write) to simulate multilevel data (subjects in locations) used in power calculations. It uses lmer to fit a mixed logistic model to the simulated data based on inputs of means, variances, slopes and proportions:   (fitmodel <- lmer(modelformula,data,famil

[R] minor allele frequency comparison

2011-12-08 Thread Scott Harding
C/T 75/13 1460/280 chromosome 1 14282752 A/G 69/19 1380/360 chromosome 1 595434T/G 83/5 1680/60 I really appreciate any guidance. Thanks Scott

Re: [R] R/parallel

2011-12-08 Thread Scott Raynaud
recompiling with BLAS and pnmath would be a better option.   My main question is how to handle the random number generation when the child processes are spawned.  That's a problem no matter what method I choose to create the threads.   From: Tal Galili To: Scott Raynaud Cc: &quo

  1   2   3   4   5   >