Ruofei,
Ben's suggestion is simple and gets you close:
require(MASS)
nsim <- 100
rho <- -.9
Z <- mvrnorm(nsim, mu=c(0,0),Sigma = cbind(c(1,rho),c(rho, 1)))
U <- pnorm(Z);
a <- Z[,1]
b <- qunif(U[,2])
cor(a,b)
Pearson correlation characterizes the linear relationship between normal
r.v.'s, bu
alue, not a data frame. However, I am not
> sure. Any advice?
>
>
> Thanks for your help.
>
>
> Mike
>
>
> --
> *From:* Dalthorp, Daniel
> *Sent:* Wednesday, March 2, 2016 3:50 PM
> *To:* Michael
> *Cc:* r-help@r-project.org
> *Sub
Hi Michael,
If you are working in Windows:
# You can put the matrix directly into the clipboard
write.table(PRdist, file = 'clipboard', sep = '\t', row.names = F,
col.names = F)
The "sep" argument tells what character to use for separating columns.
Default for Excel is tab (i.e. '\t')
Default fo
Or, if you want easy labels, you can play around with contour graphs.
?contour # will give you info on how to make contour plots
The basic idea is to construct a matrix of z-values...one z for every
combination of x and y
contour(x,y,z)
The x's would then be the x-values you want in
(0.37273*log
A simple solution that will give you an idea of some of the plot parameters:
x<-seq(1,10,length=1000) # values for x-axis
x0<-c(0.4,0.5,0.6,0.7)
miny<-(log(min(x0))-(0.37273*log(max(x))-1.79389))/0.17941 # minimum
y-value to show on graph
maxy<-(log(max(x0))-(0.37273*log(min(x))-1.79389))/0.17941
Hi Val,
There are probably more elegant ways to do it, but the following is fairly
transparent:
# input data arranged as an array:
indat<-cbind(c(1,2,2,1),c(1,2,1,1),c(2,2,2,2),c(2,2,2,2),c(2,2,2,1),c(2,2,2,2),c(2,2,2,1),c(2,2,2,2),c(1,2,1,1),c(1,2,1,2))
indat
outdat<-array(dim=c(dim(indat)[1],di
Hi,
I'm currently working on an econometrics project on banking and looking to
merge a dataframe of bank specific data with dataframes of macro variables.
I am then going to transform the data set into a plm dataframe using the plm
package. The bank specific observations are indexed across time wh
Duncan Murdoch writes:
> On 29/01/2016 10:35 AM, Daniel Bastos wrote:
>> Here's how I plot a graph.
>>
>>plot(c(1,2,3), main = "graph ç")
>>
>> The main-string has a UTF-8 character "ç". I believe I'm using the
>> windows
Here's how I plot a graph.
plot(c(1,2,3), main = "graph ç")
The main-string has a UTF-8 character "ç". I believe I'm using the
windows device. It opens up on my screen. (The window says ``R
Graphics: Device 2 (ACTIVE)''.) How can I tell it to use my encoding of
choice?
I looked around the
Gwennaël,
Does the %% operator work for you?
It gives x mod y (or the remainder after dividing x into y...result is
guaranteed to be <=0 and >y)
E.g.
-150 %% 360 # 210
570 %% 360# 210
https://stat.ethz.ch/R-manual/R-devel/library/base/html/Arithmetic.html
-Dan
On Thu, Jan 28, 2
I'm finding it very difficult to figure out how to read the value of
"celltag" for a given cell in a tktable.
I'm sure it's something like:
tcl(classTable, "get", "celltag", row, column)
but of the dozens of variations of names, options, args, and formats I've
tried, nothing is working. Any sugg
I'd like to allow users to edit data in tcltk tables and to use vcmd to
validate data entry, e.g., not allowing non-numbers to be entered in
numeric cells and not allowing '\n' to be entered in text cells.
The problem is that I can't figure out how to "see" their data entry before
it is entered, a
a good one.
>
> Greetings,
>
> Adrian
> On Jan 21, 2016 10:36 PM, "Dalthorp, Daniel" wrote:
>
>> > Once you're up to speed on those issues...
>>
>> Any suggestions for getting up to speed on those issues?
>>
>>
>>
>&g
ethz.ch/pipermail/r-help/2009-January/378558.html
>
> -pd
>
> > On 21 Jan 2016, at 02:12 , Dalthorp, Daniel wrote:
> >
> > I know it should not be difficult to write the string:
> >
> > i<-4
> > j<-17
> > lbl<-paste0("[", i, &quo
> Once you're up to speed on those issues...
Any suggestions for getting up to speed on those issues?
On Thu, Jan 21, 2016 at 11:46 AM, peter dalgaard wrote:
>
> > On 21 Jan 2016, at 00:25 , Dalthorp, Daniel wrote:
> >
> > Thanks, Peter.
> >
> &g
I know it should not be difficult to write the string:
i<-4
j<-17
lbl<-paste0("[", i, ", ", j, "]")
# to a table, but I'm having a devil of a time trying to figure out how to
do it.
# the following gives lbl surrounded by braces.
tt<-tktoplevel()
tfr <- tkframe(tt)
tkgrid(tfr)
junk<-tclArray()
j
help.
E.g. with a table (called table1) with 3 columns and want to set widths to
30, 5, and 5:
colwidths<-c(30, 5, 5)
for(i in 1:3) {
tcl(table1, "width", i - 1, colwidths[i])
}
On Wed, Jan 20, 2016 at 3:07 PM, peter dalgaard wrote:
>
> > On 19 Jan 2016, at 20:48 , Da
Does anyone know a simple way to create a tcltk table with columns of
varying widths?
-Dan
--
Dan Dalthorp, PhD
USGS Forest and Rangeland Ecosystem Science Center
Forest Sciences Lab, Rm 189
3200 SW Jefferson Way
Corvallis, OR 97331
ph: 541-750-0953
ddalth...@usgs.gov
[[alternative HT
Very general question...try searching "R SQL" on google for a start.
On Fri, Jan 15, 2016 at 9:59 AM, Amoy Yang via R-help
wrote:
> Hi All,
> I am new here and a beginner for R. Can I use SQL procedure in R
> environment as it can be done in SAS starting with PROC SQL;
> Thanks for helps!
>
> A
What has happened?!
I get the following error message when I try to load rjags package (w/
Windows 7)...
> require(rjags)
Loading required package: rjags
Error in get(method, envir = home) :
lazy-load database 'C:/Program Files/R/R-3.2.3/library/rjags/R/rjags.rdb'
is corrupt
In addition: Warni
t;-matrix(runif(25),nrow=5,ncol=5)
# populating your 3-d array:
arrD3[,,1]<-mat1
arrD3[,,2]<-mat2
# to see what it looks like:
arrD3
On Fri, Jan 15, 2016 at 2:57 PM, Matteo Richiardi <
matteo.richia...@gmail.com> wrote:
> Hi Daniel,
> thanks for your answer. How can I
How about: D<-array(dim=c(d1, d2, d3))?
On Fri, Jan 15, 2016 at 2:20 PM, Matteo Richiardi <
matteo.richia...@gmail.com> wrote:
> What is the best way to store data in a cube? That is, I need to create a
> data structure D with three indexes, say i,j,h, so that I can access each
> data point D[i
How about this:
plot(0,0,xlab='',ylab='')
mtext(side=1,line=3,text='x axis',col=4)
mtext(side=2,line=3,text='y axis',col=2)
-Dan
On Mon, Jan 11, 2016 at 7:44 AM, Evan Cooch wrote:
> Consider a simple plot of X vs Y. There are elements on the plot that
> represent X, or Y, that are presented i
dev.off(which) can be used to close a specific graphics device where
"which" is the index of the device, but is there a way to assign a custom
number (or name) to a windows device so that specific window can be later
closed via dev.off (or some other method) if it is open?
The following does NOT w
eate a pdf document and the warning is cluttering
up my pdf document. Is there anyway to alter the plot statement to avoid
the warning? If not, is there a way to get knitr not to print out the
warning to the markdown document so that it doesn't end up in the pdf
document?
Thanks,
Thank you. This was helpful. I'll try it.
Daniel Wiegert
> On Nov 12, 2015, at 1:11 AM, Jim Lemon wrote:
>
> Hi Daniel,
> As David said, these problems are almost always due to fancy quotes. If you
> want plain old ASCII 34 (double) or 39 (single) quotes, try using
times for R to
accept them, or type them in manually in R.
Is there anyway to make this easier? What font in Word/Excel etc. does R prefer
and will function using quotes? Commas, letters, numbers, all work, just the
full double quotes won't work.
Thank you,
D
address I subscribed to the list with. When I resent it, I failed to
copy the subject line. My apologies.
- Daniel
On 10/30/2015 03:12 PM, Kevin Wright wrote:
Maybe you want
summary(aov(Y ~ A + Error(A:B)))
Kevin
On Fri, Oct 30, 2015 at 9:32 AM, Wagenaar, Daniel (wagenadl)
mailto:wagen...@ucm
`Sum Sq`[1] / a$`Sum Sq`[2]
P <- 1 - pf(F, a$Df[1], a$Df[2])
(They are not R programmers (yet).) And to be honest, I would find it
easier to read those results directly from the table as well.
Thanks,
Daniel Wagenaar
--
Daniel A. Wagenaar, PhD
Assistant Professor
Department of Biologica
On 10/30/2015 11:17 AM, Mark Leeds wrote:
Daniel: Just to complete my solution, here's the code for doing the
mean. Didn't expect this to take 3 emails !!! Have a good weekend.
temp <- tapply(f$value, f$justtimes, mean)
finalDF <- data.frame(chrontimes = times(rownames(temp)
uot;2015-10-30 00:50:00",
"2015-10-30 09:30:00", "2015-10-30 21:10:00", "2015-10-31 00:50:00",
"2015-10-31 10:30:00"), class = "factor"), value = c(88L, 17L,
80L, 28L, 23L, 39L, 82L, 79L)), .Names = c("date", "value"), ro
ow to get the "textbook" result in the most straightforward way?
(I'd like to be able to give me students a simple procedure...)
Thanks,
Daniel Wagenaar
--
Daniel A. Wagenaar, PhD
Assistant Professor
Department of Biological Sciences
McMicken College of Arts and Sciences
Univers
o convert strings to factors,
so you can modify as needed. In addition, if your files aren't as
regular as I inferred, you can increase the number of rows to read in
the first line to ensure getting the classes right.
Hope this is helpful,
Dan
--
Daniel Nordlund
Bothell, WA USA
___
x27;, believeNRows=FALSE, colQuote=NULL)
Hope this is helpful,
Dan
--
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
I am wondering why for generalized linear models with Gamma, Poisson and
Negative Binomial distributions that there appears to be no discussion
about estimating the medians or the modes of the distributions. For example
in clinical trials for count data where a log link is used it is the
quantity
times does not even
present itself here?
Thanks
Dan
On Fri, Jul 31, 2015 at 5:06 PM, Greg Snow <538...@gmail.com> wrote:
> Daniel,
>
> Basically just responding to your last paragraph (the others are
> interesting, but I think that you are learning as much as anyone and I
>
on). So I don't think so much that you are doing
> things wrong, just that you are learning that the models are complex.
>
> Another approach to simulation that you could try is to simulate the
> event time and censoring time using copulas (and therefore they can be
> correlated to give
dependent data and then transform the
> values using your Weibul. Or you could generate your event times and
> censoring times based on x1 and x2, but then only include x1 in the
> model.
>
> On Wed, Jul 22, 2015 at 2:20 AM, Daniel Meddings
> wrote:
> > I wish to simulate event
ecklenburg-Vorpommern",
"Niedersachsen",
"Nordrhein-Westfalen",
"Rheinland-Pfalz",
"Saarland",
"Sachsen",
"Sachsen-Anhalt",
"Schleswig-Holstein",
"Thüringen"), value = 1:16)
G3 <- gvisGeoMap(GR, locationvar
I wish to simulate event times where the censoring is informative, and to
compare parameter estimator quality from a Cox PH model with estimates
obtained from event times generated with non-informative censoring. However
I am struggling to do this, and I conclude rather than a technical flaw in
my
equal to
zero!
On Jul 13, 2015 5:31 PM, "Daniel Nordlund" mailto:djnordl...@frontier.com>> wrote:
On 7/13/2015 3:01 PM, Lida Zeighami wrote:
Hi there,
I have a matrix which its elements are 0, 1,2,NA
I want to remove the columns which the colsums are equ
ces
bootstrap coefficients so clearly I don't want to print the bootstrap
output but I do want users to be able to access it.
Many thanks,
Daniel
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE a
orld
problem. There is no error checking, and for large samples it may not
scale well.
Hope this is helpful,
Dan
--
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/l
m contain identical letters 1 and 3, AND ~50% again contain identical
letters 2 and 3 (except in this example as it is not possible from the choices).
Can multiple probability weightings be combined in such a manner?
Ben,
If I correctly understand your requirements, you can't do what you
lowing:
gaussianKernel <- function(u) exp(-u^2/2)/(2*pi)^.5
densityFunction <- function(x, df, ker, h){
difference = t(t(df) - x)/h
W = sum(apply(difference, 1, ker)) / (nrow(df)*h)
}
If you are wanting to do density estimation for real world work, I would
get help from someone i
;- rnorm(100,200,50)
# create deciles for binning
decile <- as.numeric(cut(rt, quantile(rt,0:10)/10),include.lowest=TRUE))
# collect into a dataframe (not really necessary)
df <- data.frame(rt=rt, decile=decile)
#compute the bin means
aggregate(rt,list(decile),mean,data=df)
This sho
On 5/16/2015 1:19 PM, Daniel Nordlund wrote:
On 5/16/2015 12:32 PM, li li wrote:
Hi all,
I wrote the following code and have two questions:
(1) As you can see, I would like different colors for different
types. It does not come out that way in the graph from this code.
Anyone know how to
over the colors chosen. I created a vector, color, specifying
the colors and the order they will be applied.
2. the scales parameter is what you are looking for. The parameter is
looking for a list which consists of name=value pairs. You can rotate
the labels by a specified number of deg
index=3)
Error in boot.ci(results, type = "bca", index = 3) :
object 'results' not found
A reproducible example means that when I run your code on my machine, I
get the same results / warnings / errors that you get. I got something
different.
Dan
--
Daniel Nordl
Has anybody noticed that if field URL in DESCRIPTION contains a uri with
66 or more characters, then file DESCRIPTION in the resulting package
includes a line break at the beginning?
So this (source DESCRIPTION):
URL:
http://ecdc.europa.eu/en/data-tools/seroincidence-calculator-tool/Pages/def
e help says of the formula
argument
formula
formula for restricted MIDAS regression or midas_r object. Formula must
include fmls function
your formula does not include the fmls() function, it uses mls(). So I
think your problem may have to do with how you are calling the midas
assuming you want to sample without replacement.
Generalizing it to other data structures is left as an exercise for the
reader.
replicate(100,mean(sample(yourdata,30, replace=FALSE)))
hope this is helpful,
Dan
--
Daniel Nordlund
Bothell, WA USA
__
R
-- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Daniel Nordlund
Bothell, WA USA
_
AS you could use PROC FCMP to
turn this into a function.
If you are interested contact me offline and I will send you a PROC FCMP
implementation.
Dan
--
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and mor
specific correlation between the series.
I have found some help here:
http://stats.stackexchange.com/questions/71211/how-to-generate-normal-random-variable-vector-which-is-spatially-auto-correlated
but I have a hard time implementing it in R, being a R-novice.
Thanks,
Daniel
[[alternative HTML
;FirstKey"))
rmatch.pos(app, "attr2", rmatch.pos(app, "SecondKey"))
rmatch.pos(app, "attr1", rmatch.pos(app, "ERROR"))
rmatch.pos(app, "ERROR", rmatch.pos(app, "attr1"))
On Tue, Dec 30, 2014 at 1:08 PM, Daniel Gabrieli
wrote:
>
uncan Murdoch
wrote:
> On 29/12/2014 4:41 PM, Daniel Gabrieli wrote:
> > I cross posted this on Stack Overflow:
> >
> http://stackoverflow.com/questions/27694466/chaining-multiple-replacement-functions-in-r
> >
> >
> > I am using R to work with a large JS obje
I cross posted this on Stack Overflow:
http://stackoverflow.com/questions/27694466/chaining-multiple-replacement-functions-in-r
I am using R to work with a large JS object (using the library rjsonio). As
such, I have a lot of nested lists, which are getting somewhat cumbersome
to work with. I hav
Subject: PLEASE HELP ON R-LANGUAGE
I am writing my project on "The effect of Bootstrapping on Time Series
Data". My problem is how I will write r programme to run a non stationary
time series data. I am familia to r language already.
Please help me.
[[alternative HTML version deleted]]
Good Morning Sir/Ma,
I am POPOOLA DANIEL a Forest Biometrician in making from the University of
Ibadan, Ibadan Nigeria. Please Sir/Ma I am having issues on performing
Non-linear mixed model on R (using maximum likelihood approach). I am trying to
input four different measured variables which
-04 2011-11-04
Pradip K. Muhuri, PhD
SAMHSA/CBHSQ
1 Choke Cherry Road, Room 2-1071
Rockville, MD 20857
Tel: 240-276-1070
Fax: 240-276-1260
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Daniel Nordlund
Sent: Sunday, November 09, 2
rg/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
I am not familiar with the mutate() function from dplyr, but you can get
your wanted results as follows:
data2 <- within(data1, oidflag <- apply(data1[,-1], 1, max, na.rm=TRUE))
Hope this is h
Hello, all!
So, as stated in the title, the Lme4 package used to output p-values for the
fixed effects. What happened?!
Literally 2 weeks ago, I ran code, got output with no errors, and had
p-values listed for my fixed effects.
Now, running THE SAME CODE with THE SAME DATASET (nothing at all ha
,list(Population), mean))
with(your_data_frame,aggregate(R,list(Population), var))
hope this is helpful,
Dan
Daniel Nordlund
Bothell, WA USA
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Maybe, you can add a parameter with a list of the column names and inside
the function an if statement, if the list lenght is 1 use the function with
1 column else the 3 columns. I am sure you can find better solutions.
Daniel Merino
2014-11-04 9:30 GMT-03:00 Frederic Ntirenganya :
> Hi Dan
Frederic,
Check the lubridate library.
install.packages("lubridate")
library(lubridate)
yday(Sys.Date())
d <- 4
m <- 11
y <- 2014
yday(as.Date(paste(y,m,d,sep="-")))
Daniel Merino
2014-11-04 7:01 GMT-03:00 Frederic Ntirenganya :
> Dear All,
>
> I would l
0+100, order.by=seq(as.Date("2014-10-22"),
as.Date("2014-10-31"), by="day"))
str(df)
Return.calculate(df, "log")
If you are going to use the PerformanceAnalytics package I highly recommend
you to checkl the xts package.
Daniel Merino
2014-10-31 15:24 GMT-
Bart,
Check if the following could help you.
library(xts)
y <- c(rnorm(10,25), rnorm(10,32),rnorm(10,25), rnorm(10,20),
rnorm(10,25)); x <- seq(c(ISOdate(2000,3,20)), by = "hour", length.out =
length(y))
z <- xts( y, order.by=as.POSIXct(x))
limit <- ifelse( lag(z) < 22
You are separating the objects using commas, the correct should be semicolons.
Example:
xi1=matrix(data=rep(0.0,600),ncol=3);xi2=matrix(data=rep(0.0,600),ncol=3))
Best regards,
Daniel Miquelluti
Em S�bado, 11 de Outubro de 2014 18:06, Uwe Ligges
escreveu:
On 11.10.2014 06:09, thanoon
re the rest un itl 1?
qreg1<-crq(Surv(TIME,EVENT,type="right")~VAR1+VAR2,
data=DATA_TRAIN,method = "Portnoy",grid=seq(0.1,1,by=0.1))
gives only one value for tau=0.0
I want to control the parameters tau or grid to get the results for
the quantiles I want, but it
] "="
# is trying to give me a pretty blatant hint... {: -)
On Mon, Sep 22, 2014 at 1:24 PM, Duncan Murdoch
wrote:
> On 22/09/2014 11:34 AM, Daniel Fuka wrote:
>>
>> Howdy Duncan,
>>
>> Thanks for the quick reply! I must be missing something
>> simple/obvi
release of the package.
>
> Bill Dunlap
> TIBCO Software
> wdunlap tibco.com
>
>
> On Mon, Sep 22, 2014 at 8:34 AM, Daniel Fuka wrote:
>> Howdy Duncan,
>>
>> Thanks for the quick reply! I must be missing something
>> simple/obvious. I need to have the &
my name is muka,\", x)"
return(song)
}
Thanks again for the quick reply and help you are giving me!
dan
On Mon, Sep 22, 2014 at 10:37 AM, Duncan Murdoch
wrote:
> On 22/09/2014 9:16 AM, Daniel Fuka wrote:
>>
>> Howdy,
>>
>> I have searched the list
Howdy,
I have searched the lists and can not seem to find a solution to my
problem. I need to be able to dynamically modify a string inside a
function to build a new function. "sub" replaces with a quoted
string... and "parse" of "sub" returns expression... How can I get an
unquoted string from a
I'm trying to set up an AR(2) model in the dlm context. I've generated a time
series utilizing the code:
am = 800; #sample size
des = 200; #initial values to be discarded
V = 0.5
v = rnorm((am+des+1),0,sqrt(V))
W = 0.9
w = rnorm((am+des+1),0,sqrt(W))
U = 0.9
R CMD check does not object to this code when checking a package:
foo1 <- function (bar) {
with(bar, {
x })
}
but produces a warning:
foo2: no visible binding for global variable 'x'
in response to this:
foo2 <- function (bar) {
within(bar, {
x })
}
Is this an R bug, o
uot;, but I need wt to be an argument in quotes, not an
object. It seems the issue is related to ?lm = "All of weights, subset
and offset are evaluated in the same way as variables in formula, that
is first in data and then in the environment of formula., but I can't
figure it out. Can you
results
closely match those from Epicure.
Using the data here: http://dwoll.de/err/dat.txt
The stan model fit below replicates the results from Epicure here:
http://dwoll.de/err/epicure.log
Of course I am still interested in learning about other options or approaches.
Daniel
##---
## rstan
ndows 8.1
Do I need any other installations o packages or something else?
The code was working in older Windows XP and Linux
Thanks
--
**
Juan Daniel GarcÃa Villabrille, Ingeniero de Montes, colegiado nº
, p=probs, replace=TRUE)
}
new_sample(x)
Daniel Nordlund
Bothell, WA USA
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Dan Abner
> Sent: Monday, June 23, 2014 3:19 PM
> To: Greg Snow
> Cc: r-help@r-p
I have coded qualitative data with many (20+) different codes from a survey
in an excel file. I am using the with() function to output the codes so we
know what's there. Is it possible to direct the output from with() to an
excel file? If not, what's another function that has the same, er,
function
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Daniel Nordlund
> Sent: Thursday, May 01, 2014 1:10 AM
> To: r-help@r-project.org
> Subject: Re: [R] A combinatorial assignment problem
>
>
r,m,max_iter=120) {
n <- 0
cmb <- combn(r,m)
repeat {
n <- n+1
tbl <- table(map<-cmb[,sample(1:choose(r,m),k)])
if(min(tbl) == max(tbl)-1) break
if(n > max_iter) break
}
return(t(map))
}
a <- assignment(10,7,3)
Dan
Daniel Nordlund
freezes R on my Win 7 Pro x64 box using either 64-bit R-3.0.3 or R-3.1.0. You
might try switching to shell() instead of system()
> command <- paste(aa, fnm)
> shell(command)
However, it all depends on what programs you are trying to run and what
behavior you expect.
Dan
Daniel
sp3 = structure(c(1L, 1L, 1L, 1L), .Label = "mean", class =
> "factor")), .Names = c("watershed",
> "year", "sp1", "sp2", "sp3"), class = "data.frame", row.names = c(NA,
> -4L))
>
> Any s
Thanks. This helped greatly. I’m sorry about sending the html message, and
hope this one is plain-text. Dan.
On Feb 21, 2014, at 7:00 AM, Prof Brian Ripley wrote:
> On 21/02/2014 10:40, Daniel Kelley wrote:
>> I’m wondering whether anyone can help me with a translation exercise. I
Im wondering whether anyone can help me with a translation exercise. I have a
package named oce, which does oceanographic processing, and Id like to make
it produce graphs with labels that work in different languages. For example,
in English I write Depth and in Spanish Id like to write
as cut-n-paste from your email) and got a plot
of the data with regression line without any error message. What version of R
and ISwR are you running? Have you tried closing R and starting a clean
session?
Since the example worked for me, there is not much else I can help with.
Dan
Dan
e of 1.5 hours:
>
> > as.POSIXct("2010-04-04 03:00:00")- as.POSIXct("2010-04-04 02:30:00")
> Time difference of 1.5 hours
>
>
>
>
>
>
> Any suggestions?
> Thanks,
> David Fox.
>
Daylight savings time change in Australia?
Dan
Da
rrect and fluid communication"
--
****
Daniel Patón Domínguez
Numerical Ecology. Ecology Unit
Department of Plant Biology, Ecology and Earth Sciences
Faculty of Sciences. University of Extremadura
Avda. Elvas s/n 06071 Badajoz (Spain)
http://unex.academ
Dear all:
I want to predict a presence/absence vector using a presence/absence matrix of
events. What library can do this in R?
Many thanks
--
Daniel Patón Domínguez
Numerical Ecology. Ecology Unit
Department of Plant Biology
read up on
?sample
Dan
Daniel Nordlund
Bothell, WA USA
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Alaios
> Sent: Monday, January 20, 2014 10:12 PM
> To: Alaios; R-help@r-project.org
> Subject
Thanks for the reply. Another great option would be "missing" (like in SAS),
especially for factors. I'm struggling to figure out how to do this with
"tables".
Daniel Cher, MD
djc...@gmail.com
+1-650-269-5763
This message and its attachments are confidential.
-O
All 8 100
>
I'm looking for percents to have gender=1 or gender=2 as the denominator.
I.e.,
Gender Race n Percent
1 32 *50*
42 *50*
2 30 *0*
44 *100*
All 8 100
>
event and offset pyears.
Many thanks, D
> -Original Message-
> From: David Winsemius [mailto:dwinsem...@comcast.net]
> Sent: Thursday, January 09, 2014 4:33 AM
> To: Wollschlaeger, Daniel
> Cc: r-help@r-project.org
> Subject: Re: AW: [R] Linear relative rate / excess re
fset.
Unfortunately, I won't have comparable baseline rate tables. And while I could
fit a separate model only to the unexposed group for expected counts, I'd
prefer to fit both factors (lambda0 and 1+ERR) simultaneously - as it is
typically done in the existing literature.
o me (eg., gnm,
timereg) but I currently don't see how to correctly specify the model.
Any help on how to approach ERR models in R is highly appreciated!
With many thanks and best regards
Daniel
[1] Preston DL. Beyond Dose Response: Describing Long-Term Health Effects of
Radiation Expos
a try, but could not find the problem there either.
The code below was pasted though textEdit and converted to plain text. I
hope this takes care of any embedded characters.
\documentclass{article}
\begin{document}
<>=
library(knitr)
library(ggplot2)
@
\title{Knitr and ggplot2}
\author{
I am forwarding to r-help just in case anyone else where wondering or have
the same problem running the example.
Cheers
D
-- Forwarded message --
From: Daniel Haugstvedt
Date: Sun, Dec 22, 2013 at 11:54 PM
Subject: Re: [R] Knitr, ggplot and consistent fonts
To: John Kane
Hi
work.) There has to be others facing the same problem and someone
must have found a nice solution.
Additional attempts from my side which failed are not included in the
example. I have tested the Google results i could find without any luck.
Cheers
Daniel
PS. I know the example plots could have
101 - 200 of 1303 matches
Mail list logo
