That's the only one I pull up with a quick search, but it seems pretty
good -- what do you need from it?
Best,
Michael
On Fri, May 18, 2012 at 7:16 PM, JWDougherty j...@surewest.net wrote:
Can someone point me to an implementation of the runs up and runs
down test, or does such beast exist in
On May 19, 2012, at 01:45 , Noah Silverman wrote:
Hello,
I need to design a fairly simple front-end for someone to use an R script
system that I've built. My thought was to just use the text based menus
available in the base R package, perhaps in some kind of loop.
How have other
JWDougherty jwd at surewest.net writes:
Can someone point me to an implementation of the runs up and runs
down test, or does such beast exist in R? From web searches the
runs up and runs down test is commonly used for testing pseudo-random
number generators and in simulations. John C. Davis
Hi all,
I am a beginner R user and need some help with a simple loop function.
Currently, I have seven datasets (TOWER1,TOWER2...TOWER7) that are all in
the same format (same # of col and headers). I am trying to add a new column
(factor) to each dataset that simply identifies the dataset.
Hi R-listers,
Q1) What is the difference between the scatterplot and plot function?
Q2) I am able to make a graph with the scatterplot function:
scatterplot(DevelopIndex ~ Veg,
+ data = Turtle,
+ xlab = Vegetation border (m),
+ ylab = Embryonic development
On May 18, 2012, at 16:15 , David Costantini wrote:
Sorry about that.
Basically I have two fixed factors (experimental group and sex) and a random
factor
(brood). I need the random factor to control for pseudoreplication because I
have siblings,
which are pseudoreplicates because share
Dear Peter
thank you for your suggestions.
the design is balanced in terms of birds per group or sex, but not for brood.
I can't do averages because (i) siblings from a same brood are randomly
allocated to
different experimental groups, (ii) this would drastically reduce the power of
analysis.
Thanks so much! I thought R places NA for missing values. I'll have to
read up on it more. Thanks again!
On May 18, 2:23 pm, Marc Schwartz marc_schwa...@me.com wrote:
On May 18, 2012, at 2:26 PM, Giggles wrote:
I am a newbie and can't figure out how to recode a numeric value. In
my data
Thank you very much!!
On May 18, 2:22 pm, David Winsemius dwinsem...@comcast.net wrote:
On May 18, 2012, at 3:26 PM, Giggles wrote:
I am a newbie and can't figure out how to recode a numeric value. In
my data (pharm311), I have a column called explain and I need to
find all the 6's and
On Sat, May 19, 2012 at 1:18 AM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
I had such good luck with my previous question to r-help, (a few minutes
ago) that I thought I would try again with the following query:
Suppose I have
xNm - gamma
I would like to be able to do
Hi,
I want to do an interactive plot in R:
I have a time series and I can plot it in the normal way:
plot(modifieddate,timeseries,type=l,xlab=Date,ylab=values, main=title)
(modifieddate is the values for the date on the x axis)
The aim now is, that I want to do a Latex Presentation, where if
Hello everyone,
I want to get a 1km by lkm grid raster image using my csv data. If I call
latitude=a, longitude=b and preciptation=c.
a-(1,2,3,4,5)
b-(6,7,8,9,10)
c-(10,20, 30,40, 50)
Then I found an example in r help which goes like
pts = read.table(file.csv,..)
library(sp)
library(rgdal)
You might want to give the packages animation a try.
On May 19, 2012, at 11:47 AM, c...@mail.tu-berlin.de wrote:
Hi,
I want to do an interactive plot in R:
I have a time series and I can plot it in the normal way:
plot(modifieddate,timeseries,type=l,xlab=Date,ylab=values, main=title)
First, what was going wrong:
#this creates a character vector, not an object containing your data sets
TOWERS - c(TOWER1,TOWER2,TOWER3,TOWER4,TOWER5,TOWER6,TOWER7)
for(i in 1:7){
# this creates a variable named TOWER.i (not TOWER.1, TOWER.2,
..., TOWER.7)
# which will get overwritten
provided you get the call to read.table (or perhaps read.csv) right
and presuming that file contains only the image data, you should
be able to say:
r - raster(as.matrix(read.csv(file.csv)))
extent(r) - extent(xmin, xmax, ymin, ymax)
and not worry about the projection (if it is plain old decimal
I presume you mean car::scatterplot
A1) plot is a generic function for plotting 'things' and scatterplot
is a specific tool provided by that library
A2) from the documentation one infers that
1) the straight line is a regression line
2) the other lines are produced using loess
On May 18, 2012, at 20:19 , Patrick Callier wrote:
Hi all,
I am running 64-bit R 2.15.0 on windows 7. I am trying to use read.delim
to read from a file that has 2-byte unicode (CJK) characters.
Here is an example of the data (it is tab-delimited if that gets messed up):
HITId HITTypeId
This already in your question:
+init=epsg:4326
That is equivalent to +proj=longlat +ellps=WGS84 just looked up via the
EPSG code. A good site to explore this topic is spatialreference.org
On Saturday, May 19, 2012, Belay Gebregiorgis wrote:
Hello everyone,
I want to get a 1km by lkm grid
mh, I have read about the package before, but it does not help me a
much I dont know how to do this?
Zitat von Erich Neuwirth erich.neuwi...@univie.ac.at:
You might want to give the packages animation a try.
On May 19, 2012, at 11:47 AM, c...@mail.tu-berlin.de wrote:
Hi,
I want to do an
Hi R users,
I have a dataset with multiple variables and i'm trying to weigh average
depths for fish species per year by their abundance (CPUE. I have tried the
weighted.mean function but as i have two columns for the x value the lenghts
differ with the w (CPUE). How do I solve this problem?
So
Thanks Simon. I have other vectors in the CSV file so after naming my
latitude and longitude x and y, I tried this. The first two lines run fine
but I get error in the second line.
r - raster(as.matrix(F1))
extent(r) - extent(xmin, xmax, ymin, ymax)
F1-data.frame(x, y,z)
Error in function
I am using the following:
library(RODBC)
chan = odbcConnectExcel(rats-lda)
rats.lda = sqlFetch(chan, data)
close(chan)
And getting the following error message:
library(RODBC)
Error in library(RODBC) : there is no package called ‘RODBC’
chan = odbcConnectExcel(rats-lda)
Error: could not find
Hi All,
I have a question regarding contingency tables. I would like to calculate the
mean and standard deviation of a continuous variable from my own dataset based
on the percentages of a contingency table I obtained from a scientific article.
dataset- data.frame(cbind(case=rep(0:1,5),
i have time series as
1.3578511
0.5119648
1.3189847
0.9214787
1.2272616
4.9167998
1.2272616
1.2272616
0.8854192
2.3386331
1.6132899
0.2030302
0.8426226
1.2277843
NA
1.3189847
1.3578511
0.8530141
2.3386331
1.0541099
0.7747481
0.5764672
1.3189847
1.2160533
1.2272616
0.6715839
0.9651803
1.6132899
On 05/19/2012 07:44 AM, SteveQ wrote:
I am using the following:
library(RODBC)
chan = odbcConnectExcel(rats-lda)
rats.lda = sqlFetch(chan, data)
close(chan)
And getting the following error message:
library(RODBC)
Error in library(RODBC) : there is no package called ‘RODBC’
The above error
David,
See corr.test in the psych package.
Bill
On May 18, 2012, at 9:35 AM, David Studer wrote:
Hi everybody!
Does anyone know how to obtain a inter-item-correlation-table
(with p-values or significance-levels)? (as SPSS does, either spearman or
pearson)
Repeatedly using cor.test()
Hello,
The error is that you are trying to use the number of rows of a
character string.
TOWERS[1], TOWERS[2], etc, are not data frames. Use a print statement
before the line that throws the error to check it.
Your problem can be solved along the lines of what follows.
Note that I've put all
Hello,
You are not using dates.
2011-01-01
[1] 2009
2011-07-01
[1] 2003
as.Date(2011-07-01)
Error in as.Date.numeric(2011 - 7 - 1) : 'origin' must be supplied
as.Date(2011-07-01)
[1] 2011-07-01
The error is because as.Date is seeing 2003 as a number, the number of
days since an
... and here is another incantation that may be informative.
xnm- as.name(gamma') ## This does the parsing
plot(0, xlab =bquote(.(xnm))
The initial puzzle is that if you just set
xnm - gamma
bquote will insert the string gamma rather than the symbol. After
all, that's what plotmath sees for
?install.packages
This functionality may also be available from a menu in an R GUI, e.g.
on WIndows. But as you failed to provide your system details -- as
requested by the posting guide -- we can't say.
HOWEVER ... something appears to be rotten in Denmark (with apologies
to Peter D.). IIRC,
Hello,
Post a data example, it will give us something to work on, and give you
better answers.
dput( head(allspecies, 20) )
Then copy the output of this and paste it in your post.
As for your problem, considering it's description, and untested (no
data), maybe
result - by(allspecies,
I don't think it's part of the standard install -- I certainly didn't
get it with Simon's [Mac] CRAN build and when I did download it it's
marked with priority NA -- so no need to worry there.
Best,
Michael
On Sat, May 19, 2012 at 10:36 AM, Bert Gunter gunter.ber...@gene.com wrote:
..and a final addendum:
xnm - quote(gamma)
## makes xnm the name gamma not the string gamma
plot(0,xlab = bquote( .(xnm))
-- Bert
On Sat, May 19, 2012 at 7:23 AM, Bert Gunter bgun...@gene.com wrote:
... and here is another incantation that may be informative.
xnm- as.name(gamma') ## This
library (zoo)
?na.approx
Note that you need to define an index (time base) to go along with your data,
but that could be as simple as a sequence of integers.
---
Jeff NewmillerThe .
On Sat, May 19, 2012 at 11:09 AM, Jeff Newmiller
jdnew...@dcn.davis.ca.us wrote:
library (zoo)
?na.approx
Note that you need to define an index (time base) to go along with your data,
but that could be as simple as a sequence of integers.
Actually na.approx in the zoo package has a default
Hello,
Thanks.
The revised 'by' I had posted works with me.
An alternative, more complicated, is using split/sapply.
sapply(split(allspecies, allspecies$Year), function(x)
weighted.mean(x$Depth, x$CPUE))
Rui Barradas
__
R-help@r-project.org
Hello,
Can someone help me on how to make a subset of my dataframe using two dates.
This is what I have tried, but no data is being excluded in the new frame.
six_months-subset(Two_years, vdate2011-01-01|vdate2011-07-01)
or all data is excluded
six_months-subset(Two_years,
Hi all,
I am a student and I ma writing my master thesis about oil and gasoline
options. In my master thesis I would like to estimate some parameter of the
oil market. I base my pricing models according to Schwartz (1997). In order
to get those data a I need to run a Kalman filter algorith with
This is the output that i got
dput( head(allspecies, 20) )
structure(list(Year = c(1980L, 1980L, 1980L, 1980L, 1980L, 1980L,
1980L, 1980L, 1980L, 1980L, 1980L, 1980L, 1981L, 1981L, 1981L,
1981L, 1981L, 1981L, 1981L, 1981L), Quarter = c(1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
Dear R Users,
I am facing a problem analyzing an incomplete block design with two
replicates. As you can see in the attached .xls file, the factor2 (6
levels) nested within factor1(two levels) nested within replicates all were
chosen as random effects in the statistical model (see below). Note
It works now :)
thank you!
--
View this message in context:
http://r.789695.n4.nabble.com/weighted-averages-for-two-variables-tp4630589p4630610.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Hey,
I'm doing my thesis and I need to produce a lot of plots. With 2 I have a
problem that I can't get past:
http://r.789695.n4.nabble.com/file/n4630611/question.png
The left one has lines instead of dots. All I want is a basic line (it's a
time series). So I used the type=l command, but that
R has four Kalman functions -- type apropos(Kalman) to see them --
then you can get info about any of them in particular by typing ?
funcname or running example(funcname) to see worked examples.
Beyond the scope of your question, you also might be interested in
this:
On Sat, May 19, 2012 at 9:31 AM, peter dalgaard pda...@gmail.com wrote:
Best practice is a bit contentious, but several people have found that the
tcltk package offers a path of low resistance.
Additionally check gWidgets.
Liviu
__
Hi Teo:
On May 19, 2012, at 6:57 AM, mala10af wrote:
Hi all,
I am a student and I ma writing my master thesis about oil and gasoline
options. In my master thesis I would like to estimate some parameter of the
oil market. I base my pricing models according to Schwartz (1997). In order
to
On 19.05.2012 15:34, Allen wrote:
Hello,
Can someone help me on how to make a subset of my dataframe using two dates.
This is what I have tried, but no data is being excluded in the new frame.
six_months-subset(Two_years, vdate2011-01-01|vdate2011-07-01)
Convert the *string* 2011-01-01 to
On 19.05.2012 17:40, bets wrote:
Hey,
I'm doing my thesis and I need to produce a lot of plots. With 2 I have a
problem that I can't get past:
http://r.789695.n4.nabble.com/file/n4630611/question.png
The left one has lines instead of dots. All I want is a basic line (it's a
time series). So
parse(text=paste(...)) works in simple cases but not in others. The
fortune about it is there because it is tempting to use but if you bury it
in a general purpose function it will cause problems when people
start using nonstandard names for variables. bquote(), substitute(),
call(), and
On Sat, May 19, 2012 at 11:07 AM, William Dunlap wdun...@tibco.com wrote:
parse(text=paste(...)) works in simple cases but not in others. The
fortune about it is there because it is tempting to use but if you bury it
in a general purpose function it will cause problems when people
start using
On May 19, 2012, at 7:46 AM, Ross, Stephanie wrote:
Hi All,
I have a question regarding contingency tables. I would like to
calculate the mean and standard deviation of a continuous variable
from my own dataset based on the percentages of a contingency table
I obtained from a scientific
Sorry, I hoped this was already clear enough.
Oké, here are the first rows of the relevant part of the data:
Datum week T.local
29/06/111 19.272420
5/07/11 2 19.305034
13/07/113 17.766046
20/07/114
On Sat, 19 May 2012 09:50:27 +0200
Hans W Borchers hwborch...@googlemail.com wrote:
JWDougherty jwd at surewest.net writes:
Can someone point me to an implementation of the runs up and runs
down test, or does such beast exist in R? From web searches the
runs up and runs down test is
On Sat, May 19, 2012 at 2:16 PM, Paul Johnson pauljoh...@gmail.com wrote:
On Sat, May 19, 2012 at 11:07 AM, William Dunlap wdun...@tibco.com wrote:
parse(text=paste(...)) works in simple cases but not in others. The
fortune about it is there because it is tempting to use but if you bury it
in
i use download.file() ,it works good for package source code
t-available.packages(type=source) #all packages(list) from repository in
source code .tar.gz format
length(t[,2]) #total packages-- 3796
#package source file download
for( i in seq_along(t[,1]) ) {
download.file(
On May 19, 2012, at 12:36 PM, bets wrote:
Sorry, I hoped this was already clear enough.
Below you offer two examples of console output without naming the
objects from which they came or the code you are attempting to us.
(Also in the all-to-typical manner of those posting through Nabble
You need to specify that you want a binary download - make sure to
read help(download.file), particular that of argument 'mode'.
Alternatively, use downloadFile() in the R.utils package which
download as binary by default.
/Henrik
On Sat, May 19, 2012 at 7:22 AM, sagarnikam123
Hi there,
I'm working on a for loop to perform some simple calculations on data.
However, I'm having trouble extracting components from my data frame using
names pulled off a list to reference headers. I'm relatively new to R, but
I've done some programming in various other languages, and have
On May 19, 2012, at 3:32 PM, Chris North wrote:
Hi there,
I'm working on a for loop to perform some simple calculations on data.
However, I'm having trouble extracting components from my data frame
using
names pulled off a list to reference headers. I'm relatively new to
R, but
I've done
That means it calls compiled C code so you'll need to download the
source code from CRAN to see it. It will be somewhere in the
rugarch/src/ directory, but you'll have to find it manually or with a
good editor/IDE -- you're out of R now.
The primitive() means that .C is itself a primitive
Dear R Helpers,
I am getting an error message from the try function that I don't
understand so I am hoping that someone can help.
I am scraping from web pages, but sometimes they disappear. When that
happens I need to control for it with some sort of function.
This web page is parsed without a
On May 19, 2012, at 5:00 PM, Sparks, John James wrote:
Dear R Helpers,
I am getting an error message from the try function that I don't
understand so I am hoping that someone can help.
I am scraping from web pages, but sometimes they disappear. When that
happens I need to control for it
Paul Johnson wrote
power - gamma ; x - Waist Size (cm) ; y - Weight (kg) # invalid R
names
plot(0, main=bquote(.(as.name(x))^.(as.name(power))/.(as.name(y
It appears to me that 2/3 of the as.name usages are not needed, at
least not in R 2.15 on Debian Linux.
This works just as well
On 12-05-18 12:56 PM, jaimie villanueva wrote:
hi
someone can show me how can i get the source code of a function. Is a S4
class or Method. (I'm not an expert in R environment)
Exactly, Function ugarchsim from library (rugarch).
I need to know (in detailed ) how the variance and mean ecuation
Many thanks to Baptiste Auguie (off list), Bill Dunlap, Bert Gunter,
Gabor Grothendieck,
and Paul Johnson for their interesting and informative responses to my post.
The main item to note is that I was piling on unnecessary baggage by using
plot(1:10,xlab=
On 12-05-18 2:23 PM, David Winsemius wrote:
On May 18, 2012, at 10:52 AM, Duncan Murdoch wrote:
On 12-05-18 5:15 AM, gholi bahrami wrote:
Dear Sir/Madam
I am a BSc student studying statistics in Ferdowsi University of
Mashhad,Iran.
Translating R messages to Farsi is my thesis and I have
Hi,
Is it possible to use x and y below to produce the vector shown in
desired.result (which represents the values in y ordered according to
the order specified by x)? Ideally the solution must be efficient in
terms of timing as I'm dealing with *very* long vectors.
x - c(7, 6, 5, 9, 4, 14, 8,
Hi all,
I have been struggling with ANOVAs on R. I am new to R, so I created a simple
data frame, and I do some analyses on R just to learn R and then check them on
SPSS to make sure that I am doing fine. Here is the problem that I've run into:
when we use the aov function, it uses SS Type I as
Hey,
Thanks for the reply. I tried tod do change this in two ways:
-with the formula, as you recommended (and after reading ?strptime some
similar formulas as well), but this did not help.
-by replacing the date to the form dd/mm/yy. But this didn't change anything
to the plot layout as well.
The
Hello,
Could you post the output of
dput(head(settlepw, 10))
dput(head(hatch, 10))
Just copy the output and paste it in a post. These functions give the first
10 lines (head) and the structure of your datasets.
I'm asking this because with the two example datasets you gave in an earlier
post,
Will
y[order(match(y,x))]
[1] 9 8 11 2 1
do?
Regards
Søren
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Lars Bishop
Sent: 20. maj 2012 01:05
To: r-help@r-project.org
Subject: [R] Help with ordering values
Hi,
Is it
On May 19, 2012, at 6:58 PM, jacaranda tree wrote:
Hi all,
I have been struggling with ANOVAs on R. I am new to R, so I created
a simple data frame, and I do some analyses on R just to learn R and
then check them on SPSS to make sure that I am doing fine. Here is
the problem that I've
Hello, All:
The writeFindFn2xls function in the sos package tries to
write an Excel file with 3 sheets ('PackageSum2', 'findFn', 'call').
Unfortunately, it is often unable to do this because of configuration
problems that are not easy to fix. I've found 3 contributed packages
that
Thank you Simon,
Yes this was in the library(car). There were no negative values in the
scatterplot graph or negative limits displayed.
Q-2 Are you sure the lines from the scatterplot are not this?:
solid green line: linear regression
solid red line: non-linear regression
upper dashed red line:
I have been using XLConnect to write multisheet Excel without any problems.
Sent from my iPad
On May 19, 2012, at 21:32, Spencer Graves
spencer.gra...@structuremonitoring.com wrote:
Hello, All:
The writeFindFn2xls function in the sos package tries to write an
Excel file with 3
74 matches
Mail list logo