Displaying 20 results from an estimated 417 matches similar to: "Rounding of problem with sum command in R"

2017 Aug 22

1

Rounding of problem with sum command in R

On 2017-08-22 9:26 AM, niharika singhal wrote:
> Hello I have a vector
> v=c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.1134463,0.2074027)
> when i do
> sum(v)
> or
> 0.08844446+0.1744455+0.1379778+0.1209769+0.1573065+0.1134463+0.2074027
> i am getting output as 1
No:? That's only the display:
> sum(v)-1
[1] 1.6e-07
????? hope this helps.? Spencer

2017 Aug 22

1

Rounding of problem with sum command in R

... and following up on Spencer's answer, see the "digits" argument of
?options and ?print.default.
Cheers,
Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Tue, Aug 22, 2017 at 7:30 AM, Spencer Graves
<spencer.graves

2017 Aug 26

1

Find maxima of a function

Please keep the list in cc.
Sorry, it didn't work as expected. Maybe someone else have an appropriate
solution.
Best,
Ulrik
On Sa., 26. Aug. 2017, 12:57 niharika singhal <niharikasinghal1990 at gmail.com>
wrote:
> Hi
>
> Thanks for you mail,
> I really appreciate your time on my problem
>
> I have posted this problem on
>
>
>

2017 Aug 26

2

Fwd: Find maxima of a function

Hi,
Thanks for your mail, and time
It is not working for some arguments, when mean value is like >6.
case
mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.
1134463,0.2074027)
rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133),
Norm(664.0713, 3.674773),
Norm(669.0484, 4.101381),

2017 Sep 20

1

How to use depmix for HMM with intial parameters

