Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal random
sample with missing values and ran into an unexpected result from my call to
dmvnorm() within the likelihood function. Particular details are provided
below. It appears that dmvnorm() makes a call to
I have a simple table, 2 columns and 1994 rows. First row, Crime is how
many crimes happen every month per 10 inhabitants and second row is
Income which contains the average income recorded in a city.
here's the head(dataset):
Crime Income
1 356.5152 4285.720
2 734.5625 4114.291
3
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal random
sample with missing values and ran into an unexpected result from my call to
dmvnorm()
There are at least 5 different version of dmvnorm. None of them
glsnow wrote
From your question it is not clear what your question/concerns really are,
and from what we can see it could very well be that you do not understand
the statistics that you are computing (not just the R implementation). We
ask for a reproducible example because that helps us to
On Oct 18, 2013, at 08:37 , David Winsemius wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal
random sample with missing values and ran into an unexpected result from my
call to dmvnorm()
There
Hi
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Manuel Figueroa
Sent: Friday, October 18, 2013 6:24 AM
To: r-help@r-project.org
Subject: [R] (no subject)
I have a simple table, 2 columns and 1994 rows. First row, Crime
Thank you Greg
I will likely go with a list of data frames as each data frame will
represent independent entities and as the syntax will probably be
better that way.
I did find it puzzling that multidimensional data frame or the ability
to hold a vector in a cell of a data frame was not added in
I don't get an error message. However, I had to make up fstatus and all the
rest because your example is not reproducible. Please provide the list with
more information and you'll likely get a more helpful answer from someone.
-Original Message-
From: Elan InP
Hi all
I do not understand why I am getting the following error message. I am showing
all my code this time. Can
anybody help me with this? Thanks in advance.
install.packages(cmprsk)
library(cmprsk)
ndim = 5
# two covariates standard normal variates
z1 = rnorm(ndim)
z2 = rnorm(ndim)
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Subject: [R] Lattice xyplot: Fill Legend Points
When I specify pch = 19 for a scatter plot the points are filled
circles.
Deapite reading ?points and
On Fri, 18 Oct 2013, S Ellison wrote:
Using trellis.par.set seems more successful than specifying par.settings
in the call, however much I'd prefer to avoid relying on globals. Example:
ss-trellis.par.get(superpose.symbol)
ss$pch=rep(19,7)
ss$col=rainbow(7)
trellis.par.set(superpose.symbol=ss)
Dear Duncan,
Thank you for your quick reply. I've got the basic version of what I'm
looking for now (see below). My next step will be your rgl::persp3d suggestion
for the hidden lines control.
Best,
Colin
library(MASS)
X - mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,1),2))
X.kde -
-Original Message-
... the kindest guide I can give is to
read an Introduction to R (ships with R) or a R web tutorial of your choice
No quibble with the advice, but it prompted me to look again at the R Intro.
Interestingly, the Intro doesn't mention subset() at all; the subsetting
-Original Message-
mydat
basel_asset_class defa_frequency
1 2 0.150
2 8 0.070
3 8 0.030
4 8 0.001
I need to get the subset of this data.frame where no of records for the
On 18-10-2013, at 14:06, Elan InP englishinpa...@aim.com wrote:
Hi all
I do not understand why I am getting the following error message. I am
showing all my code this time. Can
anybody help me with this? Thanks in advance.
install.packages(cmprsk)
library(cmprsk)
ndim = 5
#
Thanks for your help David!
I was running the same code the other day and it worked fine although it
took a while as well. You are right that dff shud be df1 and maybe it's a
portion of my data so it have an error of length =0.
About CPU usage, I got it by clicking ctrl+alt+delete and it showed
Hi,
I just got started with R. I am trying to load some data into PostgreSQL
using RPostgreSQL. Everything went quite smoothly except for that the table
created by using dbWriteTable does not have a primary key.
The ideal solution is to ask PostgreSQL to do a auto-increment with the
rows I am
Thanks. This was helpful.
Hermann
2013/10/17 Bretschneider (R) brets...@xs4all.nl
On 17 Oct 2013, at 13:44 , Hermann Norpois wrote:
Hello,
my dots of 0 and 2 are quite close to the marging. So I would like to move
the 0 and the 2 both towards the 1. I wish to be my dots more centered.
Dear sir,
Thanks a lot for your guidance. I have been benefited immensely by this
discussion. Thanks again.
Regards
Katherine
On Friday, 18 October 2013 2:50 AM, Bert Gunter gunter.ber...@gene.com wrote:
Thanks, Bill.
But ?ave specifically says:
ave(x, ..., FUN = mean)
Arguments:
x
A
Hi,
I'm trying to combine 15 different data frames but all have time field as
common. But the thing is time intervals are different (like min's and
sec's).
So I have to check the common time from all the 15 data frames and merge it
to a single data frame and then export it to excel or csv.
You might want to use the profiler (Rprof) on a subset of your code to
see where time is being spent. Find a subet that runs for a minute,
or so, and enable profiling for the test. Take a look and see which
functions are taking the time. This will be a start. You can also
watch the task monitor
You at least need to provide a subset of what the data looks likes.
Since you are merging on 'time', if there are not the same times in
each set, then you will have missing values. For example, in the
merges that you are doing, have you looked at the size of the
resulting object compared to the
Hi,
Assuming that you provided the sample data from the file.
temp - readLines(textConnection(#Hogd/met, Temp, 005[M], Value
#Hogd/met, Difftemp, 051[M], Value
BA0+
1 MTEMP005 1 [deg.C]
2 MDTMP051 1 [deg.C]
EOH
891231, 2400, -1.5, -0.21,
900101, 0100, -1.4, -0.25,
900101, 0200, -1.6, -0.28,
On Oct 17, 2013, at 11:37 PM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal
random sample with missing values and ran into an unexpected result from my
Dear all,
I have a list that is created like that
Spans-list( c(837e6,842e6),
c(832e6,837e6),
c(930.1e6,935.1e6)
)
I would like to include a second list that will contain the string that would
correspond to the numbers at the left side.
I would like thus
On Oct 18, 2013, at 1:12 AM, peter dalgaard pda...@gmail.com wrote:
On Oct 18, 2013, at 08:37 , David Winsemius wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal
random sample with missing
Hi all,
I am using the ... argument to parmeterize a user define fn to accept
multiple input objects. I subsquently save all these data as a list.
Question: what is the best way to recover or extract the original object
names that were fed to the fn?
Thanks,
Dan
[[alternative HTML
I'm not exactly sure what you mean by names. Does the following meet
your needs?
f - function(...)names(list(...))
f(a=2,b=3)
[1] a b
f(a=2,3)
[1] a
If not, a reproducible example of what you want might be helpful.
Cheers,
Bert
On Fri, Oct 18, 2013 at 9:05 AM, Dan Abner
On 10/18/2013 9:10 AM, Berend Hasselman wrote:
On 18-10-2013, at 14:06, Elan InP englishinpa...@aim.com wrote:
Hi all
I do not understand why I am getting the following error message. I am
showing all my code this time. Can
anybody help me with this? Thanks in advance.
Hi,
Not sure this is what you wanted.
lstNew - list(Spans, lapply(lapply(Spans,`[`,1),as.character) )
str(lstNew)
#List of 2
$ :List of 3
..$ : num [1:2] 8.37e+08 8.42e+08
..$ : num [1:2] 8.32e+08 8.37e+08
..$ : num [1:2] 9.30e+08 9.35e+08
$ :List of 3
..$ : chr 8.37e+08
..$ : chr
Hello all.
I have this code:
myLOOCV - function(myformula, data) {
Y - all.vars(myformula)[1]
Scores- numeric(length(data[,1]))
for (i in 1:length(data[,1])) {
train - data[-i,]
test - data[i,]
myModel - lrm(myformula, train)
Scores[i] -
Dear list,
I am using the hurdle model for modelling the habitat of rare fish species.
However I do get an error message when I try to model my data:
test_new1-hurdle(GALUMEL~ depth + sal + slope + vrm + lat:long +
offset(log(haul_numb)), dist = negbin, data = datafit_elasmo)
Error in
Hello all
I've been performing a series of multidimensional optimizations (3 variables)
using the optima() function. Recently, I noticed that the solution is rarely
unimodal. Is there a package or function that handles multimodal
multidimensional optimizations? I really appreciate any
Hi everyone,
If I compute a Ordered Logistic or Probit Regression with the polr
function from MASS package. the summary give me : coefficients, Standard
error and Tvalue.. but not directly the p.value.
I can compute manualy the Pvalue, but Is there a way to directly obtain
the pa.value, and I
Thanks for your advice Jim!
I tried Rprof but since the code just freezes the system, I am not able to
get results so far as I had to close R after waiting for a long time. I am
confused that the same code would work differently on the same system.
I tried out foreach package as well but didnt
1. Always cc to the list unless it is truly a private offlist reply.
This is to get help from a wider audience, as may well be required
here. Translation: Take my solution with a grain of salt. It is
fragile at best.
2. I think ?match.call and ?deparse are what you're looking for:
f -
1. This is not an R question. Post on a statistics or optimization list.
2. However, I would note that if what you mean by multimodalsolution
is that there is more than one value that gives the same optimal
value (within numerical error?) then, probably no: An optimization
function by definition
Hi Bert,
Thank you for the code.
However, I don't see what I am doing different, but my output is different.
I would much rather have output similar to yours where only the input
objects are returned (instead of the fn name, the encapsulating
parentheses, etc.):
That's because I screwed up! I gave you the wrong function, f instead of g .
Here's g:
g - function(...){
sapply(as.list(match.call())[-1],deparse)
}
and the example should now work.
Sheepishly,
-Bert
On Fri, Oct 18, 2013 at 10:21 AM, Dan Abner dan.abne...@gmail.com wrote:
Hi Bert,
I am using the ... argument to parmeterize a user define fn to accept
multiple input objects. I subsquently save all these data as a list.
Question: what is the best way to recover or extract the original object
names that were fed to the fn?
The following function, ellipsisInfo, returns
On Oct 18, 2013, at 1:12 AM, peter dalgaard wrote:
On Oct 18, 2013, at 08:37 , David Winsemius wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal
random sample with missing values and ran into
Many thanks Bert! You are a big help!!
On Fri, Oct 18, 2013 at 1:25 PM, Bert Gunter gunter.ber...@gene.com wrote:
That's because I screwed up! I gave you the wrong function, f instead of
g .
Here's g:
g - function(...){
sapply(as.list(match.call())[-1],deparse)
}
and the example
Yes, similar, but better, as match.call() will get unwanted named
arguments, too.
However, I do not understand the
substitute(...())
idiom. Would you care to explain it? (No is an acceptable answer!).
I would have expressed it as:
as.list(substitute(list(...)))[-1]
to convert the parse tree
On Oct 18, 2013, at 8:26 AM, Steven LeBlanc wrote:
On Oct 17, 2013, at 11:37 PM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate normal
random sample with
On Fri, 18 Oct 2013, Lauria, Valentina wrote:
Dear list,
I am using the hurdle model for modelling the habitat of rare fish
species. However I do get an error message when I try to model my data:
test_new1-hurdle(GALUMEL~ depth + sal + slope + vrm + lat:long + offset(log(haul_numb)),
dist
However, I do not understand the
substitute(...())
idiom. Would you care to explain it? (No is an acceptable answer!).
I don't completely understand it either, I treat it as an idiom. I saw it on
this list once.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original
On Oct 18, 2013, at 10:31 AM, David Winsemius wrote:
On Oct 18, 2013, at 1:12 AM, peter dalgaard wrote:
On Oct 18, 2013, at 08:37 , David Winsemius wrote:
On Oct 17, 2013, at 9:11 PM, Steven LeBlanc wrote:
Greets,
I'm using nlminb() to estimate the parameters of a multivariate
When the system locks up, what do you see in the Task Manager? Is it
consuming CPU and memory? On the example data you sent, you won't get
a match on the time since there is not match for the first entry in
df1 in the 'b' dataframe. This leads to an error that you are not
checking for. Have
when the system locks up it is consuming CPU to ~50%-60%, physical memory
is about ~33%.
I did tried on a subset of the data, e.g i in 1:100 and it works although
takes about 5-10 min.
Thanks for your help!
On Fri, Oct 18, 2013 at 12:49 PM, jim holtman jholt...@gmail.com wrote:
When the
Hi,
May be this helps:
dat1 - structure(list(Position = c(15L, 22L, 38L, 49L, 55L, 61L, 62L,
14L, 29L, 63L, 46L, 22L, 18L, 24L, 22L, 49L, 42L, 38L, 29L, 22L,
29L, 23L, 42L), Count = c(15L, NA, NA, 5L, NA, 17L, 18L, NA,
NA, NA, 8L, NA, 20L, NA, NA, 16L, 19L, NA, NA, NA, 13L, NA, 33L
)), .Names
Wanted to re-start this thread a bit, since I'm still not exactly sure
the best approach to my problem -- basically, the parameters I'm try
to make persistent are installation locations of a particular command
line program that is not installed along with an R package I'm working
on (GDAL, for
1. I do not recall saying any such thing.
2. HOWEVER, no matter. There is no choice. Follow Brian's advice.
-- Bert
On Fri, Oct 18, 2013 at 3:46 PM, Jonathan Greenberg j...@illinois.edu wrote:
Wanted to re-start this thread a bit, since I'm still not exactly sure
the best approach to my
On 13-10-18 6:46 PM, Jonathan Greenberg wrote:
Wanted to re-start this thread a bit, since I'm still not exactly sure
the best approach to my problem -- basically, the parameters I'm try
to make persistent are installation locations of a particular command
line program that is not installed
Hello all,
I'm not really sure how to search for this in google/Rseek so there is probably
a command to do it. I also know I could write an apply loop to find it but
thought I would ask all you lovely R gurus.
I have a very long vector (length=1855190) it looks something like this
Hi Earl
Unfortunately, the code works for me, i.e. indents _and_ displays the accented
vowels correctly.
Can you send me the output of the function call
libxmlVersion()
and also sessionInfo(), please?
D.
On 10/18/13 10:27 AM, Earl Brown wrote:
Thanks Duncan. However, now I can't get the
I have a very long vector (length=1855190) it looks something like this
......etc so it would be something equivalent of doing:
rep(c(1,2,3,4,5), c(10,30,24,65,3))
How can I find the index of where the step/jump is? For example using the
above I would
get an index of 0,
Hint: ?diff
-- Bert
On Fri, Oct 18, 2013 at 5:17 PM, Benton, Paul
hpaul.bento...@imperial.ac.uk wrote:
Hello all,
I'm not really sure how to search for this in google/Rseek so there is
probably a command to do it. I also know I could write an apply loop to find
it but thought I would ask
Thanks Duncan. However, now I can't get the Spanish and Portuguese accented
vowels to come out correctly and still keep the indents in the saved document,
even when I set encoding = UTF-8:
library(XML)
concepts - c(español, português)
info - c(info about español, info about português)
doc -
Merge has an option 'all', when set to 'TRUE' merge will return a
data.frame of every row in both data.frames. You can also specify all.x,
and all.y.
I do not suggest taking this path. Rather, standardize time columns for
all 15 data.frames with as.POSIXct()/as.POSIXlt() depending on the input
Hi,
Found a bug in the function when tested. So, try this (added one more line):
#Modified function
fun1 - function(dat,n) {
rl - rle(is.na(dat[,Count]))
indx -
which(is.na(dat[,Count]))[rep(rl$lengths[rl$values],rl$lengths[rl$values])==n]
lst1 -
Hi,
Try
dat - read.table(text=
a b c
0 89 4
1 93 3
0 88 5,sep=,header=TRUE)
datNew - within(dat,{a-a*(100/1);c - c*(100/5)})
#or
vec1 - c(1,100,5)
as.data.frame(t(t(dat*(100))/vec1))
#or
sweep(dat*100,2,vec1,`/`)
A.K.
Hi, i am trying to get all the colums of my data
On 19-10-2013, at 02:17, Benton, Paul hpaul.bento...@imperial.ac.uk wrote:
Hello all,
I'm not really sure how to search for this in google/Rseek so there is
probably a command to do it. I also know I could write an apply loop to find
it but thought I would ask all you lovely R gurus.
In addition to Bill's method, you may also use:
vec1 - rep(c(1,2,3,4,5), c(10,30,24,65,3))
c(0,which(diff(vec2)!=0))
#or
indx - cumsum(rle(vec2)$lengths)
c(0,indx[-length(indx)])
#Bill's method was found to be the fastest
vec3 - rep(vec1,1e4)
system.time( res - c(0,which(diff(vec3)!=0)))
#
Hi,
The code is based on what you pasted on the page. In the original post, you
used sep=, while the pasted data shows , as delimiter.
dat1 - subset(head(read.table(text=#Hogd/met, Temp, 005[M], Value
#Hogd/met, Difftemp, 051[M], Value
BA0+
1 MTEMP005 1 [deg.C]
2 MDTMP051 1 [deg.C]
EOH
891231,
64 matches
Mail list logo