In this particular case it is part of the C99 stanadrd (7.12.7.4) for
the 'pow' function R uses.
On Wed, 27 Oct 2010, Berwin A Turlach wrote:
G'day Gregory,
On Tue, 26 Oct 2010 19:05:03 -0400
Gregory Ryslik wrote:
Hi,
This might be me missing something painfully obvious but why does the
cu
Hello useRs,
I'm pleased to announce the general availability of the R package
futile.paradigm, which is a language extension that implements
functional dispatching in R. This is an alternative to the current
object-oriented styles, replacing them with a functional programming
style that provides
On Tue, 26 Oct 2010, Kjetil Halvorsen wrote:
I have been searching the documentation for some time nmow, but
cannot find it. It must be possible to download packages (many), but
only install the help system?
I am not sure what you actally want: R CMD INSTALL does not install
the Rd files --
Actually I want to see how close are some point to a line segment so I want
to
use some zoom lenses and zoom in and out into different parts of the plot and
see how some places look like.
Best Regards
From: Greg Snow
Sent: Tue, October 26, 2010 6:47:32 P
?prop.test
--
View this message in context:
http://r.789695.n4.nabble.com/clinical-data-analysis-tp3014811p3014856.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
This probably fixes it, although an example would be nice :
write.csv(dfr, row.names=FALSE)
you change the names with the ?names function before writing it to a csv
file.
remko
--
View this message in context:
http://r.789695.n4.nabble.com/writing-a-table-tp3014821p3014839.html
Sent from the
Hi,
I have tab-delimited data with an unequal number of entries per column, of
the sort:
A B C
1 2 2
3 4 1
5 2 2
6 2
5 2
3
6
2
I would like to make a histogram of the frequencies of each represented
number in a "stacked" histogram, where you can see the contribution of eac
When ever I write a table of columns to a csv file the columns are offset by
one column and do not match the column headers. How do I align my columns
with the column headers? Also, how do I change the name of my column
headers.
--
View this message in context:
http://r.789695.n4.nabble.com/wri
Hi R-users
I am trying to use the GAM function of the mgcv package. But I am having
problem trying to specify the k parameter.
Although I managed to run some models by giving to the parameter some
(random) value, and it is explained by Wood (2006) that it does not seem
to "really" affect th
I have two sets of data
national wide: total patient number=500,000, 8400 died, mortality rate=1.7%
Hospital1: total patient number=230, 6 died, mortality rate=2.6%
Is the mortality rate of hospital1 different from mortality rate of nationl
wide?
Which function or package should I use to solve
On Tue, Oct 26, 2010 at 8:37 PM, Matt Curcio wrote:
> Hi All,
> I am learning R and having a little trouble with the usage and proper
> definitions of data.frames vs. matrix vs vectors. I have read many R
> tutorials, and looked over ump-teen 'cheat' sheets and have found that
> no one has articul
G'day Bill,
On Wed, 27 Oct 2010 10:34:27 +1100
wrote:
[...]
> It is no surprise that this does not work when working in the real
> domain, except "by fluke" with something like
>
> > -4^(1/3)
> [1] -1.587401
> >
>
> where the precedence of the operators is not what you might expect.
> Now th
Hi All,
I am learning R and having a little trouble with the usage and proper
definitions of data.frames vs. matrix vs vectors. I have read many R
tutorials, and looked over ump-teen 'cheat' sheets and have found that
no one has articulated a really good definition of the differences
between 'data.
Because it is implemented as
antilog((1/3)*log(-4))
most likely using base 2 for the log/antilog functions.
"Gregory Ryslik" wrote:
>Hi,
>
>This might be me missing something painfully obvious but why does the
>cube root of the following produce an NaN?
>
>> (-4)^(1/3)
>[1] NaN
>>
>
>As we ca
install.packages('sos')# if you don't have it already
library(sos)
rs <- ???roots
# 216 matches
summary(rs)
# in 106 packages
rs # opens a web browser with all 216 matched in a table
# listing the package with the most matches first.
# This included roots{signal}, which referenced polyroot{base},
G'day Gregory,
On Tue, 26 Oct 2010 19:05:03 -0400
Gregory Ryslik wrote:
> Hi,
>
> This might be me missing something painfully obvious but why does the
> cube root of the following produce an NaN?
>
> > (-4)^(1/3)
> [1] NaN
1/3 is not exactly representable as a binary number. My guess is tha
I have been searching the documentation for some time nmow, but cannot find it.
It must be possible to download packages (many), but only install the
help system?
How?
Kjetil
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-
Try this:
cut(DF$A, c(-Inf, 0, 4, 8, 12, 14, 16, 20))
On Tue, Oct 26, 2010 at 7:01 PM, LCOG1 wrote:
>
> Hi everyone,
>
> I have been using R too long to have to ask this but here i am. How do i
> create a separate bin for the 0 value? So for:
>
> #Create data frame
> DF<-data.frame(A=0:20)
On Tue, Oct 26, 2010 at 7:17 PM, Gabor Grothendieck
wrote:
> On Tue, Oct 26, 2010 at 3:28 PM, Susanta Mohapatra
> wrote:
>> Hi,
>>
>> I am working with a dataset for sometime and I need some help in parsing
>> some data.
>>
>> There is a column called "Duration" which has data like following:
>>
To take it one step further:
> x <- as.complex(-4)
> cx <- x^(1/3)
>
> r <- complex(modulus = Mod(cx), argument = Arg(cx)*c(1,3,5))
> r
[1] 0.793701+1.37473i -1.587401+0.0i 0.793701-1.37473i
> r^3
[1] -4+0i -4+0i -4+0i
>
So when you ask for "the" cube root of -4, R has a choice of three p
hmm interesting. When I did -4^(1/3) got the correct answer, but then again
that's because it processes the negative later. i.e. -4^(1/2) gave me -2
instead of the 2i I expected. Also when I did (-4+0i)^(1/3) it gave me
0.793701+1.37473i. Possible bug?
Sachin
--- Please consider the environment b
On Tue, Oct 26, 2010 at 3:28 PM, Susanta Mohapatra
wrote:
> Hi,
>
> I am working with a dataset for sometime and I need some help in parsing
> some data.
>
> There is a column called "Duration" which has data like following:
>
> 2 minutes => 120
> 2 min => 120
> 10 seconds =>10
> 2 hrs =>7200
> 2
Look at this:
> x <- as.complex(-4)
> x
[1] -4+0i
> x^(1/3)
[1] 0.793701+1.37473i
> (-4)^(1/3)
[1] NaN
It seems that R gives you the principal root, which is complex, and
not the real root.
Kjetil
On Tue, Oct 26, 2010 at 8:05 PM, Gregory Ryslik wrote:
> Hi,
>
> This might be me missing somethi
On Tue, 26 Oct 2010, Marc Schwartz wrote:
On Oct 26, 2010, at 2:59 PM, Ben Hunter wrote:
Hi everyone,
Why am I having such a tough time finding a way to put an mlogit summary
table into latex? Everywhere I read says that using Sweave and latex is the
most sophisticated, dynamic way to get out
Hi,
This might be me missing something painfully obvious but why does the cube root
of the following produce an NaN?
> (-4)^(1/3)
[1] NaN
>
As we can see:
> (-1.587401)^3
[1] -4
Thanks!
Greg
__
R-help@r-project.org mailing list
https://stat.ethz.c
> 1. What is everyone else using? The network effect is important since
> you want people to be able to access your repository and you want to
> leverage your knowledge of the version control system for other
> projects' repositories. To that extent Subversion is the clear choice
> since its used
> The caret package has answers to all your questions.
>> 1) How to obtain a variable (attribute) importance using
>> e1071:SVM (or other
>> svm methods)?
I haven't implemented a model-specific method for variables importance
for SVM models. I know of one package (svmpath) that will return the
re
On Tue, Oct 26, 2010 at 12:16 PM, Tal Galili wrote:
> Hello all,
>
> I wish to learn a version control system for managing my R (data analysis)
> projects.
>
> I know of SVN and github, and wonder if there is any reason for which I
> should prefer the one over the other (or any other platform). An
On Tue, Oct 26, 2010 at 3:13 PM, Daisy Englert Duursma
wrote:
> There are several ways to do this but the package ggplot2
>
> library(ggplot2)
> qplot(displ,hwy,data=mpg,colour=factor(cyl))
>
>
That can of course be done also using the standard plot command
(substitute variable names as necessar
There are several ways to do this but the package ggplot2
library(ggplot2)
qplot(displ,hwy,data=mpg,colour=factor(cyl))
On Wed, Oct 27, 2010 at 9:06 AM, elaine kuo wrote:
> Dear List,
>
>
>
> I am using the command plot to present the relationship
>
> between bird richness (Y axis) and elevatio
On Oct 26, 2010, at 2:59 PM, Ben Hunter wrote:
> Hi everyone,
>
> Why am I having such a tough time finding a way to put an mlogit summary
> table into latex? Everywhere I read says that using Sweave and latex is the
> most sophisticated, dynamic way to get output, but it appears very limited
> i
Dear List,
I am using the command plot to present the relationship
between bird richness (Y axis) and elevation (X axis).
However, I would like to observe
the distributions of bird richness in different administrative areas (A, B,
C,
., G) in this plot.
For example, the dots in area A mig
I don't know why I forgot that you can do this as well :
area.poly(intersect(p1,p2))
... a bit more straightforward.
greetings,
Remko
--
View this message in context:
http://r.789695.n4.nabble.com/calculate-area-between-intersecting-polygons-tp3012980p3014581.html
Sent from the R help mailin
On 10/26/2010 2:33 PM, Bond, Stephen wrote:
How can I expose the code behind plot.survfit??
Uwe Ligges. R Help Desk: Accessing the sources. R News, 6(4):43-45,
October 2006.
http://cran.r-project.org/doc/Rnews/
Thanks a lot.
Stephen B
[[alternative HTML version deleted]]
--
Bond, Stephen cibc.com> writes:
>
> How can I expose the code behind plot.survfit??
library(survival)
survival:::plot.survfit
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://w
Peter Francis me.com> writes:
>
> Dear List,
>
> I am looking to plot error bars on a line using dispersion.
>
> I have values for the upper value and for the lower values, however
> i am unsure how to plot different values
> for the upper CI and the lower CI?
>
> I have been using
>
> disp
Thanks a lot, Douglas. It's very heplful.
A clarification question about specifying the model in lmer. You said
it should be:
mix.lmer <- lmer(DV ~a+b+c+d+(e+f+g+h+ii|group), mydata)
I assume it was a typo and you meant that the last predictor in
brackets should be i (rather than ii), right?
Als
On 26/10/2010 5:17 PM, Xiuquan Wang wrote:
Hi,
I am now using R to implement a stepwise algrithom which includes a
recursion function.
e.g:
a = 1
*f_recursion* = function(id)
{
b = a + id; #: row A
if (...){ a = a +1;* f_recursion*(b) } #: row B
How can I expose the code behind plot.survfit??
Thanks a lot.
Stephen B
[[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 http://www.R-proje
Hi,
I am now using R to implement a stepwise algrithom which includes a
recursion function.
e.g:
a = 1
*f_recursion* = function(id)
{
b = a + id; #: row A
if (...){ a = a +1;* f_recursion*(b) } #: row B
else{ }
}
--
Hi,
I am working with a dataset for sometime and I need some help in parsing
some data.
There is a column called "Duration" which has data like following:
2 minutes => 120
2 min => 120
10 seconds =>10
2 hrs =>7200
2-3 minutes => 150 or 120
5 minutes (when i arrived => 300
Flyby approx 20 sec. =
> Date: Tue, 26 Oct 2010 12:53:14 -0400
> From: mike...@gmail.com
> To: j...@bitwrit.com.au
> CC: r-help@r-project.org
> Subject: Re: [R] runtime on ising model
>
> I have an update on where the issue is coming from.
>
> I commented out the code for "
Hi everyone,
I have been using R too long to have to ask this but here i am. How do i
create a separate bin for the 0 value? So for:
#Create data frame
DF<-data.frame(A=0:20)
#Create label vector
labs<-1:6
#Create buckets and label
DF$Cut<-cut(DF$A,c(0,4,8,12,14,16,20),labels=labs,include.l
On 10/25/2010 8:56 PM, Daisy Englert Duursma wrote:
Hello,
If I have a dataframe:
example(data.frame)
zz<-c("aa_bb","bb_cc","cc_dd","dd_ee","ee_ff","ff_gg","gg_hh","ii_jj","jj_kk","kk_ll")
ddd<- cbind(dd, group = zz)
and I want to divide the column named group by the "_", how would I do this?
On 2010-10-26 11:48, Jonathan P Daily wrote:
?loess
use this instead:
fit<- loess(b~a)
lines(a, predict(fit))
I don't think that will work when there are incomplete cases,
in which case 'a' and predict(fit) may not correspond.
I think that it's always best to define a set of predictor
values
Thanks David!
After setting the LD_LIBRARY_PATH variable(to /usr/local/lib) R was able to
successfully install the ncdf4 package.
On Tue, Oct 26, 2010 at 2:24 PM, David Pierce [via R] <
ml-node+3014258-455613790-200...@n4.nabble.com
> wrote:
> shaticus wrote:
>
> >
> > Hello all,
> >
> > I coul
Hi everyone,
Why am I having such a tough time finding a way to put an mlogit summary
table into latex? Everywhere I read says that using Sweave and latex is the
most sophisticated, dynamic way to get output, but it appears very limited
in this respect. I'm just starting out with Sweave and LaTeX
On Tue, Oct 26, 2010 at 12:27 PM, Dimitri Liakhovitski
wrote:
> Hello,
> and sorry for asking a question without the data - hope it can still
> be answered:
> I've run two things on the same data:
> # Using lme:
> mix.lme <- lme(DV ~a+b+c+d+e+f+h+i, random = random = ~ e+f+h+i|
> group, data = m
Hello
I am trying to use sbrier in ipred but got an error message (see below).
Can someone help?
---
I. function()
{
library(ipred)
library(survival)
set.seed(12345)
age <- rnorm(30, 50, 10)
stime <- rexp(30)
cens <- runif(30,.5,2)
sevent <- as.numeric(stime <= cens)
stime <- pmin(stime,
Is it possible to get survfit to produce the survival line for a single strata
like
preddow <- survfit(modall,newdata=newdat,se.fit=F,strata=2)
# the strata argument is being ignored in the call above
Or even get a more economical/faster calculation of the hazard directly from
the coxph object
On Tue, Oct 26, 2010 at 11:55 AM, Dennis Murphy wrote:
> Hi:
>
> When it comes to split, apply, combine, think plyr.
>
> library(plyr)
> ldply(split(afvtprelvefs, afvtprelvefs$basestudy),
> function(x) coef(lm (ef ~ quartile, data=x, weights=1/ef_std)))
Or do it in two steps:
models <- d
On Oct 26, 2010, at 8:57 AM, sr500 wrote:
Hello,
I was wondering if anyone knew of a function that fits haplotype
data into a
cox proportional hazard model. I have computed my Haplotype
frequencies
using the haplo.stats package. I have also been using the haplo.glm
function
but this is
?loess
use this instead:
fit <- loess(b~a)
lines(a, predict(fit))
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a room when its empty? Does the room,
the thing itself
On Oct 26, 2010, at 11:15 AM, Gavin Simpson wrote:
On Fri, 2010-10-22 at 15:39 +0200, Claudia Beleites wrote:
On 10/22/2010 03:15 PM, DrCJones wrote:
Being a chemist, it seemed natural to me to put the i after the
concentration
brackets into a subscript - though you didn't say you want
> git is where the world is headed. This video is a little old:
> http://www.youtube.com/watch?v=4XpnKHJAok8, but does a good job
> getting the point across.
And lots of R users are using github already:
http://github.com/languages/R/created
Hadley
--
Assistant Professor / Dobelman Family Juni
Hi:
When it comes to split, apply, combine, think plyr.
library(plyr)
ldply(split(afvtprelvefs, afvtprelvefs$basestudy),
function(x) coef(lm (ef ~ quartile, data=x, weights=1/ef_std)))
.id (Intercept) quartile
1 CBP090802020.92140 3.38546887
2 CBP090802129.31632 0.013
I have an update on where the issue is coming from.
I commented out the code for "pos[k+1] <- M[i,j]" and the if statement for
time = 10^4, 10^5, 10^6, 10^7 and the storage and everything ran fast(er).
Next I added back in the "pos" statements and still runtimes were good
(around 20 minutes).
So
Hello,
I was wondering if anyone knew of a function that fits haplotype data into a
cox proportional hazard model. I have computed my Haplotype frequencies
using the haplo.stats package. I have also been using the haplo.glm function
but this is a linear regression and is not quite what I am look
Hello Masters,
I run the loess() function to obtain local weighted regressions, given
lowess() can't handle NAs, but I don't
improve significantly my situation.., actually loess() performance leave
me much puzzled
I attach my easy experiment below
#--SCRIPT---
Hi all,
I have a problem with this code... as it generates an error in R.
z<-predict(dat.fit, newdata=grd)
I saw a post on R forum about this ["chfactor.c", line 130: singular matrix in
function LDLfactor() ] error and tried pretty much anything I could read about
it and still hav
On Fri, 2010-10-22 at 05:54 -0700, Penny Adversario wrote:
> I am doing cluster analysis on 8768 respondents on 5 lifestyle
> variables and am having difficulty constructing a dissimilarity matrix
> which I will use for PAM. I always get an error: “cannot allocate
> vector of size 293.3 Mb” even
shaticus wrote:
>
> Hello all,
>
> I could use some help installing the ncdf4 package in R (under CentOS
> 5.4).
>...
> When I issue an 'install.packages(c('ncdf4')) in R, compilation succeeds
> but
> I receive the following error during the loading phase of installation:
>
> ** building package in
On Fri, 2010-10-22 at 15:39 +0200, Claudia Beleites wrote:
> On 10/22/2010 03:15 PM, DrCJones wrote:
>
> Being a chemist, it seemed natural to me to put the i after the concentration
> brackets into a subscript - though you didn't say you want that.
>
> A more "correct" expression would be:
>
Hello all,
I could use some help installing the ncdf4 package in R (under CentOS 5.4).
I've installed R-2.12.0, zlib-1.2.5, hdf5-1.8.4-patch1 and NetCDF4.1.1 from
source. Make check reports all tests have passed for all of these programs.
When I issue an 'install.packages(c('ncdf4')) in R, comp
Try this:
sapply(by(x, x$basestudy, lm, formula = ef ~ quartile), coef)
On Tue, Oct 26, 2010 at 1:08 PM, Small Sandy (NHS Greater Glasgow & Clyde) <
sandy.sm...@nhs.net> wrote:
> Hi
>
> I need some help getting results from multiple linear models into a
> dataframe.
> Let me explain the problem.
Marc is exactly right about people having strong opinions.
R-forge is really the _only_ reason to consider using svn.
git is where the world is headed. This video is a little old:
http://www.youtube.com/watch?v=4XpnKHJAok8, but does a good job
getting the point across.
Hg is a good alternative,
On Oct 26, 2010, at 10:22 AM, Small Sandy (NHS Greater Glasgow &
Clyde) wrote:
Thanks David
That's great
As a matter of interest, to get a data frame by studies why do you
have to do
fitsdf <- as.data.frame(t(as.data.frame(fits)))
The apply family of functions often return results rot
On Oct 26, 2010, at 9:27 AM, David Smith wrote:
Many thanks for the help.
You could express your thanks by including context the next time you
present a follow-up (as requested in the Posting Guide). Only a
minority of readers view this list on Nabble, so we don't see the web
delivered
Hello,
and sorry for asking a question without the data - hope it can still
be answered:
I've run two things on the same data:
# Using lme:
mix.lme <- lme(DV ~a+b+c+d+e+f+h+i, random = random = ~ e+f+h+i|
group, data = mydata)
# Using lmer
mix.lmer <- lmer(DV
~a+b+c+d+(1|group)+(e|group)+(f|grou
Thanks David
That's great
As a matter of interest, to get a data frame by studies why do you have to do
fitsdf <- as.data.frame(t(as.data.frame(fits)))
Why doesn't
fitsdf <- as.data.frame(t(fits))
work?
Sandy Small
From: David Winsemius [dwinsem...@comc
On Oct 26, 2010, at 11:22 AM, Duncan Murdoch wrote:
> On 26/10/2010 12:16 PM, Tal Galili wrote:
>> Hello all,
>>
>> I wish to learn a version control system for managing my R (data analysis)
>> projects.
>>
>> I know of SVN and github, and wonder if there is any reason for which I
>> should pref
On Tue, Oct 26, 2010 at 12:56 PM, Manta wrote:
>
> Thanks for your help Gabor. That would be exactly what I am looking for. If I
> use your code I get the nice representation I am looking for. However, when
> I try to apply the code in the same fashion to my case, it does not produce
> the x-axis.
Thanks for your help Gabor. That would be exactly what I am looking for. If I
use your code I get the nice representation I am looking for. However, when
I try to apply the code in the same fashion to my case, it does not produce
the x-axis. I believe the problem hinges on the following warning me
The caret package has answers to all your questions.
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Neeti
> Sent: Tuesday, October 26, 2010 10:42 AM
> To: r-help@r-project.org
> Subject: [R] to determine the variable importa
For a quick exploration of the plot you can use the zoomplot function in the
TeachingDemos package. But for production graphs it is better to explicitly
set the xlim and ylim parameters in creating the plot up front.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthc
We can further generalize this:
Suppose we want to constrain parameters such that :
b2 = a * b1
b3 = a * b1
We can do the following:
fit.a <- glm( y ~ I(x1 + a* x2 + a * x3), data= , ... )
For a fixed value of `a', we compute the log-likelihood of `fit.a'. This is
the profile likeli
Thanks a lot - it's very helpful.
On Tue, Oct 26, 2010 at 11:37 AM, Ista Zahn wrote:
> Hi Dimitri,
> The lme function is not in the lme4 package, so there is some
> confusion there. But you can use weights with the lmer function in
> lme4. ?lmer tells you that weights are specified the same way
Many thanks for the help.
I assumed that I would need to account for the variables in the model, even
though I wish to assign a zero coefficient to them. I've looked at the
offset function, but does this not just assign the value 1 to the variables?
How would I specify a zero coefficient to more
I would still recommend
vector_of_column_number <- apply(yourdata, 1, match, x=1)
as the simplest way if you only want the number of the
column that has the first 1 or "1" (the call works as is
for both numeric and character data). Rows which have no
1s will return a value of NA.
Anything wron
On 26/10/2010 12:16 PM, Tal Galili wrote:
Hello all,
I wish to learn a version control system for managing my R (data analysis)
projects.
I know of SVN and github, and wonder if there is any reason for which I
should prefer the one over the other (or any other platform). An example for
a reason
On 26/10/2010 11:30 AM, Alaios wrote:
in a simple plot. When i do plot is it possible to zoom in or out or this is not
possible at all?
There's no general support for that, but you could conceivably write it
yourself using getGraphicsEvent. The example code there adjusts xlim
and ylim accord
Hello all,
I wish to learn a version control system for managing my R (data analysis)
projects.
I know of SVN and github, and wonder if there is any reason for which I
should prefer the one over the other (or any other platform). An example for
a reason could be if it will make it easier for me t
Try this:
read.table('don.5.clusters.txt', header = TRUE, comment.char = '', quote =
'')
On Tue, Oct 26, 2010 at 1:15 PM, Donald Braman wrote:
> That's one of the things I tried, but which didn't work. I get the
> following
> error when I do that:
>
> Error in read.table(file = "don.5.clusters.
On Tue, Oct 26, 2010 at 8:38 AM, Jonathan P Daily wrote:
> I think that this would be possible if you save the graph to a scalable
> format.
This is true to an extent. I have not checked on postscript or xfig,
but at least for PDF, even though you can "zoom"/blow the picture up,
you still have f
On Oct 26, 2010, at 8:30 AM, Alaios wrote:
in a simple plot. When i do plot is it possible to zoom in or out or
this is not
possible at all?
Zoom? Do you mean restrict the region plotted to specific ranges?
xlim and ylim arguments provide that facility.
--
David.
Best Regards
Alex
Hi Dimitri,
The lme function is not in the lme4 package, so there is some
confusion there. But you can use weights with the lmer function in
lme4. ?lmer tells you that weights are specified the same way as in
the lm function, and refers you to ?lm for details.
HTH,
Ista
On Tue, Oct 26, 2010 at 11
I think that this would be possible if you save the graph to a scalable
format.
Try looking into:
?postscript
?xfig
?pdf
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a
Hi Alex,
After you have created the plot I do not know of a way to zoom (in
base graphics), but you can always use the xlim and ylim arguments to
focus on a particular area (or effectively zoom) when you are create
the plot. For instance,
plot(x = 1:10, y = 1:10)
plot(x = 1:10, y = 1:10, xlim
On Oct 26, 2010, at 8:08 AM, Small Sandy (NHS Greater Glasgow & Clyde)
wrote:
Hi
I need some help getting results from multiple linear models into a
dataframe.
Let me explain the problem.
I have a dataframe with ejection fraction results measured over a
number of quartiles and grouped
in a simple plot. When i do plot is it possible to zoom in or out or this is
not
possible at all?
Best Regards
Alex
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-hel
Hello!
I am sorry if it's a naive/wrong question. But can one run a
regression with weights using lme?
Thank you!
--
Dimitri Liakhovitski
Ninah Consulting
www.ninah.com
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On Oct 26, 2010, at 7:50 AM, Viechtbauer Wolfgang (STAT) wrote:
The constraint b1=b2 in a model such as b0 + b1 x1 + b2 x2 + b3 x3
implies that b0 + b1 (x1 + x2) + b3 x3, so just add x1 and x2 (call
this x12) and fit the model b0 + b1 x12 + b3 x3 and you have imposed
the constraint that b1
That's one of the things I tried, but which didn't work. I get the following
error when I do that:
Error in read.table(file = "don.5.clusters.txt", header = TRUE, comment.char
= "", :
more columns than column names
If I remove the hashes by other means, I don't get that error.
On Tue, Oct 26
Right, I forgot to mention to use header = T.
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a room when its empty? Does the room,
the thing itself have purpose? Or do we
Hi
I need some help getting results from multiple linear models into a dataframe.
Let me explain the problem.
I have a dataframe with ejection fraction results measured over a number of
quartiles and grouped by base_study.
My dataframe (800 different base_studies) looks like
> afvtprelvefs
base
Try this:
Lines <- "SampleIDDisease
E-CBIL-28-raw-cel-1435145228.cel1
E-CBIL-28-raw-cel-1435145451.cel2
E-CBIL-28-raw-cel-1435145479.cel2
E-CBIL-28-raw-cel-1435145132.cel3
E-CBIL-28-raw-cel-1435145417.cel3
E-CBIL-28-raw-cel-1435145301.cel2
E-
"A call to read.table(..., sep = "", ...) reads in any length of whitespace
as the delimiter. On your sample text it read in a 2 column dataframe. "
Thanks! That works for the file, but when I enter in my next line of code it
doesn't work. I'm not sure if this is the problem, or if the next line
If I try that I get this:
Error in scan(file, what = "", sep = sep, quote = quote, nlines = 1, quiet =
TRUE, :
invalid 'sep' value: must be one byte
?
--
View this message in context:
http://r.789695.n4.nabble.com/Reading-in-a-tab-delimitated-file-tp3013620p3013771.html
Sent from the R he
Caveats and disclaimers:
I am quite happy to undertake self-teaching if directed to a relevant
prior posting and welcome such
direction. I have programming and statistical training/experience which
I would characterize as Masters level.
Thank you for reading and replying to this post. It is very
I have a vector of monthly log asset price returns. I would like use the boot
package to sample one-year returns, and then calculate confidence intervals on
the loss distribution only. More concretely, I would like to say something like
"99% of LOSSES (not RETURNS) are above cutoff X."
If the v
1 - 100 of 143 matches
Mail list logo
