kchkchkch wrote
>
> Heh. Yes,
> Berend Hasselman wrote
>>
>> Shouldn't the line with uniroot be (why try as function?)
>>
>> sol = uniroot(fcn, lower = p1, upper = 1)
>>
>> Before the line with for(i in 1:n) insert the following
>>
>> fcn = function(p2) p1*f(p1) + (.20/5.66)*(exp(5.66*(p2
Dear R helpers,
I have a question about drawing random numbers from many categorical
distributions.
Consider n individuals, each follows a categorical distribution defined
over k categories.
Consider a simple case in which n=4, k=3 as below
catDisMat <-
rbind(c(0.1,0.2,0.7),c(0.2,0.2,0.6),c(0.1,
Hi,
I use the function fpot of packages evd. If I call the fit that I obtain
"fit", I want to modify the name of the axis and the main title that is
produced by plot(fit1). To do this, I want to create a new function where I
would have the names modified. The problem is that I can't find the funct
I am still interested in this. Is there no way to pay directly online? via
paypal or other?
Thanks,
Xu
--
View this message in context:
http://r.789695.n4.nabble.com/Re-Monetary-support-to-the-R-project-Was-Re-Executable-for-Production-Use-tp1585186p4198369.html
Sent from the R help mailing lis
Dear R-users,
I am a relatively new R user and I have a simple question. Is there a
command attached to a R function, which performs matrix pooling?
Thanks in advance,
Jana
--
--
"The world is full of mysteries. Life is one. The curious limitations of
finite minds are another."(J.B.S. Halda
Hello,
I have a data frame with hourly or sub-hourly weather records that span
several years, and from that data frame I'm trying to select only the
records taken closest to noon for each day. Here's what I've done so far:
#Add a column to the data frame showing the difference between noon and th
I tried using unzip but it gave the following error "Invalid zip name
argument"
options(HTTPUserAgent="Firefox/8.0.1")
u =
"http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip";
o = getURLContent(u, verbose = TRUE, useragent = getOption("HTTPUserAgent"))
afte
On 2011-12-14 18:32, Carl Witthoft wrote:
That only puts different colors on different axes. I was wondering
about the (silly) possibility of having a range of colors for the ticks
on a given axis.
I think you would have to either hack the code or issue
several axis() calls, e.g.
plot(1:9,
Bert,
The OP may be looking for recordPlot/replayPlot.
We were not told the OS (x11 exists even on Windows). The standard
screen devices on Windows (windows() ) and Mac OS X (quartz()) have plot
histories that can be used to flick between plots. And 'flick' is the
operative word here: it is
That only puts different colors on different axes. I was wondering
about the (silly) possibility of having a range of colors for the ticks
on a given axis.
From: Hans W Borchers
Date: Tue, 13 Dec 2011 14:21:31 +
Carl Witthoft witthoft.com> writes:
>
> Hi,
> So far as I can tell, the '
library(openair)
?windRose
--
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600FAX:(360) 407-7534
Olympia,
There is also a windrose (though not with that name) function in the
openair package.
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 12/12/11 1:20 PM, "Adrienne Wootten" wrote:
>Greetings!
>
>I'm having an issue with t
Please follow the posting guide and provide some example code.
Bottom line, most functions that "draw" a graph return a graphic
object that, like any other R object, can be given a name in the
global environment, put into a list, etc.
Alternatively, depending on what you mean -- it was not entire
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of William Dunlap
> Sent: Wednesday, December 14, 2011 4:19 PM
> To: Trying To learn again; r-help@r-project.org
> Subject: Re: [R] Omit Inf
>
> > range(a, finite=TRUE)[2] # [1] fo
> range(a, finite=TRUE)[2] # [1] for the min
[1] 2
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Trying To learn
> again
> Sent: Wednesday, December 14, 2011 3:55 PM
Hi all
> a<-c(2,Inf)
> max(a)
[1] Inf
>
How can I say...omit Inf?
I cannot find anything...
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guid
Dear R members,
I would like to store 18 graphics at the end of a loop (as we do for data
frames and arrays). For each iteration, I use the function x11() and I want
to keep my graphs in a single window to re-use them later. Do you know a
function for storing graphics ?
Many thanks
Julie
--
Vi
Hi everyone,
I have been following these forums closely for the last few months but this
is my first time posting. Basically I am trying to get an Oracle SQL
Developer Binary Large Object(BLOB/CLOB) file in as an R object in R. For
those not familiar with a BLOB/CLOB file, it is basically lik a ta
Thanks David and Patrick,
really appreciate tour help...
m.
Il giorno 14 Dec 2011, alle ore 18:33, David Winsemius ha scritto:
>
> On Dec 14, 2011, at 11:43 AM, matteo wrote:
>
>> set.seed(3)
>> mydata <- data.frame(var = rnorm(100,20,1),
>> temp = sin(sort(rep(c(1:10),10))
Hello,
I am writing a test function, and would like to have it return an
appropriate 'htest'
object. This is very easy, but the test can be run for several orders, a
frequent situation
in time series tests.
It would be nice to return a data.frame with order, params, test statistic,
p.value.
Afte
Hi all,
I have solved my question
Output<- matrix(0, length(x), 1)
for (i in 1:length(x)){
files <- paste('KT', i, '.csv', sep = '')
da1 <- read.csv(files)
Output[i,1]<-da1[1,2]
}
2011/12/14 Trying To learn again
> Hi all,
>
> I have 100 csv files always with this information (I Attac
First PRoblem solved
> read.csv("KT80.csv")
X x
1 1 0.01331361
>
I ommitt "" in calling the file name...
2011/12/14 Trying To learn again
> Hi all,
>
> I have 100 csv files always with this information (I Attach two example
> excels)
>
> KT80.csv contains:
> ,"x"
> 1,127.188
On Wed, Dec 14, 2011 at 4:43 PM, Greg Snow wrote:
> There is also the problem that SWords license does not allow for
> commercial use. R2wd, write.table, and odfWeave dont have this
> restriction.
>
The free license for SWord does have the commercial use restriction.
There is also a commeric
There is also the problem that SWord's license does not allow for commercial
use. R2wd, write.table, and odfWeave don't have this restriction.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
From: Richard M. Heiberger [mailto:r..
Richard,
I have looked at SWord before, but to my knowledge it does not deal directly
with the tabular objects created by the tables package (please correct me if I
am wrong). These objects do have a matrix as the main data, but the attributes
are different from the usual dimnames. There are
On 14.12.2011 19:43, Cable, Sam B Civ USAF AFMC AFRL/RVBXI wrote:
Am trying to install package JGR and, by necessity, JavaGD on a few
Linux platforms. R rev. 2.14.0. Had some success on one platform, but
am running into a problem on a new platform. I installed rJava, JavaGD,
iplots, and JGR.
Do unz / unzip get the job done?
Michael
On Wed, Dec 14, 2011 at 8:19 AM, himanshuy wrote:
> Hi
>
> I am trying to download and unzip a zip file from the following location.
>
> http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip
>
> I initially there was issu
On Wed, Dec 14, 2011 at 7:04 AM, agent dunham wrote:
> Thank you,
>
> I guess it didn't work for me, maybe is not possible?
>
> I've tried:
>
>> con<- odbcDriverConnect("Driver=SQL Server;
>> Server=...\\...;Database=...;Uid=...;Pwd=... ;")
>
>> v1=sqlQuery(con, "select v1 from sqltable where v3 =
Heh. Yes,
Berend Hasselman wrote
>
> Shouldn't the line with uniroot be (why try as function?)
>
> sol = uniroot(fcn, lower = p1, upper = 1)
>
> Before the line with for(i in 1:n) insert the following
>
> fcn = function(p2) p1*f(p1) + (.20/5.66)*(exp(5.66*(p2 - p0)) -
> exp(5.66*(p1 - p0)))
I have a big SQL Server application that needs some help from R, and this R
needs to load a large workspace in order to service the database request. I'd
like to keep an R process running continually because loading the workspace
takes about 1 minute. The data in the workspace doesn't change ver
On Wed, 14-Dec-2011 at 09:55AM -0800, Peter Langfelder wrote:
|> On Wed, Dec 14, 2011 at 4:35 AM, Chris Mcowen wrote:
|> > Dear List,
|> >
|> >
|> >
|> > I am unsure if this is the correct list to post to, if it isn't I
apologise.
|> >
|> >
|> >
|> > I am using SSH to access a Linux version of R
On Wed, 14-Dec-2011 at 08:43AM -0800, matteo wrote:
|> Dear all,
|>
|> here is a simple problem that surely you already come across, but is giving
|> me a big headache...
|>
|> I have a dataframe like this:
|>
|> set.seed(3)
|> mydata <- data.frame(var = rnorm(100,20,1),
|>
Am trying to install package JGR and, by necessity, JavaGD on a few
Linux platforms. R rev. 2.14.0. Had some success on one platform, but
am running into a problem on a new platform. I installed rJava, JavaGD,
iplots, and JGR. Things went - I thought - without a hitch. Now, when
I try to use J
On Dec 14, 2011, at 10:01 AM, Tony Stocker wrote:
All,
Given the following commands:
ag.m35<-read.table("m35.txt",header=TRUE,sep=",")
ag.m35.lp<-subset(ag.m35, race=="lp")
aov.m35.lp=aov(year~time,data=ag.m35.lp)
anova.m35.lp=anova(aov.m35.lp)
anova.m35.lp
Analysis of Variance Table
Respo
On Dec 14, 2011, at 11:43 AM, matteo wrote:
set.seed(3)
mydata <- data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
subj = as.factor(rep(c(1:10),5)))
and I need to make a scatter plot for each subj, not a problem, but...
what i want is
Greg,
Please look at the SWord package. This package integrates MS Word with R
in a manner similar to the SWeave integration of LaTeX with R.
Download SWord from rcom.univie.ac.at
If you have a recent download of RExcel from the RAndFriends installer, then
you will already have SWord on your mach
1) Use dput() to submit data.
2) Would this work? (It requires your data are evenly spaced, but I
think that's it) d[seq(1, nrow(d), by = 3), ]
Michael
On Wed, Dec 14, 2011 at 7:17 AM, abcdef ghijk wrote:
> Good Morning ,
>
> I want to sample the following time series for every third hour. For
kchkchkch wrote
>
> I have one equation, two unknowns, so I am trying to build the solution
> set by running through possible values for one unknown, and then using
> uniroot to solve for the accompanying second solution, then graphing the
> two vectors.
>
> p0 = .36
> f = function(x) 0.29 * exp
On Wed, Dec 14, 2011 at 4:35 AM, Chris Mcowen wrote:
> Dear List,
>
>
>
> I am unsure if this is the correct list to post to, if it isn't I apologise.
>
>
>
> I am using SSH to access a Linux version of R on a remote computer as it
> offers more memory and processing power. The model will take 1-2
What I did in the past (not with R scripts) is to start my jobs using
at (start the job at a specified time e.g. now) or batch (start the
job when the cpu drops below ?%)
at now "R CMD BATCH yourscript.R"
or
batch "R CMD BATCH yourscript.R"
something like that, you'll have to look at the m
On Wed, 14 Dec 2011, Chris Mcowen wrote:
As a Linux newbie, I was wondering if anybody here knew how to keep R
running and interactive and return to it on a later date?
Chris,
Screen will retain the session and keep the process going. Read 'man
screen'.
Whether this works from a remote M
Colors probably are not the best for so many levels and combinations. Look at
the symbols function (or the my.symbols and subplot functions in the
TeachingDemos package) for ways to add symbols to a map showing multiple
variables.
-Original Message-
From: r-help-boun...@r-project.org [
Duncan,
If you are taking suggestions for expanding the tables package (looks great)
then I would suggest some way to get the tables into MS products. If I create
a full output/report myself then I am happy to work in LaTeX, but much of what
I do is to produce tables and graphs to clients tha
All,
Given the following commands:
> ag.m35<-read.table("m35.txt",header=TRUE,sep=",")
> ag.m35.lp<-subset(ag.m35, race=="lp")
> aov.m35.lp=aov(year~time,data=ag.m35.lp)
> anova.m35.lp=anova(aov.m35.lp)
> anova.m35.lp
Analysis of Variance Table
Response: year
Df Sum Sq Mean Sq F value
Good Morning ,
I want to sample the following time series for every third hour. For example at
00:00,03:00,06:00,09:00 etc.
2011-01-01 00:00:00 0.00e+00
2011-01-01 01:00:00 1.471667e+01
2011-01-01 02:00:00 1.576667e+01
2011-01-01 03:00:00 0.00e+00
2011-01-01 04:00:
Dear all,
here is a simple problem that surely you already come across, but is giving
me a big headache...
I have a dataframe like this:
set.seed(3)
mydata <- data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
subj = as.factor(rep(c(1:1
I have one equation, two unknowns, so I am trying to build the solution set
by running through possible values for one unknown, and then using uniroot
to solve for the accompanying second solution, then graphing the two
vectors.
p0 = .36
f = function(x) 0.29 * exp(5.66*(x - p0))
f.integral = integ
Hi
I am trying to download and unzip a zip file from the following location.
http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip
I initially there was issue with useragent so i ran the following code
options(HTTPUserAgent="Firefox/8.0.1")
u =
"http://www.nse
Hi R users,
I'm using the 'bam' function in mgcv to examine trends in a remotely sensed
vegetation index. I have one random effect variable, 'cons' (with six
levels) which identifies different subjects within this analysis.
My model is specified as follows:
rm4<-bam(trend~factor(zone)+s(cons,bs
Sorry for some typo last mail.
I corrected it,and resent.Sorry for it.
Sir:
I find out that for 2 level factor, if I set it to "factor", then I'll get
error reply.
For the instance last mail, if I use:
result1<-glm(y ~ factor(gender),family = binomial);#gender has 2 levels
logistic.display
Thank you,
I guess it didn't work for me, maybe is not possible?
I've tried:
> con<- odbcDriverConnect("Driver=SQL Server;
> Server=...\\...;Database=...;Uid=...;Pwd=... ;")
> v1=sqlQuery(con, "select v1 from sqltable where v3 =cte and v2 in (select
> v2 from R_dataframe) order by (select v2
Dear List,
I am unsure if this is the correct list to post to, if it isn't I apologise.
I am using SSH to access a Linux version of R on a remote computer as it
offers more memory and processing power. The model will take 1-2 days to
run, I am accessing R through Putty and when I close the
Emma,
That is because generate_unit expects a data.frame with one row and
columns id and size:
generate_unit(data.frame(id=1, size=10))
Jan
Emma Thomas schreef:
Dear Jan,
Thanks for your reply.
The first solution works well for my needs for now, but I have a
question about the sec
Actually, scratch that, sorry!
I put the second part of your second solution code into a function and get the
right data frame in the end. So:
generate_unit<- function(unit) {
pid<- 1:unit$size
senior<- rep(0, unit$size)
senior[sample(unit$size, 2)] <- 1
return(data.frame(unit_id=unit$id, pid=p
Dear Jan,
Thanks for your reply.
The first solution works well for my needs for now, but I have a question about
the second. If I run your code and then call the function:
generate_unit(10)
I get an error that
Error in unit$size : $ operator is invalid for atomic vectors
Did you experience
Forgot to cc the list.
M
On Wed, Dec 14, 2011 at 12:19 PM, R. Michael Weylandt
wrote:
> It's a bit of a hack-job, but this might get you started:
>
> d <- structure(c("A", "A", "A", "B", "A", "A", "C", "C", "A", "D",
> "C", "C"), .Dim = 3:4)
>
> plot(0, type = "n", xlim = c(0,ncol(d)), ylim = c(
Thank you very much, Prof Ripley. The problem is solved and my understanding is
improved.
Happy holidays!
- Alex
From: Prof Brian Ripley
To: Alex Zhang
Cc: "r-help@R-project.org"
Sent: Wednesday, December 14, 2011 11:39 AM
Subject: Re: [R] A Question Re
On 14/12/2011 16:24, Alex Zhang wrote:
Dear all,
I am using the ISOdatetime function like this:
test.info$TradeTime = with(test.info, mapply(FUN = ISOdatetime, Year, Month, Day, Hour,
Minute, 60, "EST"))
Where Year, Month etc are all numeric.
I think ISOdatetime should return a POSIXct obje
I may be confused, but this seems ill-defined since, in the example,
the color corresponds to a numerical value: characters are only there
as labels. Do you intend to map the characters onto some numerical
spectrum? If so, how?
As a side note, one of the best places to get visualization ideas in
m
Dear Jean,
That is a simple and effective solution. Thank you!
Best,
Shige
On Tue, Dec 13, 2011 at 4:16 PM, Jean V Adams wrote:
>
> Shige Song wrote on 12/12/2011 07:49:39 PM:
>
>
>> Dear All,
>>
>> I am fitting a simple GAM model using the "gam" function in the "mgcv"
>> package. The only inde
Hi, I am looking for options to plot the following type of matrices:
"A" "B" "C" "D"
"A" "A" "C" "C"
"A" "A" "A" "C"
as a image like this: http://www.phaget4.org/R/image002.jpg
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
Dear all,
I am using the ISOdatetime function like this:
test.info$TradeTime = with(test.info, mapply(FUN = ISOdatetime, Year, Month,
Day, Hour, Minute, 60, "EST"))
Where Year, Month etc are all numeric.
I think ISOdatetime should return a POSIXct object. However, the result I
obtained from
> Have a left-censored dataset, attempting to use a Tobit model and am
working
> with the censReg package. I like how easy it is to move from glm
models to
> predictions with 'predict' and wanted to ask if there was a way to do
so
> with objects output from 'censReg'.
You can fit tobin regres
arunkumar gmail.com> writes:
> Please any one help in finding the predicted value for lmer function
>
> model<- lmer(formula =formula,data=data,verbose=TRUE,family = "gaussian")
>
> I need to get predicted value for this model.
Please see the code at http://glmm.wikidot.com/faq , and
di
Emma,
If, as you say, each unit is the same you can just repeat the units to
obtain the required number of units. For example,
unit_size <- 10
n_units <- 10
unit_id <- rep(1:n_units, each=unit_size)
pid <- rep(1:unit_size, n_units)
senior <- ifelse(pid <= 2, 1, 0)
pop <- d
Hi robin,
I'm not sure is what you need, but that's an esthetically nice
solution (one single line without any loop :) )
matrix(apply(log(cbind(as.numeric(a),as.numeric(b),as.numeric(c),as.numeric(d))),1,sd),3)
hope it could help,
PF
On Mon, Dec 12, 2011 at 5:15 PM, Robin Cura wrote:
> Hello,
Hi,
I am trying to use the system2 function to execute external
applications(actually it is psql) inside R 2.14.0 on X64 Windows, but
the psql command has escape character inside it's full command line,
can you help to figure out the correct parameters for the system2
function,
The working comman
Hi Sacha,
I was referring to the BRugs package, where I usually proceed as follows:
result = BRugsFit (...)
param_samples = samplesSample("param")(1)
But (1) takes a long time for big chains. So my question was specific to
BRugs, if there an alternative way to get the samples.
openbugs() m
68 matches
Mail list logo