Amitoj S. Chopra wrote:
>
> I am trying to write a program that uses R and takes a pdb file, and
> converts it to a pqr file. This task is simple generally, using the
> website, http://pdb2pqr-1.wustl.edu/pdb2pqr/. How do you use R to input a
> pdb file (that is on hand) into the upload pdb file
RagingJim wrote:
>
>
> ... it will ask the user for the name of the csv file it wants to open,
> and then prompt for what the user wants to save it as.
>
http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/FileOpenSave.html
RagingJim wrote:
>
>
> Secondly, I need to create a running mean
So I have to use this table of min, max, and mean temps for certain
years http://www.stat.berkeley.edu/classes/s133/data/january.tab. I am supposed
to figure out which year had the hottest January and which had the coldest. But
I dont know how to!
Nick Manginelli
[[alternative
Dear R helpers
I have some variables say ABC, DEF, PQR, LMN and XYZ. I am choosing any three
varaibles at random at a time for my analysis and name these files as
input1.csv, input2.csv and input3.csv. So if I choose variables say ABC, DEF
and PQR, I am passing the specifications of these vari
Steve Murray-2 wrote:
>
> Dear all,
>
> I have a grid (data frame) dataset at 0.5 x 0.5 degrees spatial resolution
> (720 columns x 360 rows; regular spacing) and wish to coarsen this to a
> resolution of 2.5 x 2.5 degrees. A simple calculation which takes the mean
> of a block of points to for
G'day, I am new user to R, and have been thrown in the deep end with a
something my company want me to write.
my code is as follows:
kenttemp=read.csv("mnowak.11.1.csv")
rows=nrow(kenttemp)-5
kent=kenttemp[1:rows,] #have to remove the last 5 lines of the
graph as they interfere with
bluesky...@gmail.com wrote:
I got the following output when I try the code at the end.
str(attr(terms1, 'variables'))
language list(skips, Panel, Opening)
class(attr(terms1, 'variables'))
[1] "call"
I checked ?call, where 'call' is a function. Could somebody let me
know what a 'call' class
Sunita Patil wrote:
Hello Peter
Sorry I forgot to paste the error for tittle it was as follows
*Error in title(main = "Year Vs StudentsPassed") : plot.new has not been
called yet*
Well, that makes sense. You can't put a title on a
plot that doesn't exist.
-Peter
Regards
Our Thoughts h
On Tue, Feb 9, 2010 at 6:09 PM, John Fox wrote:
> Dear bluesky315,
>
> There are several ways in R to determine regressors associated with factors.
> One way is to set the global contrasts option. To get "deviation"
> regressors, use options(contrasts=c("contr.sum", "contr.poly")), and see
> ?opti
Hi Hadley,
Thank you for taking the time to help me with this - I've constructed
the following example to illustrate my problem:
require(ggplot2)
data.set <- data.frame(
Time = c(rep(1, 4),rep(2, 4), rep(3, 4), rep(4, 4)),
Type = rep(c('a', 'b', 'c', 'd'), 4),
Value = c(10, 12, 14, 16, 14, 14
Hi Dennis,
Thank you for your response. No, NB is not a matrix, and I have no covariates.
Here's a very small sample of the data:
effect NB
-0.003200 1
-0.120800 3
-0.003200 2
-7.69 1
-1.442100 2
-0.000900 1
-0.014200 3
-5.015000 0
-0.0014
On 02/10/2010 01:20 AM, Karin Lagesen wrote:
I have a data set which I would like to plot as a set of concentric
circles. The data represent a count of the number of characteristics
shared by various elements - an example would look like this:
1 100
2 75
3 50
4 25
I.e. all four sets share 25 ch
I got the following output when I try the code at the end.
> str(attr(terms1, 'variables'))
language list(skips, Panel, Opening)
> class(attr(terms1, 'variables'))
[1] "call"
I checked ?call, where 'call' is a function. Could somebody let me
know what a 'call' class is and what 'language' is?
#
Hello Peter
Sorry I forgot to paste the error for tittle it was as follows
*Error in title(main = "Year Vs StudentsPassed") : plot.new has not been
called yet*
Regards
Our Thoughts have the Power to Change our Destiny.
Sunita
Sent from Pune, MH, India
On Tue, Feb 9, 2010 at 11:22 PM, Peter E
Dear Sir/Madam,
May I know if there is any function that estimates the ARCH or GARCH
models based on truncated likelihood?
Thanks,
Helen
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting gui
It depends on the type of file and your system. 'count.fields()' is
impractical for large files because it generates a matrix with the same
number of dimensions as the file. It would be easier to use scan() with the
delimiter argument set up to read to the end of line marker, "\n" I believe,
and th
One possibility I can see is to replace - by NA and use mean with
na.rm=TRUE.
--- On Wed, 10/2/10, Steve Murray wrote:
> From: Steve Murray
> Subject: [R] Resampling a grid to coarsen its resolution
> To: r-help@r-project.org
> Received: Wednesday, 10 February, 2010, 3:20 AM
>
> Dear all,
Use "\\^" instead.
For example:
df<-data.frame(c(1,2,4),c(4,3,2))
names(df)<-c("Amt","Resp")
df.form<-formula(Resp~0+Amt+I(Amt^2), data=df)
strsplit(as.character(df.form)[3], "\\^", perl=T)
should work just fine.
Sincerely,
KeithC.
-Original Message-
From: Gabor Grothendieck [mailto:ggr
Hi R-users,
I have this code here:
library(numDeriv)
fprime <- function(z)
{ alp <- 2.0165;
rho <- 0.868;
# simplified expressions
a <- alp-0.5
c1 <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)
c2 <- sqrt(rho)/(1-rho)
t1 <- exp(-z/(1-rho))
t2 <- (z/(2*c2))^a
bes1
Can one write a zip file from R?
I want to create a file with a name like "dat.zip", being a zip
file containing "dat.csv". I can create "dat.csv", then call
"system('zip -r9 dat.zip dat.csv')". Is there a better way?
I can use "gzfile" to write a gz file, but I don't kno
Dear Sir/Madam,
May I know if there is any function that estimates the ARCH or GARCH
models based on truncated likelihood?
Thanks,
Helen
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting gu
Thanks so much for your help. I am realizing that I may be
over-complicating things for myself. I have learned a ton about creating
methods, but I feel like I am trying to reinvent the data.frame class.
Basically, I am trying to create a data.frame type object where I can
enforce the header names
Hi,
Is your code reproducible?
I can't even find the package "adapt" on the CRAN repository. I am not sure
what exactlt happened to that package, but do remember seeing something about
it relatively recently in R-help. Are you using an older version of adapt?
Ravi.
_
Not exactly sure what you are asking: is this it
> DF <- data.frame(length=c(1,2,3,4),
+ Tanks=c("a3","a4","c4","h4"))
> DF
length Tanks
1 1a3
2 2a4
3 3c4
4 4h4
> DF$length[DF$Tanks == 'h4'] - DF$length[DF$Tanks == 'a3']
[1] 3
>
On Tue, Feb 9, 2010 at 9:37 P
Tena koe Keith
?which
perhaps. As in:
which(Tanks=='h4')-which(Tanks=='a3')
HTH ...
Peter Alspach
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Marlin Keith Cox
> Sent: Wednesday, 10 February 2010 3:37 p.m.
> To: r-hel
OK, this is very elementary, but I need help. I have looked in Verzani,
past postings etc.
Problem: I need to subtract the "length" date between "h4" and "a3" #which
would be 4-1
I would rather not convert the two columns into four columns (with headings
being "a3","a4","c4","h4").
DF <- data
Note that the first set of coefficients minus the second set of coefficients
is constant (12.78262)
--
View this message in context:
http://n4.nabble.com/linear-predictors-and-survreg-function-tp1474388p1475332.html
Sent from the R help mailing list archive at Nabble.com.
__
WIll this do it for you:
> x <- read.table(textConnection("ColA ColB
+ 10
+ 30
+ 21
+ 20
+ 10
+ 40
+ 11
+ 91
+ 30
+ 50
+ 21"), header=TRUE)
> closeAllConnections()
> x.s <- split(x, cumsum(x$ColB))
> x.l <- do
Dear R Help:
I am trying to sum data from one column in a dataframe based on a value in
another. I do not know how to do this easily in R.
For example:
Col A Col B
10
30
21
20
10
40
11
91
30
50
21
I would l
Hello All,
I have a dataset with a continuous response variable and an ordered factor
predictor. I am very interested in using the package OrdFacReg to run my
analysis, but I am having a difficult time deciphering the code and making it
work for my dataset. Given that this is a new package, I
> I am trying to figure out how to control table width and alignment on the
> page for a table generated by odfTable. Based on reading odfWeave
> documentation (including formattingOut.odt), here is how I manipulate the
> styles:
>
st = getStyleDefs()
# modify the table style
tab = getStyles()$ta
Tena koe Nilza
I think merge() will do what you want. For example:
> completeDates
V1
1 01-12-09
2 02-12-09
3 03-12-09
4 04-12-09
5 05-12-09
6 06-12-09
7 07-12-09
8 08-12-09
> testData
V1 V2
1 01-12-09 26.8
2 03-12-09 27.9
3 04-12-09 25.6
4 05-12-09 20.8
5 08-12-09 20.8
> merge
On Tue, Feb 9, 2010 at 6:46 PM, Martin Batholdy wrote:
> hi,
>
> I have a vector full of strings like;
>
>
> xy_100_ab xy_101_ab xy_102_ab xy_103_ab
>
>
> I want to seperate each string in three pieces and the separator should be
> the "_"
>
> at the end I want a data.frame like
Hi:
# Create a sequence of consecutive dates and assign to a data frame
dates <- data.frame(day = seq(as.Date('2009-12-01'), by = 'days', length =
10))
# Let's call your data frame (not vector) 'df'; to make sure that things are
# compatible,
names(df) <- c('day', 'y')
df$day <- as.Date(df$day)
Dear all,
I am trying to identify the days absent in a vector so as to be able to
display missing data in the graphics in .It should be done automatically
since the graphics will be generated for different periods of the year.
I have been trying to do a function below but it did not work. I re
I'm interested in using R's plotting capabilities to try to generate a graph
showing the relationship/pairing frequency of words appearing in a block of
unstructured text. I don't have a specific algorithm or approach in mind, just
looking to portray text in an interesting fashion. The output I
Dear bluesky315,
There are several ways in R to determine regressors associated with factors.
One way is to set the global contrasts option. To get "deviation"
regressors, use options(contrasts=c("contr.sum", "contr.poly")), and see
?options and ?contrasts for details. Also see Section 11.1.1 of t
On Feb 9, 2010, at 6:46 PM, Martin Batholdy wrote:
hi,
I have a vector full of strings like;
xy_100_ab xy_101_ab xy_102_ab xy_103_ab
I want to seperate each string in three pieces and the separator
should be the "_"
at the end I want a data.frame like:
column1 column
Hi Martin,
Here is a sugestion:
string <- c("xy_100_ab", "xy_101_ab","xy_102_ab","xy_103_ab")
out <- data.frame( do.call( rbind, strsplit( string, '_' ) ) )
names(out) <- paste('column',1:3,sep="")
out
HTH,
Jorge
On Tue, Feb 9, 2010 at 6:46 PM, Martin Batholdy <> wrote:
> hi,
>
> I have a vect
On Feb 9, 2010, at 6:33 PM, bluesky...@gmail.com wrote:
The model matrix for the code at the end the email is shown below.
Since the model matrix doesn't have -1, I think that it is made of
dummy regressors rather than deviation regressors. I'm wondering how
to make a model matrix using deviati
I'd like to superimpose a ksmooth() onto a barplot().
My data is:
> d
2009-06-20 2009-06-21 2009-06-22 2009-06-23 2009-06-24
2009-06-25 2009-06-26 2009-06-27 2009-06-28 2009-06-29 2009-06-30
2009-07-01 2009-07-02
Same Breed (B) 12.64 21.08 13.52 12.51
hi,
I have a vector full of strings like;
xy_100_ab xy_101_ab xy_102_ab xy_103_ab
I want to seperate each string in three pieces and the separator should be the
"_"
at the end I want a data.frame like:
column1 column2 column3
xy 100 ab
xy
The model matrix for the code at the end the email is shown below.
Since the model matrix doesn't have -1, I think that it is made of
dummy regressors rather than deviation regressors. I'm wondering how
to make a model matrix using deviation regressors. Could somebody let
me know?
> model.matrix(a
Hi Liam,
Yes, that's what that code should do. Could you please send a small
reproducible example?
Hadley
On Tue, Feb 9, 2010 at 4:21 PM, Liam Blanckenberg
wrote:
> Hadley,
>
> Thanks for the pointing that error out to me. Unfortunately, revising
> my code to:
>
> p <- ggplot2(~, aes(x = ~, y
I received also this message that did not show up here:
This sounds like a problem in (linear) mixed effects models. Look into
packages
nlme (with the lme() function) and lme4 (function lmer()) for starters. What
you
have falls under the heading of longitudinal or repeated measures data. Both
pack
Thanks
hadley wickham wrote:
Hi Paul,
That's a bug in the current version of ggplot. I'm working on update
for later this week.
Hadley
On Mon, Feb 8, 2010 at 5:56 PM, Paul Sutcliffe wrote:
In ggplot2 how do you justify the legend text ?
In the example below the opts(legend.text = them
Hello all,
I am trying to minimize a function which contains a double integral, using
"nlminb" for the minimization and "adapt" for the integral. The integral is
over two variables (thita and radiusb)
and the 3 free parameters I want to derive from the minimization are
counts0, index and radius_
Hello,
I am experiencing problems viewing the help files in R. The easiest
way to describe this is to just post an image. You can see what I'm
talking about here: http://www.flickr.com/photos/matthewbk/4344722268/
If I copy and paste the "weird" characters, they show up like they
should.
Hadley,
Thanks for the pointing that error out to me. Unfortunately, revising
my code to:
p <- ggplot2(~, aes(x = ~, y = ~, colour = Type)) +
geom_line(aes(colour = Type), position = 'stack')
still does not generate what I'm after. I'm essentially after a line
plot where each 'Type' ('series' in
## Continuing with the same example, I now show the full set of
contrasts and their
## associated regression coefficients. Remember that regression coefficients in
## models with factors are meaningful only when the dummy variables
are also displayed.
## construct full set of dummy variables to d
Thanks Prof. Ripley. I had forgotten that the OP had specified internet2. I
don't use it, so my suggestions only apply when it is turned off.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From:
You can use the history function to see a list of previous commands (but it
does not automatically run them).
Do
> history(pat='plot')
To get a list of all the commands that you ran that included the string 'plot'
somewhere in it. You can then copy and paste to rerun (under the windows GUI
(
Version 3.20 of the survey package is percolating through CRAN
The major additions since the last announcement on this list (3.18, in
September) are
- database-backed designs can now use replicate weights
- some multivariate statistics: principal components, factor analysis.
The NEW
It is not clear what exactly you want, but here are some ideas.
Lattice and ggplot2 allow you to save graph objects that can then be
modified/printed/plotted, there is information in there on types of axes etc.
For basic graphics (and others) you can use the history command to see recent
plotti
On Tue, Feb 9, 2010 at 4:38 AM, Demirtas, Hakan wrote:
> Does lmer do three-level mixed-effects models?
Yes.
> What forms of outcome
> variables can it handle (continuous, survival, binary)? I'd appreciate any
> help.
Continuous outcomes for lmer. glmer can handle binary outcomes.
__
Googling for: contrasts for interaction terms produces useful results. For
example:
http://www.ats.ucla.edu/stat/spss/faq/test_subcommand_mixed_mod.htm
Rich's example gave you an interpretation of a three-way interaction using
contrasts. Maybe a bit hidden, but they are there. For example, excl
On 02/08/2010 11:05 PM, Jiiindo wrote:
>
> Thanks for your answer,
> I can run it with an application Java normal. In my web-application, when i
> need, i start R et calling R function (new REvaluator).eval("1+1") to
> calcule an expression "1+1"example.
> And output of your code:
>
> Loading RIn
As others have said, z[z[, 1] < 0, ] does it. Just in case
you're wondering why your subset command won't work,
str() is your friend (as is so often the case):
str(z)
str(as.data.frame(z)) ## (I don't think that R has 'as.data.set')
So z is a matrix with column *names* "x1", etc;
as.data.frame(
This is useful. But it is not exactly what I am looking for. I can use
'contrast' to infer linear combination factor level means as discussed
in Statistical Inference by Casella and Berger. What is the equivalent
'contrast' for interaction terms? Is there a book that discuss this?
On Tue, Feb 9, 2
Hi:
> z = matrix(rnorm(n*m),n,m,dimnames =list(rnames,cnames))
> z
x1 x2 x3 x4 x5
1 -0.3942900 -0.61202639 -1.804958629 -0.1351786 -1.27659221
2 -0.0593134 0.34111969 1.465554862 1.1780870 -0.57326541
3 1.1000254 -1.12936310 0.153253338 -1.523
Try this:
z[z[,1] < 0,]
On Tue, Feb 9, 2010 at 6:12 PM, DonDiego wrote:
>
> Hi,
>
> I have a matrix of data values like the example bellow. I would like to
> extract a subset of the matrix for the values where the first column is
> negative. I am using the subset function. However, I am getting
On Tue, 9 Feb 2010, Greg Snow wrote:
Timing can be an issue. You need to set the environment variables before
accessing the internet for the 1st time in each R session. If you get the
error, then try setting the variables, it will not work.
For my situation I have http_proxy and http_proxy_
On Feb 10, 5:21 am, mnstn wrote:
> Hello All,
> In bash, to re-execute a command, say, 'gnuplot plot.gnu', one can type
> !gnu. Is there a similar feature in the R console? For example:
> plot(g1$mean,g2$mean)
> .
> .
> .
> Lines, many lines of commands
> ->Here I want to repeat the previous plo
Hi,
I have a matrix of data values like the example bellow. I would like to
extract a subset of the matrix for the values where the first column is
negative. I am using the subset function. However, I am getting an error
message that the conditional variable doe snot exist. For some reason, the
Timing can be an issue. You need to set the environment variables before
accessing the internet for the 1st time in each R session. If you get the
error, then try setting the variables, it will not work.
For my situation I have http_proxy and http_proxy_user set for the entire
computer (in wi
What do you mean by 'reset'? Why not put your processing in a
function that you call so that each time you invoke it, any variables
local to the function start in an 'initialized' state (depends on how
you are using them). If you need them to have some specific value,
then put in the code necessa
Hi ,
I have several commands that I want to run, I was wondering if there is an
easy way to reset my variables before the start of each run
for ( i in 1:200){
reset the variables
run the commands
}
My solution right now is to set all my vectors to 0 and my data frames to
NULL after the start o
i also having the problem with R .
R version 2.10.0 (2009-10-26)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()
Thanks David,
I guess I have to settle with the upper-arrow thing I have been doing all
along. My figures do not return any value so ggplot, lattice are out.
--
View this message in context:
http://n4.nabble.com/Re-execute-previous-command-tp1474629p1474813.html
Sent from the R help mailing list
## Artificial data with all interactions significant.
## The interaction2wt plot shows all main effects and all pairwise
## interactions. We see in the "Y ~ A|B" panel (or in the
## interaction.plot) that Y goes uphill for levels 1 and 2 of B and
## goes down and then up for level 3 of B. This i
You might do better posting this on the R-sig-mixed-models list.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ista Zahn
Sent: Tuesday, February 09, 2010 9:19 AM
To: Weber, Sam
Cc: r-h
Here is a simple 3-step solution:
1. type ?barplot
2. find the section labelled 'Arguments'
3. carefully read what each argument means/does
Voila, problem(s) solved.
(It was good to include some data, but not so good
to say that title(...) is giving "an error". No need
to be secretive about the
Hello all; My question is part statistical and part R.
I have performed model II major axis regression in R using both the smatr()
and lmodel2() packages, but neither offers an option to statistically weight
my regression. I have a vector of weights which I would like to apply to
each of my regr
Classification trees (rpart package), but know too little about your
problem...
Best, Andrej
On Feb 9, 5:50 pm, Robert Lonsinger wrote:
> Can anybody recommend a nonparametric alternative to linear discriminant
> function that may be available as a package in R?
>
> Cheers,
> Rob
>
> [[a
Jeff Harring wrote:
>
> I am looking for R code to be able to fit a linear-linear piecewise
> model with person-specific changepoint. I have searched the web, but
> have not been able to locate any code.
>
> ..
> I would appreciate any advice on this code or guidance as to where to
> search
Simone Vincenzi-2 wrote:
>
> I have a statistical problem with the comparison of short time-series,
> representing densities of fish in different streams. For each stream (6 in
> total, here below showed only part of the dataset) I have 8 years of
> density
> data, as follows:
>
>
>year
Hello
(this might be a very simple question)
My data is as follows (table name is student)
YearStudentsPassed
1 2000300
2 2001360
3 2002450
4 2003450
5 2004270
6 2005280
7 2006400
8 2007270
I want to plot a barplot so for this
If I use barplot(s
On 09/02/2010 12:18 PM, bluesky...@gmail.com wrote:
In ?eval, it says
... If ‘envir’ is not
specified, then the default is ‘parent.frame()’ (the environment
where the call to ‘eval’ was made).
I tried the following example with "eval(expr)" and "eval(expr,
parent.frame())" in f().
Dear everyone,
I'm coding the Horowitz-Spokoiny (2001) test [1], and I would be very
grateful or some advice regarding the Kernel density (apologies
beforehand if my terminology is not fully correct). I have looked into
ksmooth and npreg, but with no success.
Given a (n x p) matrix of covariate
On Feb 9, 2010, at 12:13 PM, David Winsemius wrote:
On Feb 9, 2010, at 11:21 AM, mnstn wrote:
Hello All,
In bash, to re-execute a command, say, 'gnuplot plot.gnu', one can
type
!gnu. Is there a similar feature in the R console? For example:
plot(g1$mean,g2$mean)
In most system typeing
In ?eval, it says
... If ‘envir’ is not
specified, then the default is ‘parent.frame()’ (the environment
where the call to ‘eval’ was made).
I tried the following example with "eval(expr)" and "eval(expr,
parent.frame())" in f(). The results are different, which are not
consistent w
Dear R and odfWeave users,
I am trying to figure out how to control table width and alignment on the
page for a table generated by odfTable. Based on reading odfWeave
documentation (including formattingOut.odt), here is how I manipulate the
styles:
st = getStyleDefs()
# modify the table style
tab
Hi Sam,
Good question. I originally guessed that the "simple effect" (I know
some people on this list don't seem to care for that term, but it's
always made sense to me) coefficients were in the same direction, such
that the effect if Origin at Treat=hot was significantly different
from zero, but n
On Feb 9, 2010, at 11:21 AM, mnstn wrote:
Hello All,
In bash, to re-execute a command, say, 'gnuplot plot.gnu', one can
type
!gnu. Is there a similar feature in the R console? For example:
plot(g1$mean,g2$mean)
In most system typeing will prompt the console to echo the
last command (a
On Feb 9, 2010, at 11:24 AM, Alex Levitchi wrote:
Hello
I am recently began to work with R, so I am not so experienced.
But anyway I cannot find a clear way to process my dataframe which
is a bigger one.
It shows similar to this
name=c("A","B","C","B","C","C","C","B","C")
nicknames=c("A1",
try this:
> x <- read.table(textConnection("name nicknames value
+ 1 A A1 4
+ 2 B B1 5
+ 3 C C1 9
+ 4 B B2 2
+ 5 C C2 7
+ 6 C C3 6
+ 7 C C4 3
+ 8 B B3 6
+ 9 C C5 7"), header=TRUE)
> closeAllConnections()
> result <- do.call(rbind, lapply(split(x, x$name), function(.name){
+ data.frame(name=.na
Dear R-list,
I have a statistical problem with the comparison of short time-series,
representing densities of fish in different streams. For each stream (6 in
total, here below showed only part of the dataset) I have 8 years of density
data, as follows:
year density stream
1 2000 0.51 s
Hello All,
In bash, to re-execute a command, say, 'gnuplot plot.gnu', one can type
!gnu. Is there a similar feature in the R console? For example:
plot(g1$mean,g2$mean)
.
.
.
Lines, many lines of commands
->Here I want to repeat the previous plot command.
!plot and plot(+TAB) do not work.
This q
Can anybody recommend a nonparametric alternative to linear discriminant
function that may be available as a package in R?
Cheers,
Rob
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listin
Hello all,
Here is my solution, in case someone else needs it.
I have a dataframe consisting of two columns.
col1<-factor(c("a","a","b","b","c","c"))
col2<-factor(c("a","b","c","d","e","f"))
somedf<-data.frame(col1,col2)
somedf
col1 col2
1 a d
2 a e
3 b f
Hello
I am recently began to work with R, so I am not so experienced.
But anyway I cannot find a clear way to process my dataframe which is a bigger
one.
It shows similar to this
> name=c("A","B","C","B","C","C","C","B","C")
> nicknames=c("A1","B1","C1","B2","C2","C3","C4","B3","C5")
> valu
Dear list users,
I was wondering if there is any way to know, for an open graphical peripheric,
the instruction that was used to call it. For example, if I create a plot using
plot(runif(10))
is there any way to have the call returned to me? That would be especially
usefull to know it there ar
Dear all,
I have a grid (data frame) dataset at 0.5 x 0.5 degrees spatial resolution (720
columns x 360 rows; regular spacing) and wish to coarsen this to a resolution
of 2.5 x 2.5 degrees. A simple calculation which takes the mean of a block of
points to form the regridded values would do the
If it is not standardized, I have the following script that I use to
calculate the r^2 for each variable of a model. I do wonder, however,
if there is a way to calculate a model-wide r^2 like index.
-Jarrett
##
#returns the r^2 for endogenous variables and error path coefficients
#
# Last up
I have difficulties in interpreting high order interaction terms in
high-way ANOVA.
According to Introductory Statistics with R by Peter Dalgaard (Section 12.5),
"The exact definition of the interaction terms and the interpretation of their
associated regression coefficients can be elusive. Some
Thank you for your rapid answer!
--
View this message in context:
http://n4.nabble.com/step-and-glmer-tp1474390p1474580.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/lis
thank you all for your help, the get() function worked perfectly.
--
View this message in context:
http://n4.nabble.com/performing-the-same-commands-on-two-different-data-sets-tp1474452p1474515.html
Sent from the R help mailing list archive at Nabble.com.
__
I am looking for R code to be able to fit a linear-linear piecewise
model with person-specific changepoint. I have searched the web, but
have not been able to locate any code.
Below is my attempt at some code:
chgpt = function(a1,a2,a3,gam,wave){
yht=numeric(10)
y1=(wave <= gam)*(a1+(a2*wave))
On Tue, Feb 9, 2010 at 2:20 PM, Karin Lagesen wrote:
> I have a data set which I would like to plot as a set of concentric
> circles. The data represent a count of the number of characteristics
> shared by various elements - an example would look like this:
>
> 1 100
> 2 75
> 3 50
> 4 25
>
> I.e.
Hi Liam,
Your syntax is a little off. You want:
p <- ggplot2(~, aes(x = ~, y = ~, colour = Type)) +
geom_area(aes(fill = Type), position = 'stack')
Position isn't an aesthetic.
Hadley
On Sun, Feb 7, 2010 at 10:40 PM, Liam Blanckenberg
wrote:
> Hi all,
>
> I have been hunting around for hou
1 - 100 of 132 matches
Mail list logo