Hello,
I have initial parameters for HMM model and I want to use depmixS4 package.
The parameters are in the form
intial_prob_matrix=matrix(c(0.07614213, 0.45177665, 0.47208122), nrow=1,
ncol=3, byrow = TRUE)
transition_matrix=matrix(c(0.46666667,0.46666667,0.06666667,
0.06741573,0.5617978,0.37078652,
0.02173913,0.3478261,0.63043478), nrow = 3, ncol =

2017 Aug 26

1

Find maxima of a function

I have a Gaussian mixture model with some parameters
mean=(506.8644,672.8448,829.902)
sigma=(61.02859,9.149168,74.84682)
c=(0.1241933, 0.6329082, 0.2428986)
And the plot look something like below.[image: enter image description here]
<https://i.stack.imgur.com/4uUQ9.png>
Also, if I change my parameters to
mean=(2.15,2.0,2.9)
sigma=(0.1,0.1,0.1)
c=(1/3,1/3,1/3)

2017 Oct 10

1

"Time Series Plotting"

Hello,
I need some help in plotting time series.
I have dataframe Data with two column and thousands of row, I want wherever
the sequence corresponding to Energy column is changed the color change
should be reflected in Time Series plot, some rows of dataframe are below
Energy sequence
13536 3
14335 3
14638 3
25363 3
18511 2
18371 2
14555 3
I am able

2017 Aug 27

1

Fwd: Find maxima of a function

Dear Niharika,
As I said before, the problem is basically an optimization issue. You should isolate the problematic part from the rest of your study. Sometimes, more information does not help to solution. All the answers from us (Ulrik, David, me) are more or less are correct to find a maximum point. Newton?s method is also correct. But after answers, you only say, it didn?t find the right

2017 Aug 24

1

Problem in optimization of Gaussian Mixture model

Hello,
I am facing a problem with optimization in R from 2-3 weeks.
I have some Gaussian mixtures parameters and I want to find the maximum in
that
*Parameters are in the form *
mean1 mean2 mean3 sigma1 sigma2 sigma3 c1 c2 c3
506.8644 672.8448 829.902 61.02859 9.149168 74.84682 0.1241933
0.6329082 0.2428986
I have used optima and optimx to find the

2017 Aug 27

1

Fwd: Find maxima of a function

---------- Forwarded message ----------
From: niharika singhal <niharikasinghal1990 at gmail.com>
Date: Sun, Aug 27, 2017 at 11:57 AM
Subject: Re: Find maxima of a function
To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>,
"Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo
<ulrik.stervbo at gmail.com>

2017 Aug 27

1

Fwd: Find maxima of a function

I have not followed the history of this thread, but I am quite flummoxed as to why the OP is rewriting code to estimate parameters from an univariate Gaussian mixture model when alternatives such as EMCluster (which generally appears to handle initialization better than MClust) exist. Or perhaps there is more to it in which case I apologize. But I thought that I would make the OP aware of the

2017 Aug 28

1

"Improvement with the R code"

Hi,
I think you overthought this one a little bit, I don't know if this is the
kind of code you are expecting but I came up with something like that:
generate_transition_matrix <- function(data, n_states) {
#To be sure I imagine you should check n_states is right at this point
transitions <- matrix(0, n_states, n_states)
#we could improve a little bit here because at

2017 Aug 28

1

"Improvement with the R code"

Hello,
I am trying to implement a formula
aij= transition from state S_i to S_j/no of transition at state S_i
Code I have written is working with three state {1,2,3 }, but if the number
of states become={1,2,3,4,......n} then the code will not work, so can some
help me with this.
For and some rows of my data frame look like

2017 Aug 26

1

Find maxima of a function

> On 26 Aug 2017, at 16:39, niharika singhal <niharikasinghal1990 at gmail.com> wrote:
>
> Hi,
>
> Thanks for your mail, and time
>
> It is not working for some arguments, when mean value is like >6.
>
>
> case
>
> mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.
> 1134463,0.2074027)
>
> rv

2017 Aug 10

1

"Help On optim"

Hello,
I have some parameters from Mclust function. The parameters are in the form
*parametersDf *
* mu_1 mu_2 var_mc1 var_mc2 c1
c2 *
*2 1.357283 2.962736 0.466154 0.1320129 0.5258975
0.4741025 *
*21 8.357283 9.962736 0.466154 0.1320129 0.5258975
0.4741025 *
Each row in the above data frame

2002 Jul 08

4

How to freeze & thaw ext2(ext3) file system?

hi...
i intend to develop a backup application for linux
platform...
the workflow is as follows:
- freeze the file-system(ext2 or ext3)
- take the bckup
- release the file system.
i am not able to find the relevant ioctl or api for
freezing & releasing the file system...
so... i need some help...
thanks...
Satish Singhal
__________________________________________________
Do You Yahoo!?

2017 Aug 28

1

"Improvement with the R code"

Chuck (Is it fine to call you Chuck?) has far more R jutsu than I do
obviously.
I don't know much about pmin and factor but it might worth looking into if
you want to manipulate states by names (I assume this is why one might want
to use it?)
generate_transition_matrix <- function(data, states)
prop.table(table(head(data, -1), tail(data, -1)), 1)[states,]

2017 Aug 28

0

"Improvement with the R code"

Ok, I assumed you wanted to compute a matrix M for all states such that
M[i][j]= transition from state i to state j / number of transition from
state i
but from what you just answered it looks like you want to compute a matrix
M for a set of states S such that:
M[S_i][S_j]= transition from state S_i to state S_j / number of transition
from state S_i where S_i is in the set S
?

2005 Mar 12

3

How to use TIMECHECK tag in htb.init ?

Does anyone use tb.init TIMECHECK to change rate/ceil for classes based on periods over the day ? I had problems setting the right syntax. Can anyone help me with a quick example. I have tried the following with no success:
LEAF=sfq
RULE=xx.xx.xx.xx
TIME=06:00-18:00;128Kbit
TIME=18:00-06:00;256Kbit
RATE=196Kbit
_______________________________________________
LARTC mailing list

2017 Oct 11

0

"Time Series Plotting"

Hi
keep your emails to R help, I do not offer private consultance and others could have different opinion how to solve your problem.
Did you even try my suggestion? If not, why not? If yes in what respect it does not comply with your expectations.
Cheers
Petr
From: niharika singhal [mailto:niharikasinghal1990 at gmail.com]
Sent: Tuesday, October 10, 2017 2:34 PM
To: PIKAL Petr <petr.pikal at