example of using
lattice to plot Simple Feature polygons? I'm still learning sf and can't find
any advice for using it with lattice.
Thank you,
Scott Waichler, PhD
Pacific Northwest National Lab
Washington, USA
[[alternative HTML versi
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
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
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
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
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)
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
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
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
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
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
;)
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
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
_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
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
,
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
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,
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
__
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
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
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
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
_
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
?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
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
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
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
'
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
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
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
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
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
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.
_
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
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
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
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:
>>
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
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
>
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
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
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
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
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
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-
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
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
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]]
__
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
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
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.
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]
>
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.
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
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
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
;- 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[,
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
)/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]
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
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
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,
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
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
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
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
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
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.
***
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
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
__
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
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:
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
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
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?
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
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
- 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.
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
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
1] "ID" "ANID" "F_AREA" "Avg_z" "Manning"
[6] "IniWL" "IniSal" "DispCoeff" "HydRad" "veg1_DW"
[11] "veg2_IW""veg3_SW""veg4_DC&q
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
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
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
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
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
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/
-
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
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
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
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.
___
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
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
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
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
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
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
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
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
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
1 - 100 of 469 matches
Mail list logo