text:
Home » Courses » Electrical Engineering and Computer Science » Discrete Stochastic Processes » Video Lectures » Lecture 8: Markov Eigenvalues and Eigenvectors
Lecture 8: Markov Eigenvalues and Eigenvectors
Description: This lecture covers eigenvalues and eigenvectors of the transition matrix and the steady-state vector of Markov chains. It also includes an analysis of a 2-state Markov chain and a discussion of the Jordan form.
Instructor: Prof. Robert Gallager
Lecture 1: Introduction and...
Lecture 2: More Review; The...
Lecture 3: Law of Large Num...
Lecture 4: Poisson (The Per...
Lecture 5: Poisson Combinin...
Lecture 6: From Poisson to ...
Lecture 7: Finite-state Mar...
Lecture 8: Markov Eigenvalu...
Lecture 9: Markov Rewards a...
Lecture 10: Renewals and th...
Lecture 11: Renewals: Stron...
Lecture 12: Renewal Rewards...
Lecture 13: Little, M/G/1, ...
Lecture 14: Review
Lecture 15: The Last Renewal
Lecture 16: Renewals and Co...
Lecture 17: Countable-state...
Lecture 20: Markov Processe...
Lecture 21: Hypothesis Test...
Lecture 22: Random Walks an...
Lecture 23: Martingales (Pl...
Lecture 24: Martingales: St...
Lecture 25: Putting It All ...
The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high-quality educational resources for free. To make a donation or view additional materials from hundreds of MIT courses, visit MIT OpenCourseWare at ocw.mit.edu.
PROFESSOR: OK, let's get started again on finite-state Markov chains. Sorry I was away last week. It was a long-term commitment that I had to honor. But I think I will be around for all the rest of the lectures. So I want to start out by reviewing just a little bit. I'm spending a lot more time on finite-state Markov chains than we usually do in this course, partly because I've rewritten this section, partly because I think the material is very important. It's sort of bread-and-butter stuff, of discrete stochastic processes. You use it all the time. It's a foundation for almost everything else.
And after thinking about it for a long time, it really isn't all that complicated. I used to think that all these details of finding eigenvalues and eigenvectors and so on was extremely tedious. And it turns out that there's a very nice pleasant theory there. You can find all of these things after you know what you're doing by very simple computer packages. But they don't help if you don't know what's going on. So here, we're trying to figure out what's going on.
So let's start out by reviewing what we know about ergodic unit chains and proceed from there. An ergodic finite-state Markov chain has transition probabilities which, if you look at the transition matrix raised to the nth power, what that gives you is the transition probabilities of an n-step Markov chain. In other words, you start at time 0, and at time n, you look at what state you're in. P sub ij to the nth power is then the probability that you're in state j at time n, given that you're in state i at time 0. So this has all the information that you want about what happens to Markov chain as time gets large.
One of the things we're most concerned with is, do you go to steady state? And if you do go to steady state, how fast do you go to steady state? And of course, this matrix tells you the whole story there, because if you go to steady state, and the Markov chain forgets where it started, then P sub ij to the n goes to some constant, pi sub j, which is independent of the starting state, i, and independent of m, asymptotically, as n gets big. So this pi is a strictly positive probability vector. I shouldn't say so it is. That's something that was shown last time.
If you multiply both sides of this equation by P sub jk in sum over k, then what do you get? You get P sub ik to the n plus 1. That goes to a limit also. If the limit in n goes to infin-- then the limit as n plus 1 goes to infinity is clearly the same thing. So this quantity here is the sum over j, of pi sub j, P sub jk. And this quantity is equal to pi sub k, just by definition of this quantity. So P sub k is equal to sum of pi j.
Pjk, what does that say? That's the definition of a steady state vector. That's the definition of, if your probabilities of being in state k satisfy this equation, then one step later, you still have the same probability of being in state k. Two steps later, you still have the same probability of being in state k. So this is called the steady state equation. And a solution to that is called a steady state vector. And that satisfies this.
In matrix terms, if you rate this out, what does it say? It says the limit as n approaches infinity of p to the n is equal to the column vector, e of all 1s. The transpose here means it's a column vector. So you have a column vector times a row vector. Now, you know if you have a row vector times a column vector, that just gives you a number. If you have a column vector times a row vector, what happens? Well, for each element of the column, you get this whole row. And for the next element of the column, you get the whole row down beneath it multiplied by the element of the column, and so forth, day on.
So a column vector times a row vector is, in fact, a whole matrix. It's a j by j matrix. And since e is all 1s, what that matrix is is a matrix where every row is a steady state vector pi. So we're saying not only does this pi that we're talking about satisfy this steady state equation, but more important, it's this limiting vector here. And as n goes to infinity, you in fact do forget where you were. And the entire matrix of where you are at time n, given where you were at time 0, goes to just this fixed vector pi.
So this is a column vector, and pi is a row vector then. The same result almost holds for ergodic unit chains. What's an ergodic unit chain? An ergodic unit chain is an ergodic set of states plus a whole bunch of transient states. Doesn't matter whether the transient states are one class of transient states or whether it's multiple classes of transient states. It's just transient states. And there's one recurrent class. And we're assuming here that it's recurrent. So you can almost see intuitively that if you start out in any one of these transient states, you bum around through the transient states for a while. And eventually, you flop off into the recurrent class. And once you're in the recurrent class, there's no return. So you stay there forever.
Now, that's something that has to be proven. And it's proven in the notes. It was probably proven last time. But anyway, what happens then is that the sole difference between ergodic unit chains and just having a completely ergodic Markov chain is that the steady state factor is now positive for all ergodic states and it's 0 for all transient states. And aside from that, you still get the same behavior still. As n gets large, you go to the steady state vector, which is the steady state vector of the ergodic chain.
If you're doing this stuff by hand, how do you do it? Well, you start out just with the ergodic class. I mean, you might as well ignore everything else, because you know that eventually you're in that ergodic class. And you find the steady state vector in that ergodic class, and that's the steady state vector you're going to wind up with. This is one advantage of understanding what you're doing, because if you don't understand what you're doing and you're just using computer programs, then you never have any idea what's ergodic, what's not ergodic or anything else. You just plug it, you grind away, you get some answer and say, ah, I'll publish a paper. And you put down exactly what the computer says, but you have no interpretation of it at all.
So the other way of looking at this is, when you have a bunch of transient states, and you also have an ergodic class, you can represent a matrix if the recurrent states are at the end of the chain and the transient states are at the beginning of the chain. This matrix here is the matrix of transition probabilities within the recurrent class. These are the probabilities for going from the transient states to the recurrent class. And once you get over here, the only place to go is down here. And the transient class is just a t by t class. And the recurrent class is just a j minus t by j minus t matrix.
So the idea is that each transient state eventually has a transition to a recurrent state, and the class of recurrent states leads to study state as before. So that really, all that analysis of ergodic unit chains, if you look at it intuitively, it's all obvious. Now, as in much of mathematics, knowing that something is obvious does not relieve you of the need to prove it, because sometimes you find that something that looks obvious is true most of the time but not all of the time. And that's the purpose of doing these things.
There's another way to express this eigenvalue, eigenvector equation we have here. And that is that the transition matrix minus lambda times the identity matrix times the column vector v is equal to 0. That's the same as the equation p times v is equal to v. That's the same as a right eigenvector. Well, this is the equation for an eigenvalue 1. This is an equation for an arbitrary eigenvalue lambda. But p times v equals lambda times v is the same as p minus lambda i times v equals 0.
Why do we even bother to say something so obvious? Well, because when you look at linear algebra, how many of you have never studied any linear algebra at all or have only studied completely mathematical linear algebra, where you never deal with n-tuples as vectors or matrices or any things like this? Is there anyone? If you don't have this background, pick up-- what's his name?
AUDIENCE: Strang.
PROFESSOR: Strang. Strang's book. It's a remarkably simple-minded book which says everything as clearly as it can be stated. And it tells you everything you have to know. And it does it in a very straightforward way. So I highly recommend it to get any of the background that you might need. Most of you, I'm sure, are very familiar with these things. So I'm just reminding you of then. Now, a square matrix is singular if there's a vector v, such that a times v is equal to 0. That's just a definition as a singularity.
Now, lambda is an eigenvalue of a matrix p if and only if p minus lambda times i is singular. In other words, if there's some v for which p minus lambda i times v is equal to 0, that's what this says. You put p minus lambda i in for a, and it says it's singular if there's some v for which this matrix-- this matrix is singular if there's some v such that p minus lambda i times v is equal to 0.
So let a1 to am be the columns of a. Then a is going to be singular if a1 to am are linearly dependent. In other words, if there's some set of coefficients you can attach to a1 times v1 plus a2 times v2, plus up to am times vm such that that sum is equal to 0, that means that a1 to am are linearly dependent. It also means that the matrix a times that v is equal to 0. So those two things say the same thing again. So the square matrix, a, is singular if and only if the rows of a are linearly independent. We set columns here. Here, we're doing the same thing for rows. It still holds true. And one new thing, if and only if the determinant of a is equal to 0.
One of the nice things about determinants is that determinants are 0 if the matrix is singular, if and only if the matrix is singular. So the summary of all of this for a matrix which is a transition matrix-- namely, a stochastic matrix-- is lambda, is an eigenvalue of p, if and only if p minus lambda i is singular, if and only if the determinant of p minus lambda i is equal to 0, if and only if p times some vector v equals lambda v, and if and only if u times p equals lambda u for some u. Yes?
AUDIENCE: The second to last statement is actually linearly independent, you said? The second to last. Square matrix a. No, above that.
PROFESSOR: Oh, above that. A square matrix a is singular if and only if the rows of a are linearly dependent, yes.
AUDIENCE: Dependent.
PROFESSOR: Dependent, yes. In other words, if there's some vector v such that a times v is equal to 0, that means that those columns are linearly dependent. So we need all of those relationships. It says for every stochastic matrix-- oh, now this is something new. For every stochastic matrix, P times e is equal to e. Obviously, because if you sum up the sum of Pij over j is equal to 1. P sub ij is the probability, given that you start in state i, that in the next step, you'll be in state j. You have to be somewhere in the next step. So if you sum these quantities up, you have to get 1, which says you have to be some place.
So that's all this is saying. That's true for every finite-state Markov chain in the world, no matter how ugly it is, how many sets of recurrent states it has, how much periodicity it has. A complete generality, P times e is equal to e. So lambda is always an eigenvalue of a stochastic matrix, and e is always a right eigenvector. Well, from what we've just said, that means there has to be a left eigenvector also. So there has to be some pi such that pi times P is equal to pi.
So suddenly, we find there's also a left eigenvector. What we haven't shown yet is that that pi that satisfies this equation is a probability vector. Namely, we haven't shown that all the components of pi are greater than or equal to 0. We still have to do that. And in fact, that's not completely trivial. If we can find such a vector that is a probability vector, the compound in sum to 1 and they're not negative, then this is the equation for a steady state vector. So what we don't know yet is whether a steady state vector exists. We do know that a left eigenvector exists.
We're going to show later that there is a steady state vector pi. In other words, a non-negative vector which sums to 1 for all finite-state Markov chains. In other words, no matter how messy it is, just like e, the column vector of all 1s is always a right eigenvector of eigenvalue 1. There is always a non-negative vector pi whose components sum to 1, which is a left eigenvector with eigenvalue 1. So these two relationships hold everywhere.
Incidentally, the notes at one point claim to have shown this. And the notes really don't show it. I'm going to show it to you today. I'm sorry for that. It's something I've known for so long that I find it hard to say is this true or not. Of course it's true. But it does have to be shown, and I will show it to you later on. Chapter three of the notes is largely rewritten this year. And it has a few more typos in it than most of the other chapters. And a few of the typos are fairly important. I'll try to point some of them out as we go. But I'm sure I haven't caught them all yet.
Now, what is the determinant of an M by M matrix? It's this very simple-looking but rather messy formula, which says the determinant of a square matrix A is the sum over all partitions-- and then there's a plus minus here, which I'll talk about later-- of the product from i equals 1 to M. M is the number of states of A sub i. This is the component of the ij position. And we're taking A sub i. And then the partition that we're dealing with, mu sub i.
So what we're doing is taking a matrix with all sorts of terms in it-- A11 up to A1j on to Aj1 up to A sub jj. And these partitions we're talking about are ways of selecting one element from each row and one element from each column. Namely, that first sum there is talking about one element from each row. And then when we're talking about a permutation here, we're doing something like, for this row, we're looking at, say, this element. For this row, we might be looking at this element. For this row, we might be looking at this element, and so forth down, until finally, we're looking at some element down here.
Now, we've picked out every column and every row in doing this, but we only have one element in each row and one element in each column. If you've studied linear algebra and you're at all interested in computation, the first thing that everybody tells you is that this is a god-awful way to ever compute a determinant, because the number of permutations grows very, very fast with the size of the matrix. And therefore you don't want to use this formula very often. It's a very useful formula conceptually, though, because if we look at the determinant of p minus lambda i, if we want to ask the question, how many eigenvalues does this transition matrix have? well, the number of eigenvalues it has is the number of values of lambda such that the determinant of p minus lambda i is 0.
Now, how many such values are there? Well, you look the matrix for that, and you get A11 minus lambda A12 and A22 minus lambda Ajj minus lambda. And none of the other elements have lambda in it. So when you're looking at this formula for finding the determinant, one of the partitions is this partition, which is a polynomial of degree j in lambda. All of the others are polynomials of degree less than j in lambda. And therefore this whole bloody mess here is a polynomial of degree j and lambda.
So the equation, determinant of p minus lambda i, which is a polynomial of degree j in lambda, equals 0. How many roots does it have? Well, the fundamental theorem of algebra says that a polynomial of degree j, of complex numbers-- and real is a special case of complex-- that it has exactly j roots. So there are exactly, in this case, M-- excuse me, I've been calling it j sometimes and M sometimes. This equation here has exactly M roots to it. And since it has exactly M roots, that's the number of eigenvalues there are.
There's one flaw in that argument. And that is, some of the roots might be repeated. Say you have M roots altogether. Some of them appear more than one time, so you'll have roots of multiplicity, something or other. And when you add up the multiplicities of each of the distinct eigenvalues, you get capital M, which is the number of states. So the number of different eigenvalues is less than or equal to M. And the number of distinct eigenvalues times the multiplicity of each eigenvalue is equal to M.
That's a simple, straightforward fact. And it's worth remembering. So there are M roots to the equation. Determinant p minus lambda i equals 0. And therefore there are M eigenvalues of p. And therefore you might think that there are M eigenvectors. That, unfortunately, is not true necessarily. That's one of the really-- it's probably the only really ugly thing in linear algebra. I mean, linear algebra is a beautiful theory.
I mean, it's like Poisson's stochastic processes. Everything that can be true is true. And if something isn't true, there's a simple counter-example of why it can't be true. This thing is just a bloody mess. But unfortunately, if you have M states in a finite-state Markov chain, you might not have M different eigenvectors. And that's unfortunate, but we will forget about that for as long as we can, and we'll finally come back to it towards the end.
AUDIENCE: [INAUDIBLE]?
PROFESSOR: What?
AUDIENCE: Why would we care about all the eigenvectors if we are only concerned with the ones that [INAUDIBLE]?
PROFESSOR: Well, because we're interested in the other ones because that tells us how fast p to the M converges to what it should be. I mean, all those other eigenvalues, as we'll see, are the error terms in p to the M as it approaches this asymptotic value. And therefore we want to know what those eigenvalues are. At least we want to know what the second-biggest eigenvalue is. Now, let's look at just a case of two states. Most of the things that can happen will happen with two states, except for this ugly thing that I told you about that can't happen with two states. And therefore two states is a good thing to look at, because with two states, you can calculate everything very easily and you don't have to use any linear algebra.
So if we look at a Markov chain with two states, P sub ij is this set of transition probabilities. The left eigenvector equation is pi 1 times P11 times pi 2 times P21 is equal to lambda pi 1. And so this is writing out what we said before. The vector pi times the matrix P is equal to lambda times the vector pi. That covers both of these equations. Since M is only 2, we only have to write things out twice. Same thing for the right eigenvector equation. That's this.
The determinant of P minus lambda i, if we use this formula that we talked about here, you put A11 minus lambda, A22 minus lambda. Well, then you're done. So all you need is P11 minus lambda times P22 minus lambda. That's this permutation there. And then you have an odd permutation, A12 times A21. How do you know which permutations are even and which permutations are odd?
It's how many flips you have to do. But to see that that's consistent, you really have to look at Strang or some book on linear algebra, because it's not relevant here. But anyway, that determinant is equal to this quantity here. That's a polynomial of degree 2 in lambda. If you solve it, you find out that one solution is lambda 1 equals 1. The other solution is lambda 2 is 1 minus P12 minus P21.
Now, there are a bunch of cases to look at here. If the off-diagonal transition probabilities are both 0, what does that mean? It means if you start in state 0, you stay there. If you start in state 1, you stay there forever. If you start in state 2, you stay there forever. That's a very boring Markov chain, but it's not very nice for the theory. So we're going to leave that case out for the time being. But anyway, if you have that case, then the chain has two recurrent classes. Lambda equals 1, has multiplicity 2. You have two eigenvalues of algebraic multiplicity 2. I mean, it's just one number, but it appears twice in this determinant equation. And it also appears twice in the sense that you have two recurrent classes.
And you will find that there are two linearly independent left eigenvectors, two linearly independent right eigenvectors. And how do you find those? You use your common sense and you say, well, if you start in state 1, you're always there. If you start in state 2, you're always there. Why do I even look at these two states? This is a crazy thing where wherever I start, I stay there and I only look at state 1 or state 2. It's scarcely even a Markov chain.
If P12 and P21 are both 1, what it means is you can never go from state 1 to state 1. You always go from state 1 to state 2. And you always go from state 2 to state 1. It means you have a two-state periodic chain. And that's the other crazy case. The other case is not very interesting. There's nothing stochastic about it at all. So the chain is periodic. And if you look at this equation here, the second eigenvalue is equal to minus 1.
I might as well tell you that, in general, if you have a periodic Markov chain, just one recurrent class and it's periodic, a period d, then the eigenvalues turn out to be the uniformly spaced eigenvalues around the unit circle. One is one of the eigenvalues. We've already seen that. And the other d minus 1 eigenvalues are those uniformly spaced around the unit circle. So they add up to 360 degrees when you get all done with it.
So that's an easy case. Proving that is tedious. It's done in the notes. It's not even done in the notes. It's done in one of the exercises. And you can do it if you choose. So let's look at these eigenvector equations and the eigenvalue equations. Incidentally, if you don't know what the eigenvalues are, is this a linear set of equations? No, it's a nonlinear set of equations. This is a nonlinear set of equations in pi 1, pi 2, and lambda.
How do you solve non-linear equations like that? Well, if you have much sense, you first find out what lambda is and then you solve linear equations. And you can always do that. We've said that these solutions for lambda, there can only be M of them. And you can find them by solving this polynomial equation. Then you can solve the linear equation by finding the eigenvectors. There are packages to do all of these things, so there's nothing you should waste time on doing here. It's just knowing what the results are that's important.
From now on, I'm going to assume that P12 or P21 are greater than 0. In other words, I'm going to assume that we don't have the periodic case and we don't have the case where you have two classes of states. In other words, I'm going to assume that our Markov chain is actually ergodic. That's the assumption that I'm making here. If you then solve these equations using lambda 1 equals 1, you'll find out that pi 1 is the component sum to 1. First component is P21 over the sum. Second component is P12 over the sum. Not very interesting.
Why is the steady state probability weighted towards the largest of these transition probabilities? If P21 is bigger than P12, how do you know intuitively that you're going to be in state 1 more than you're in state 2? Is this intuitively obvious to-- yeah?
PROFESSOR: Because you make more transistors from 2 to 1. Well, actually you don't make more transitions from 2 to 1. You make exactly the same number of transitions, but since the probability is higher, it means you have to be in state 1 more of the time. Good. So these are the two. And this is the left eigenvector for the second eigenvalue-- namely, the smaller eigenvalue. Now, if you look at these equations, you'll notice that the vector pi, the left i-th eigenvector, multiplied by the right j-th eigenvector, is always equal to delta ij. In other words, the left eigenvectors are orthogonal to the right eigenvectors.
I mean, you can see this just by multiplying it out. You multiply pi 1 times nu 1, and what do you get? You get this plus this, which is 1. Delta 11 means there's something which is 1 when i is equal j and 0 when i is unequal to j. You take this and you multiply it by this, and what do you get? You get P21 times P12 over the square. Minus P12 times P21, it's 0. Same thing here. 1 minus 1, that vector times this vector, is 0 again. So the cross-terms are 0. The diagonal terms are 1. That's the way it is.
So let's move on with this. These right eigenvector equations, you can write them in matrix form. I'm doing this slowly. I hope I'm not boring those who have done a lot of linear algebra too much. But they won't go on forever, and it gets us to where we want to go. So if you take these two equations and you write them in matrix form, what you get is P times u, where u is a matrix whose columns are the vector nu 1 and the vector nu 2. And capital lambda is the diagonal matrix of the eigenvalues. If you multiply P times the first column of u, and then you look at the first column of this matrix, what you get-- yes, that's exactly the right way to do it.
And if you're not doing that, you're probably not understanding it. But if you just think of ordinary matrix vector multiplication, this all works out. Because of this orthogonality relationship, we see that the matrix whose rows are the left eigenvectors times the matrix whose columns are the right eigenvectors, that's equal to i. Namely, it's equal to the identity matrix. That's what this orthogonality relationship means. This means that this matrix is the inverse of this matrix. This proves that u is invertible. And in fact, we've done this just for m equals 2. But in fact, this proof is general and holds for arbitrary Markov chains if the eigenvectors span the space. And we'll see that later.
We're doing this for m equals 2 now, so we how to proceed when we have an arbitrary Markov chain. u is invertible. u to the minus 1 has pi 1 and pi 2 as rows. And thus P is going to be equal to-- I guess we should-- oh, we set it up here. P times u is equal to u times lambda. We've shown here that u is invertible, therefore we can multiply this equation by u to the minus 1. And we get the transition matrix P is equal to u times the diagonal matrix lambda times the matrix u to the minus 1.
What happens if we try to find P squared? Well, it's u times lambda times u to the minus 1. One of the nice things about matrices is you can multiply them, if you don't worry about the details, almost like numbers. Times u times lambda times u to the minus 1. Except you don't have commutativity. That's the only thing that you don't have. But anyway, you have u times lambda times u to the minus 1 times u times lambda times u to t he minus 1. This and this turn out to be the identity matrix, so you have u times lambda times lambda, which is lambda squared, times u to the minus 1. You still have this diagonal matrix here, but the eigenvalues have all been doubled.
If you keep doing that repeatedly, you find out that P to the n-- namely, this long-term transition matrix, which is the thing we're interested in-- is the matrix u times this diagonal matrix, lambda to the n, times u to the minus 1. Equation 329 in the text has a typo, and it should be this. It's given as u to the minus 1 times lambda to the n times u, which is not at all right. That's probably the worst typo, because if you try to say something from that, you'll get very confused. You can solve one in general if all the M eigenvalues are distinct as easily as for M equals 2. This is still valid so long as the eigenvectors span the space.
So now the thing we want to do is relatively simple. This lambda to the n is a diagonal matrix. I can represent it as the sum of M different matrices. And each of those matrices has only one diagonal element, non-0. In other words, for the case here, what we're doing is taking lambda 1, 0 to the n, 0 lambda 2 to the n, and representing this as lambda 1 to the n, 0, 0, 0, plus 0, 0, 0 lambda 2 to the n.
So we have those trivial matrices with u on the left side and u to the minus 1 on the right side. And we think of how to multiply the matrix u, which is a matrix whose columns are the eigenvectors, times this matrix with only one non-0 element, times the matrix here, whose elements are the left eigenvectors. And how do you do that? Well, if you do this for a while, and you think of what this one element here times a matrix whose rows are eigenvectors does, this non-0 term in here picks out the appropriate row here. And this non-0 element picks out the appropriate column here.
So what that gives you is p to the n is equal to the sum over the number of states in the Markov chain times lambda sub i-- the i-th value to the nth power-- times nu to the i times pi to the i. pi to the i is the i-th eigenvector of p. nu to the i is the i-th right eigenvector of p. They have nothing to do with n. The only thing that n affects is this eigenvalue here. And what this is saying is that p to the n is just the sum of eigenvalues which are, if lambda is bigger than 1, this is exploding. If lambda 1 is less than 1, it's going to 0. And if lambda 1 is equal to 1, it's staying constant. If lambda 1 is complex but has magnitude 1, then it's just gradually rotating around and not doing much of interest at all, but it's going away.
So that's what this equation means. It says that we've converted the problem of finding the nth power of p just to this problem of finding the nth power of these eigenvalues. So we've made some real progress.
AUDIENCE: Professor, what is nu i right here?
AUDIENCE: What is nu i?
PROFESSOR: nu sub i is the i-th of the right eigenvectors of the matrix p.
AUDIENCE: And pi i?
PROFESSOR: And pi i is the i-th left eigenvector. And what we've shown is that these are orthogonal to each other, orthonormal.
AUDIENCE: Can you please say again what happens when lambda is complex?
AUDIENCE: When lambda is complex, what exactly happens?
PROFESSOR: Oh, if lambda i is complex and the magnitude is less than 1, it just dies away. if the magnitude is bigger than 1, it explodes, which will be very strange. And we'll see that can't happen. And if the magnitude is 1, as you take powers of a complex number of magnitude 1, I mean, it start out here, it goes here, then here. I mean, it just rotates around in some crazy way. But it maintains its magnitude as being equal to 1 all the time.
So this is just repeating what we had before. These are the eigenvectors. If you calculate this very quickly using this and this, and if you recognize that the right eigenvector, nu 2, is the first part of it is pi sub 2, the second part of it minus pi sub 1, where pi is just this first eigenvector here. So if you do this multiplication, you find that nu to the 1-- oh, I thought I had all of these things out. This should be nu.
The first right eigenvector times the first left eigenvector. Oh, but this is all right, because I'm saying the first left eigenvector is a steady state vector, which is the thing we're interested in. That's pi 1, pi 2, pi 1, pi 2, where pi 1 is this and pi 2 is this. nu 2 times pi 2 is just this. So when we calculate p sub n, we get pi 1 plus pi 2 times this eigenvalue to the nth power. Pi 1 minus pi 1, lambda 2 to the nth power. pi 2 and pi 2 is what we get for the main eigenvalue. This is what we get for the little eigenvalue. This little eigenvalue here is 1 minus P12 minus P21, which has magnitude less than 1, unless we either have the situation where P12 is equal to P21 is equal to 0, or both of them are 1.
So these are the terms that go to 0. This solution is exact. There were no approximations in here. Before, when we analyzed what happened to P to the n, we saw that we converged, but we didn't really see how fast we converged. Now we know how fast we converge. The rate of convergence is the value of this second eigenvalue here. And that's a pretty general result. You converged like the second-largest eigenvalue. And we'll see how that works out.
Now, let's go on to the case where you have an arbitrary number of states. We've almost solved that already, because as we were looking at the case with two states, we were doing most of the things in general. If you have an n state Markov chain, a determinant of P minus lambda is a polynomial of degree M in lambda. That was what we said a while ago. It has M roots, eigenvalues. And here, we're going to assume that those roots are all distinct. So we don't have to worry about what happens with repeated roots. Each eigenvalue lambda sub i-- there are M of them now-- has a right eigenvector, nu sub i, and a left eigenvector, pi sub i. And we have seen that-- well, we haven't seen it yet. We're going to show it in a second.
pi super i times nu super j is equal to j for each ij unequal to i. If you scale either this or that, when you saw this eigenvector equation, you have a pi on both sides or a nu on both sides, and you have a scale factor which can't be determined from the eigenvector equation. So you have to choose that scaling factor somehow. If we choose the scaling factor appropriately, we get pi, the i-th left eigenvector, times the i-th right eigenvector. This is just a number now. It's that times that. We can scale things, so that's equal to 1.
Then as before, let u be the matrix with columns nu 1 to nu M, and let v have the rows, pi 1 to pi M. Because of this orthogonality relationship we've set up, v times u is equal to i. So again, the left eigenvector rows forms a matrix which is the inverse of the right eigenvector columns. So that says v is equal to u to the minus 1. Thus the eigenvector is nu, the first right eigenvector up to the nth right eigenvector, these are linearly independent. And they span M space.
That's a very peculiar thing we've done. We've said we have all these M right eigenvectors. We don't know anything about them, but what we do know is we also have M left eigenvectors. And the left eigenvectors, as we're going to show in just a second, are orthogonal to the right eigenvectors. And therefore, when we look at these two matrices, we can multiply them and get the identity matrix. And that means that the right eigenvectors have to be-- when we look at the matrix of the right eigenvectors, is non-singular.
Very, very peculiar argument. I mean, we find out that those right eigenvectors span the space, not by looking at the right eigenvectors, but by looking at how they relate to the left eigenvectors. But anyway, that's perfectly all right. And so long as we can show that we can satisfy this orthogonality condition, then in fact all this works out. v is equal to u to the minus 1. These eigenvectors are linearly independent and they span M space. Same here. And putting these equations together, P times u equals u times lambda. This is exactly what we did before. Post-multiplying by u to the minus 1, we get P equals u times lambda times u to the minus 1. P to the n is then u times lambda to the n times u to the minus 1.
All this stuff about convergence is all revolving down to simply the question of what happens to these eigenvalues. I mean, there's a mess first, finding out what all these right eigenvectors are and what all these left eigenvectors are. But once you do that, P to the n is just looking at this quantity, breaking up lambda to the n the way we did before. P to the n is just this sum here. Now, each row of P sums to 1, so e is a right eigenvector of eigenvalue 1. So we have a theorem that says the left eigenvector pi of eigenvalue 1 is a steady state vector if it's normalized to pi times e equals 1.
So we almost did that before, but now we want to be a little more careful about it. Oh, excuse me. The theorem is that the left eigenvector pi is a steady state vector if it's normalized in this way. In other words, we know that there is a left eigenvector pi, which has eigenvalue 1, because there's a right eigenvector. If there's a right eigenvector, there has to be a left eigenvector. What we don't know is that pi actually has non-negative terms. So that's the thing we want to show.
The proof is, there must be a left eigenvector pi for eigenvalue 1. We already know that. For every j, Pi sub j is equal to the sum over k times pi sub k times p sub kj. We don't know whether these are complex or real. We don't know whether they're positive or negative, if they're real. But we do know that since they satisfy this eigenvector equation, they satisfy this equation. If I take the magnitudes of all of these things, what do I get?
The magnitude on this side is pi sub j magnitude. This is less than or equal to the sum of the magnitudes of these terms. If you take two complex numbers and you add them up, you get something which, in magnitude, is less than or equal to the sum of the magnitudes. It might sound strange, but if you look in the complex plane-- imaginary, real-- and you look at one complex number, and you add it to another complex number, this distance here is less than or equal to this magnitude plus this magnitude. That's all that equation is saying.
And this is equal to this distance plus this distance if and only if each of these components of the eigenvector that we're talking about, if and only if those components are all heading off in the same direction in the complex plane. Now what do we do? Well, you look at this for a while and you say, OK, what happens if I sum this inequality over j? Well, if I sum this over j, I get one. And therefore when I sum both sides over j, the sum over j of the magnitudes of these eigenvector components is less than or equal to the sum over k of the magnitude. This is the same as this. This j is just a dummy index of summation. This is a dummy index of summation.
Obviously, this is less than or equal to this. But what's interesting here is that this is equal to this. And the only way this can be equal to this is if every one of these things are satisfied with equality. If any one of these are satisfied with inequality, then when you add them all up, this will be satisfied with inequality also, which is impossible. So all of these are satisfied with equality, which says that the magnitude of pi sub j, the vector whose elements are the magnitudes of this thing we started with, in fact form a steady state vector if we normalize them to 1.
It says these magnitudes satisfy the steady state equation. These magnitudes are real and they're positive. So when we normalize them to sum to 1, we have a steady state vector. And therefore the left eigenvector pi of eigenvalue 1 is a steady state vector if it's normalized to pi times e equals 1, which is the way we want to normalize them. So there always is a steady state vector for every finite-state Markov chain. So this is a non-negative vector satisfying a steady state vector equation. And normalizing it, we have a steady state vector. So we've demonstrated the existence of a left eigenvector which is a steady state vector.
Another theorem is that every eigenvalue satisfies lambda, magnitude of the eigenvalue is less than or equal to 1. This, again, is sort of obvious, because if you have an eigenvalue which is bigger than 1 and you start taking powers of it, it starts marching off to infinity. Now, you might say, maybe something else is balancing that. But since you only have a finite number of these things, that sounds pretty weird. And in fact, it is. So the proof of this is, we want to assume that pi super l is the l-th of these eigenvectors of P. Its eigenvalue is lambda sub l. It also is a left eigenvector of P to the n with eigenvalue lambda to the n. That's what we've shown before.
I mean, you can multiply this matrix P, and all you're doing is just taking powers of the eigenvalue. So if we start out with lambda to the n, let's forget about the l's, because we're just looking at a fixed l now. Lambda to the nth power times the j-th component of pi is equal to the sum over i of the i-th component of pi times Pij to the n, for all j.
Now I take the magnitude of everything is before. The magnitude of this is, again, less than or equal to the magnitude of this. I want to let beta be the largest of these quantities. And when I put that maximizing j in here, lambda to the l times beta is less than or equal to the sum over i of-- I can upper-bound these by beta. So I wind up with lambda to the l times beta is less than or equal to the sum over i of beta times Pij to the n.
I don't know what these powers are, but they're certainly less than or equal to 1. So lambda sub l is less than or equal to n. That's what this said. When you take this magnitude of the l-th eigenvalue, it's less than or equal to this number n. Now, if this number were larger than 1, if it was 1 plus 10 to the minus sixth, and you multiplied it by a large enough number n, that this would grow to be arbitrarily large. It can't grow to be arbitrarily large, therefore the magnitude of lambda sub l has to be less than or equal to 1. Tedious proof, but unfortunately, the notes just assume this.
Maybe I had some good, simple reason for it before. I don't have any now, so I have to go through a proof. Anyway, these two theorems, if you look at them, are valid for all finite-state Markov chains. There was no place that we used the fact that we had anything with distinct eigenvalues or anything. But now when we had distinct eigenvalues, we have the nth power of P is the sum here again over right eigenvectors times left eigenvectors. When you take a right eigenvector, which is a column vector, times a left eigenvector, which is a row vector, you get an M by M matrix. I don't know what that matrix is, but it's a matrix. It's a fixed matrix independent of n. And the only thing that's varying with n is these eigenvalues.
These quantities are less than or equal to 1. So if the chain is an ergodic unit chain, we've already seen that one eigenvalue is 1, and the rest of the eigenvalues are strictly less than 1 in magnitude. We saw that by showing that for an ergodic unit chain, P to the n converged. So the rate at which P to the n approaches e times pi is going to be determined by the second-largest eigenvalue in here. And that second-largest eigenvalue is going to be less than 1, strictly less than 1. We don't know what it is. Before, we knew this convergence here for an ergodic unit chain is exponential. Now we know that it's exponential and we know exactly how fast it goes, because the speed of convergence is just the second-largest eigenvalue.
If you want to know how fast P to the n approaches e times the steady state vector pi, all you have to do is find that second-largest eigenvalue, and that tells you how fast the convergence is, except for calculating these things, which are just fixed. If P is a periodic unit chain with period d, then if you read the notes-- you should read the notes-- there are d eigenvalues equally spaced around the unit circle. P to the n doesn't converge. The only thing you can say here is, what happens if you look at P to the d-th power? And you can imagine what happens if you look at P to the d-th power without doing any analysis.
I mean, we know that what happens in a periodic chain is that you rotate from one set of states to another set of states to another set of states to another set of states, and then back to the set of states you started with. And you keep rotating around. Now, there are d sets of states going around here. What happens if I take P to the d? P to the d is looking at the d-step transitions. So it's looking at, if you start here, after d steps, you're back here again, after d steps, you're back here again. So the matrix, P to the d, is in fact the matrix of d ergodic subclasses.
And for each one of them, whatever subclass you start in, you stay in that subclass forever. So the analysis of a periodic unit chain, really the classy way to do it is to look at P to the d and see what happens there. And you see that you get convergence within each subclass, but you just keep rotating among subclasses. So there's nothing very fancy going on there. You just rotate from one subclass to another. And that's the way it is. And P to the n doesn't converge. But P to the d times n does converge.
Now, let's look at the next-most complicated state. Suppose we have M states and we have M independent eigenvectors. OK, remember I told you that there was a very ugly thing in linear algebra that said, when you had an eigenvalue of multiplicity k, you might not have k linearly independent eigenvectors. You might have a smaller number of them. We'll look at an example of that later. But here, I'm saying, let's forget about that case, because it's ugly. Let's assume that whatever multiplicity each of these eigenvalues has, if you have an eigenvalue with multiplicity k, then you have k linearly independent right eigenvectors and k linearly independent left eigenvectors to correspond to that.
And then when you add up all of the eigenvectors, you have M linearly independent eigenvectors. And what happens when you have M linearly independent vectors in a space of dimension M? If you have M linearly independent vectors in a space of dimension N, you expand the whole space, which says that the vector of these eigenvectors is in fact non-singular, which says, again, we can do all of the stuff we did before. There's a little bit of a trick in showing that the left eigenvectors and the right eigenvectors can be made orthogonal. But aside from that, P to the n is again equal to the same form. And what this form says is, if all of the eigenvalues except one are less than 1, then you're again going to approach steady state.
What does that mean? Suppose I have more than one ergodic chain, more than one ergodic class, or suppose I have a periodic class or something else. Is it possible to have one eigenvalue equal to 1 and all the other eigenvalues be smaller? If there's one eigenvalue that's equal to 1, according to this formula here, eventually P to the n converges to that one value equal to 1. And right eigenvector can be taken as e.
Left eigenvector can be taken as a steady state vector pi. And we have the case of convergence. Can you have convergence to all the rows being the same if you have multiple ergodic classes? No. If you have multiple ergodic classes and you start out in one class, you stay there. You can't get out of it. If you have a periodic class and you start out in that periodic class, you can't have convergence there.
So in this situation here, where all the eigenvalues are distinct, you can only have one eigenvalue equal to 1. Here, when we're going to this more general case, we might have more than one eigenvalue equal to 1. But if in fact we only have one eigenvalue equal to 1, and all the others are strictly smaller in magnitude, then in fact you're just talking about this case of an ergodic unit chain again. It's the only place you can be.
So let's look at an example of this. Suppose you have a Markov chain which has l ergodic sets of states. You have one set of states. So we have one set of states over here, which will all go back and forth to each other. Then another set of states over here. Let's let l equal 2 in this case. So what happens in this situation? We'll have to work quickly before it gets up. Anybody with any sense, faced with a Markov chain like this, would say if we start here, we're going to stay here, if we start here, we're going to stay here. Let's just analyze this first. And then after we're done analyzing this, we'll analyze this. And then we'll put the whole thing together.
And what we will find is a transition matrix which looks like this. And if you're here, you stay here. If you're here, you stay here. We can find the eigenvalues and eigenvectors of this. We can find the eigenvalues and eigenvectors of this. If you look at this crazy formula for finding determinants, what you're stuck with is permutations within here times permutations within here. So the eigenvalues that you wind up with are products of the two eigenvalues. Or any eigenvalue here is an eigenvalue of the whole thing. Any eigenvalue here is an eigenvalue of the whole thing. And we just look at the sum of the number of eigenvalues here and the number there.
So we have a very boring case here. Each ergodic set has an eigenvalue equal to 1, has a right eigenvector equal to 1. When the steps of that state and 0 elsewhere. There's also a steady state vector on that set of states. We've already seen that. So P to the n converges to a block diagonal matrix, where for each ergodic set, the rows within that set are the same. So P to the n then is pi 1, pi 1. And then here, we have pi 2, pi 2, pi 2. So that's all that can happen here. This is limit.
So one message of this is that, after you understand ergodic unit chains, you understand almost everything. You still have to worry about periodic unit chains. But you just take a power of them, and then you have ergodic sets of states. one final thing. Good, I have five minutes to talk about this. I don't want any more time to talk about it, because I'll get terribly confused if I do. And it's a topic which, if you want to read more about it, read about it in Strang. He obviously doesn't like the topic either. Nobody likes the topic. Strang at least was driven to say something clear about it. Most people don't even bother to say something clear about it.
There's a theorem, due to, I guess, Jordan, because it's called a Jordan form. And what Jordan said is, in the nice cases we talked about, you have this decomposition of the transition matrix in P into a matrix here whose columns are the right eigenvectors times a matrix here, which is a diagonal matrix with the eigenvalues along it. And this, finally, is a matrix which is the inverse of this, and, which properly normalized, is the left eigenvectors of P. And you can replace this form by what's called a Jordan form, where P is equal to some matrix u times the Jordan form matrix j times the inverse of u.
Now, u is no longer the right eigenvectors. It can't be the right eigenvectors, because when we needed Jordan form, we don't have enough right eigenvectors to span the space. So it has to be something else. And like everyone else, we say, I don't care what that matrix is. Jordan proved that there is such a matrix, and that's all we want to know. The important thing is that this matrix j in here is as close as you can get it. It's a matrix, which along the main diagonal, has all the eigenvalues with their appropriate multiplicity. Namely, lambda 1 is an eigenvalue with multiplicity 2. Lambda 2 is an eigenvalue of multiplicity 3. And in this situation, you have two eigenvectors here, so nothing appears up there.
With this multiplicity 3 eigenvalue, there are only two linearly independent eigenvectors. And therefore Jordan says, why don't we stick a 1 in here and then solve everything else? And his theorem says, if you do that, it in fact works. So every time-- well, the eigenvalue is on the main diagonal, the ones on the next diagonal up, the only place would be anything non-0 is on the main diagonal in this form, and on the next diagonal up, where you occasionally have a 1. And the 1 is to replace the need for deficient eigenvectors. So every time you have a deficient eigenvector, you have some 1 appearing there. And then there's a way to solve for u. And I don't have any idea what it is, and I don't care.
But if you get interested in it, I think that's wonderful. But please don't tell me about it. Nice example of this is this matrix here. What happens if you try to take the determinant of P minus lambda i? Well, you have 1/2 minus lambda, 1/2 minus lambda, 1 minus lambda. What are all the permutations here that you can take? There's the permutation of the main diagonal itself. If I try to include that element, there's nothing I can do but have some element down here. And all these elements are 0. So those elements don't contribute to a determinant at all. So I have one eigenvalue which is equal to 1. I have two values at multiplicity 2, eigenvalue which is 1/2.
If you try to find the eigenvector here, you find there is only one. So in fact, this corresponds to a Jordan form, where you have 1/2. 1, and a 0, and a 1 here, and 0 everywhere else. And now if I want to find P to the n, I have u times this j times u to the minus 1 times u. All the u's in the middle cancel out, so I wind up eventually with u times j to the nth power times u to the minus 1.
What is j to the nth power? What happens if I multiply this matrix by itself n times? Well, it turns out that what happens is that this main diagonal here, you wind up with a 1/4 and then 1/8 and so forth. This term here, it goes down exponential. Well, if you multiply this by itself, eventually, you can see what's going on here more easily if you draw the Markov chain for it. You have state 1, state 2, and state 3. State 1, there's a transition 1/2 and a transition 1/2. State 2, there's a transition 1/2 and a transition 1/2, And state 3, you just stay there.
So the amount of time that it takes you to get to steady state is the amount of time it takes you-- you start in state 1. You've got to make this transition eventually, and then you've got to make this transition eventually. And the amount of time that it takes you to do that is the sum of the amount of time it takes you to go there, plus the amount of time that it takes to go there. So you have two random variables. One is the time to go here. The other is the time to go here. Both of those are geometrically decreasing random variables. When we convolve those things with each other, what we get is an extra term n. So we get an n times 1/2 to the n.
So the thing which is different in the Jordan form is, instead of having an eigenvalue to the nth power, you have an eigenvalue times-- if there's only a single one there, there's an n there. If there are two 1s both together, you get an n times n minus 1, and so forth. So worst case, you've got a polynomial to the nth power times an eigenvalue. For all practical purposes, this is still the eigenvalue going down exponentially. So for all practical purposes, what you wind up with is the second-largest eigenvalue still determines how fast you get convergence. Sorry, I took eight minutes talking about the Jordan form. I wanted to take five minutes talking about it. You can read more about it in the notes.
Markov Eigenvalues and Eigenvectors (PDF)


Harry Boyte (far right) receives the Spirit of Ghandi Award at the State Capitol. Photo courtesy of Chai News.
In recognition of his work promoting civic engagement and non-violence, Augsburg scholar Harry Boyte has been awarded the Spirit of Gandhi Award. Boyte, the senior scholar in public work philosophy with the Augsburg Sabo Center for Democracy and Citizenship, describes the spiritual, moral, and psychological aspects of non-violence in his Huffington Post article on the new nonviolence movement. The Spirit of Gandhi Award is given in celebration of Nonviolence Day, a globally-observed day on the birth anniversary of Mahatma Gandhi, at the Minnesota State Capitol.
Read more on the Spirit of Gandhi Award here.


Residents were furious when Lossiemouth Harbour Board erected bollards at each end of a lane which links the west and east ends of the port earlier this year.
The board insisted the bollards were necessary for health and safety reasons.
But local businessmen and community groups protested that blocking off the busy access road would have a dire effect the area's tourist trade.
A public meeting was later held, and various other safety measures were put forward, which Moray Council officers agreed to investigate.
But now, more than two months on and with the tourist season set to begin, campaigners are growing impatient.
John Thomson, who owns the Shore Street petrol station, has spearheaded the crusade to have the bollards removed.
He said: "Things seem to have stalled. I like to get a job done as quickly as possible and I fear the longer the bollards are there the less chance we have of getting them taken down.
"There is no question that over the summertime we will lose passing and holiday trade.
Mike Mulholland, of Lossiemouth Community Council, added: "We're still waiting for Moray Council officers to suggest ways we can make it safer that we can pass onto the harbour board.
"When we had the meeting it sounded encouraging, but we really don't want to wait much longer as the public are all desperate to know what will happen.
"I'd like to see them at least be removed for the summer, but I've a feeling we might not get that.
A harbour board spokesman yesterday said that many pedestrians had praised the recent changes, and stressed that safety would be the main priority in any decisions it makes.
Heldon and Laich councillor John Cowe has taken on residents' concerns with the local authority, but was unavailable for comment yesterday.


Melvin N. Puu, a Unit 14 member and water safety officer with the City & County of Honolulu Emergency Medical Services Department, Ocean Safety Division, has been elected to a two-year term as Unit 14 state board director.
The special election for a Unit 14 director was conducted due to the establishment of the new Bargaining Unit 14 to represent state law enforcement officers and state and county ocean and water safety officers.
The deadline for protest or challenges on the conduct of this special Unit 14 election is March 6. The State Board of Directors will certify the results at its March 29 board meeting.
Ballots were mailed to eligible HGEA members on January 28, and counted in the Honolulu headquarters building today, March 1, by N&K CPAs, Inc., HGEA's third-party elections organization.


My daughter who is three years old just got accepted as a model for Images of St. Louis.
We had to pay a fee of $300 for headshots, since she has no professional pictures. We signed a contract saying that was the only fee to come out of our pockets and that they would get 20% of what she makes on jobs only if they book it for her.
They said we were free to sign a contract with any modeling agency of our choice at any time, since they weren't guaranteeing us job placement, just auditions.
I was wondering if this company was a scam, and if so, do you have any other info on this?
Professional pictures at that price or professional pictures at any price for that age is questionable. Most people advise against professional pictures for children four years old and younger (their look changes so fast). That price is about average for adults. Non-professional snapshots should do fine for babies and small children, or much lower rates, e.g. Sears Portrait Studio.
In any case, look at the BBB record for Images of St. Louis. It shows Images of St Louis is associated with the John Robert Powers School. JRP has a history of getting students to pay for professional photos which are overpriced or aren't necessary.
It is widely established that potential models should not pay for professional pictures until they get signed with an agency. JRP also charges for classes which are are overpriced or aren't necessary. Since the owner of Images is also an owner of a JRP franchise, it is certainly possible you have paid for pictures which were overpriced or aren't necessary.
The BBB report for Images, as shown below, reveals "consumers have complained that few, if any, jobs have resulted from enrolling in the classes or the photographs."
file was opened or over the last 36 months, whichever is less.
The Bureau processed 16 complaints about this company in the last 36 months.
closed as Resolved, and 1 was closed as Disputed.


Due to circumstances beyond our control, our usual monthly date at the Old White Swan in York has moved for October and November. We will now be playing at the Swan on Friday 21 October, and Friday 18 November. Everything goes back to normal from December.
If you can't wait that long to see us, don't panic! We are playing at O'Neill's in York on Friday 7 October, and Friday 4 November. We'll see you there.


Surrounded by wildfires, an Oklahoma meteorologist was looking over his shoulder for a way out when he spotted a heifer and calf wandering in the direct line of a massive fire.
"Look at this, oh my gosh," said Mike Armstrong, swinging his camera around. "There is a baby calf right over here. Oh. My. Gosh."
Armstrong, a field meteorologist for KOCO-TV in Oklahoma City, was reporting on the Rhea wildfires near Canton, Okla., on Tuesday when he spotted the cattle. It was midafternoon, between broadcasts, so he was sharing footage on Facebook Live with chief meteorologist Damon Lane in the studio.
"Do you need to cut that fence?" Lane asked. "If you need to go, go. Oh my gosh."
Armstrong was on one side of the fence, the cow and calf were on the other, and beyond them was the quickly approaching fire line.
"Does that not drive it home how serious this is?" Armstrong said. "There's nothing I can do."
It wasn't the only thing on his mind. He'd just seen a propane tank explode only 50 feet away, and was watching homes erupt in flames nearby. Armstrong had to consider his own safety, as well as his broadcast. Not to mention the firenadoes, which are exactly as terrifying as they sound.
The cattle, however, took his attention.
He backed his truck down the road and found an unlocked gate. When he pulled into the pasture, he spotted a larger cattle herd on a hillside even closer to the fire line.
Armstrong grew up around cattle, but wasn't a rancher. After he got into the pasture, he wasn't sure how to coax the herd toward safety.
"That hill's on fire, and there's no way to get to them all," he said. "Damon, I don't know if I can do this. I can try."
Luckily, his Facebook audience came to the rescue. The commenters watching the rescue effort live all had the same advice.
"Honk your horn they will come you're way"
It's a rancher's trick to call cattle, Armstrong said. They'll think it is feeding time and come running. Sure enough, the herd followed him off the fiery hillside and toward green grass.
"Good job, Michael, everyone on Facebook is throwing up a bunch of heart signs," Lane said in the studio. "Michael, you have a lot of people who are crying now on Facebook. I'm getting a lot of comments."
"Now don't say that, you're going to make me emotional," Armstrong said, showing the cattle following his truck toward safety.
After taking the larger herd to safety, he doubled back to find the calf and cow that first caught his eye. When he used the honking technique the second time, it scared the calf. As the calf ran toward the fire line, Armstrong scrambled to open fences and follow.
Just when it seemed that the calf was going straight for the flames, it darted behind the fire line, "in the black," Armstrong said.
"As long as he stayed in that area, he'd be OK," Armstrong said. "At least it got them out of that field."
When he was reasonably sure the cattle were out of harm's way, he left the pasture, making sure to keep the gates open just in case the cattle needed to escape. He said he was surprised by how quickly the story spread.
"I've been in media now about 20 years," he said. "I've been a part of big stories, but cattle is No. 1 right now."
As for the herd, Armstrong said another KOCO reporter went to that pasture today and reports that the calf and heifer are still there, as is a large herd nearby.


text:
Let's have a political police agency
by Steve Harris on November 12, 2008
As Plato pointed out, when a civilization gets too disorganized — through pluralism, a result of democracy — it falls prey to tyrants who appeal with a lowest common denominator message.
Broun was specifically referring to a July speech by Obama, where the then-Democratic presidential nominee said he supports a civilian force helping the military when it comes to national security: "We've got to have a civilian national security force that's just as powerful, just as strong, just as well funded [as the military]," Obama said in the speech that was largely a call to national service.
Responding to those comments, Broun told the AP Monday: "That's exactly what Hitler did in Nazi Germany and it's exactly what the Soviet Union did. When he's proposing to have a national security force that's answering to him, that is as strong as the U.S. military, he's showing me signs of being Marxist."
"We can't be lulled into complacency," Broun added. "You have to remember that Adolf Hitler was elected in a democratic Germany. I'm not comparing him to Adolf Hitler. What I'm saying is there is the potential."
Those who run on the most ambiguous platforms have the greatest potential to be tyrants.
Then again, history created both men: fear of Bolshevism created Hitler, and fear of instability in Russia created Stalin.
Stupid strategies beget stupid results | Ron Paul: Republicans need to rediscover roots

redpajama_set_name: RedPajamaCommonCrawl


The president urged ministers attending the meeting to conduct a thorough and careful study of the draft law to ensure that it did not run counter to the state`s Constitution.
"This is the House of Representatives` initiative. We need to study it thoroughly and carefully," he noted in his opening remarks at the meeting, attended by Coordinating Economic Minister Darmin Nasution and several other cabinet members.
President Widodo remarked that oil and gas had become the main and strategic source of financing for Indonesia`s development, so the oil and gas law should be able to strengthen the country`s energy security.
The oil and gas draft law should also be able to help secure the country`s national production and strengthen capacities of local industries as well as boost the development of skilled human resources in the oil and gas sector, he explained.
President Widodo added that the draft law initiated by the House should be able to ensure that the country`s oil and gas management could be simpler as well as more efficient and transparent.
The limited cabinet meeting was also attended by Coordinating Minister for Political, Legal, and Security Affairs Wiranto, Finance Minister Sri Mulyani, Minister of State-Owned Enterprises Rini Sumarno, and Maritime Affairs and Fisheries Minister Susi Pudjiastuti.
Coordinating Human Development and Culture Minister Puan Maharani, Law and Human Rights Minister Yasona H. Laoly, Transportation Minister Budi Karya Sumadi, Environment and Forestry Minister Siti Nurbaya, and State Secretary Minister Pratikno also attended the meeting.
In the past, Indonesia had also become one of the members of the Organization of Petroleum Exporting Countries. However, it has currently become a net importer of oil owing to shrinking domestic production and fast-growing domestic consumption.
To deal with this condition, state-owned energy company PT Pertamina is expected to use all crude oil production of Indonesia to be fed into its refineries in 2019.
Pertamina has to spend large amount of foreign exchange for the import of crude oil to feed its oil refineries, as most of the country`s crude oil production by oil contractors is exported.


Macro environment includes a variety of political, social, technological, economic, and demographical components, which have an external impact on businesses. Thus, the main aim of this paper is to explore an article about how rivals factor of the macro environment affects BMW, Toyota Motor Corporation and the whole car manufacturing industry.
The chosen article focuses on rival factor, which is one of the core macro-environmental components. Few days ago, BMW and Toyota Motor Corporation decided to form an alliance to challenge global rivals due to extremely tough competition in the industry (Hetzner, 2012).
The news in this article is relevant to both companies, BMW and Toyota, since they are involved in cooperation with their direct competitors. The agreement might result in successful partnership between the two leading companies. The cars they are going to produce together will have the benefit of reduced weight thanks to BMW. Additionally, Toyota Motor Corporation can offer strong innovations base it has in manufacturing electric cars (Hetzner, 2012). As a result, the car producers can introduce more efficient cars to the market and increase their profits.
Moreover, the collaboration between BMW and Toyota Motor Corporation can be an external threat for other car manufacturers. The companies can improve their positions in the market and outrival other competitors thanks to the deal.
Car manufacturing sector is highly developed and influential. Changes in the segment do not only affect the companies, which participate in it. In fact, they have a large impact on the national economies, international trade, and national income indexes. Hence, the news is essential for the governments of Japan, Germany and other countries, which economies strongly depend on car manufacturing industry. There is a clear possibility of some countries' national income decrease caused by intensification of German and Japanese car producing.
Therefore, the information in this article is of a particular interest for car producers and national governments since it can result in market share and income changes.


text:
Home Newsroom News and press releases Latest from the OSCE Special Monitoring Mission to Ukraine (SMM), based on information received as of 19:30, 10 March 2019
OSCE Special Monitoring Mission to Ukraine
Latest from the OSCE Special Monitoring Mission to Ukraine (SMM), based on information received as of 19:30, 10 March 2019
KYIV 11 March 2019
Media accreditation to OSCE events
OSCE Live
This report is for the media and the general public.
Compared with the previous reporting period, between the evenings of 8 and 9 March, the SMM recorded more ceasefire violations in Donetsk region and significantly fewer in Luhansk region.
Compared with the previous 24 hours, between the evenings of 9 and 10 March, the Mission recorded fewer ceasefire violations in Donetsk and Luhansk regions.
The SMM observed fresh shrapnel damage in Zolote-5/Mykhailivka and fresh impact craters near Popasna.
The Mission recorded ceasefire violations inside the Petrivske disengagement area. It observed for the first time anti-tank mines on a road in the southern part of the Zolote disengagement area.
The SMM observed weapons in violation of the withdrawal lines on both sides of the contact line.
The Mission monitored adherence to the ceasefire to enable repairs to essential civilian infrastructure in Luhansk region.
Restrictions of the SMM's access continued in all three disengagement areas and elsewhere. Its access was also restricted in non-government-controlled areas, including at checkpoints near Novolaspa and Novoazovsk and at a heavy weapons holding area in Donetsk region.*
In Kyiv, the Mission monitored a gathering and later saw a car partially covered with soot.
Ceasefire violations[1]
In Donetsk region, between the evenings of 8 and 9 March, the SMM recorded more ceasefire violations, including, however, fewer explosions (about 120), compared with the previous reporting period (about 185 explosions). The majority of the ceasefire violations were recorded at southerly and south-westerly directions of Zaichenko (non-government-controlled, 26km north-east of Mariupol), at south-easterly and southerly directions of Chermalyk (government-controlled, 31km north-east of Mariupol) and at south-westerly directions of Svitlodarsk (government-controlled, 57km north-east of Donetsk). During the day on 9 March, positioned in Pikuzy (formerly Kominternove, non-government-controlled, 23km north-east of Mariupol), the SMM heard three explosions assessed as outgoing rounds and about 40 bursts and shots of small-arms fire, all at an assessed distance of 1km south of the nearest civilian house in the village (for recent observations of damage caused by shelling in Pikuzy, see SMM Daily Report 4 March 2019).
Between the evenings of 9 and 10 March, the Mission recorded fewer ceasefire violations, including 50 explosions, compared with the previous 24 hours. The majority of ceasefire violations were recorded at south-westerly and north-westerly directions of Yasynuvata (non-government-controlled, 16km north-east of Donetsk).
In Luhansk region, between the evenings of 8 and 9 March, the SMM recorded significantly fewer ceasefire violations (about 70 explosions), compared with the previous reporting period (about 930 explosions). The majority of the ceasefire violations were recorded on the morning of 9 March in an area between Popasna (government-controlled, 69km west of Luhansk) and Molodizhne (non-government-controlled, 63km west of Luhansk), near the forward positions of the Ukrainian Armed Forces and the armed formations (see below Fresh impact craters observed near Popasna).
Between the evenings of 9 and 10 March, the Mission recorded fewer ceasefire violations, including 11 explosions, compared with the previous 24 hours. The majority of ceasefire violations were recorded in areas north-west of Berezivske (non-government-controlled, 53km north-west of Luhansk).
Fresh shrapnel damage in Zolote-5/Mykhailivka
At 34 Poshtova Street in Zolote-5/Mykhailivka (non-government-controlled, 58km west of Luhansk), on 9 March, the SMM saw shrapnel damage to the north-north-east-facing side of a one-storey inhabited house and to an adjacent north-north-east-facing metal fence. About 4-5m north of the house, the Mission observed a fresh impact crater. The SMM was not able to assess the type of weapon used or the direction of fire. A resident of the house told the SMM that at around 19:00 on 7 March, while at home with his wife, he had heard a loud explosion nearby.
Fresh impact craters observed near Popasna
On 9 March, the SMM observed five to six fresh impact craters, assessed as caused by 82mm mortar rounds, in a field about 50m from the forward position of the Ukrainian Armed Forces in the eastern outskirts of Popasna (about 1.5km west of the nearest civilian house). The SMM heard over 60 undetermined explosions in the area earlier that day (see above).
Disengagement areas[2]
On 7 March, inside the disengagement area near Zolote (government-controlled, 60km west of Luhansk), an SMM mini-unmanned aerial vehicle (UAV) again spotted an infantry fighting vehicle (IFV) (BMP variant) near the road bridge, assessed as belonging to the Ukrainian Armed Forces.
On 9 March, the SMM observed for the first time six anti-tank mines (probable TM-62) laid out in a ditch extending across road T-1316, about 60m north of the checkpoint of the armed formations on the southern edge of the Zolote disengagement area (for similar observations in the area, see SMM Daily Report 28 February 2019).
During the day on 9 March, positioned about 1.5km north of Petrivske (non-government-controlled, 41km south of Donetsk), the SMM heard and saw two explosions assessed as impacts at an assessed range of 2-3km south-south-west (assessed as inside the disengagement area).[3] Shortly thereafter, the Mission observed two military vehicles – one IFV (BMP variant) and the other undetermined – near a previously observed position of the Ukrainian Armed Forces inside the disengagement area, about 550m south of the area's northern edge and about 1.3km west of its eastern edge.
The same day, positioned inside the disengagement area near Stanytsia Luhanska (government-controlled, 16km north-east of Luhansk), the SMM observed a calm situation.
On 10 March, positioned inside the Stanytsia Luhanska disengagement area and near the Zolote and Petrivske disengagement areas, the SMM observed calm situations.
Withdrawal of weapons
The SMM continued to monitor the withdrawal of weapons in implementation of the Memorandum and the Package of Measures and its Addendum.
Weapons in violation
Government-controlled areas
The SMM saw:
a surface-to-air missile system (9K33 Osa) moving north in Klynove (68km north-east of Donetsk) and
nine self-propelled howitzers (2S1 Gvozdika, 122mm) near Novoaidar (49km north-west of Luhansk).
Non-government-controlled areas
a surface-to-air missile system (9K35 Strela-10) east of Zaichenko with multiple fresh marks from tracked vehicles (types mostly undetermined) nearby, including those assessed as tracks of an MT-LB-based chassis, running in both directions along a road leading to Roza (36km north-east of Mariupol).
Beyond withdrawal lines but outside of designated storage sites
three self-propelled mortars (2S9 Nona-S, 120mm) near Novoaidar.
Weapons that the SMM could not verify as withdrawn[4]
At heavy weapons holding areas in government-controlled areas of Luhansk region
The SMM noted that:
four anti-tank guns (MT-12 Rapira, 100mm) were present for the first time and
31 self-propelled howitzers (28 2S3 Akatsiya, 152mm and three 2S1), 22 anti-tank guns (MT-12) and two self-propelled mortars (2S9) were again missing.
four multiple launch rocket systems (MLRS) (BM-21 Grad, 122mm) were present for the first time and
two MLRS (BM-21) were missing for the first time.
At heavy weapons holding areas beyond the respective withdrawal lines in non-government-controlled areas of Donetsk region
13 MLRS (BM-21) and eight anti-tank guns (MT-12) were again missing.
Indications of military and military-type presence inside the security zone[5]
An SMM mini-UAV spotted:
three armoured combat vehicles (ACV) (type undetermined) near Popasna.
An SMM mid-range UAV spotted:
an ACV (type undetermined) near Lebedynske (16km east of Mariupol).
- two IFVs (BMD variants) in the yard of a former school building in Voitove (33km north-west of Luhansk).
an armoured personnel carrier (APC) (BTR-70) near Berezove (31km south-west of Donetsk) and
two armoured reconnaissance vehicles (BRDM-2) and five IFVs (three BMP-2 and two undetermined variants) near Muratove (51km north-west of Luhansk).
a military truck (KamAZ-4310) moving south on a road on the western edge of Sakhanka (24km north-east of Mariupol), in a zone within which deployment of heavy armaments and military equipment is proscribed according to Point 5 of the Memorandum of 19 September 2014.
an IFV (BMP-1) and a probable APC (MT-LB) near Vesela Hora (16km north of Luhansk).
two IFVs (BMP-2) near Debaltseve (58km north-east of Donetsk) and
two IFVs (BMP-1 and BMP-2) and one APC (BTR-80) near Holubivka (formerly Kirovsk, 51km west of Luhansk).
Presence of mines
On 8 March, an SMM mini-UAV spotted for the first time at least 40 anti-tank mines (probable TM-62) in a field about 50m south-west of road H-21, west of Vesela Hora.
SMM facilitation of repairs to civilian infrastructure
On 9 March, the SMM monitored adherence to the ceasefire to enable repairs to damaged power lines near Zolote-5/Mykhailivka. On 9 and 10 March, the Mission continued to facilitate the operation of the Donetsk Filtration Station (15km north of Donetsk) and to monitor the security situation around the pumping station near Vasylivka (non-government-controlled, 20km north of Donetsk).
Border areas outside government control
On 9 March, while at a border crossing point near Uspenka (73km south-east of Donetsk) for about one hour, the SMM saw seven cars (three with Russian Federation and one with Ukrainian licence plates, and three with "DPR" plates) entering Ukraine and 25 cars (13 with Russian Federation, five with Ukrainian and one with Georgian licence plates, and six "DPR" plates), 13 covered cargo trucks (six with Ukrainian licence plates, and seven with "DPR" plates) and a bus (with "DPR" plates) exiting Ukraine.
The following day, while at the same border crossing point for 30 minutes, the SMM saw 13 cars (one with Ukrainian and one with Russian Federation licence plates, and 11 with "DPR" plates) entering Ukraine and 13 cars (eight with Russian Federation and one with Ukrainian licence plates, and four "DPR" plates) and a bus (with Russian Federation plates) exiting Ukraine.
On 9 March, while at a pedestrian border crossing point near Ulianivske (61km south-east of Donetsk) for 20 minutes, the SMM saw four pedestrians (two men and two women, aged 20-30) entering Ukraine and two pedestrians (a man and a woman in their sixties) exiting Ukraine.
On 10 March, while at a pedestrian border crossing point near Novoborovytsi (79km south of Luhansk) for 30 minutes, the SMM observed no cross-border traffic.
Gathering in Kyiv
In Kyiv, on 9 March, the SMM monitored a pre-announced gathering by National Corps at Independence Square. The Mission saw about 600-700 people (mostly young men), some of whom were holding flags of National Corps, and heard them expressing opposition to the government. Later in the day, the SMM followed up on reports that smoke grenades had been thrown in front of the presidential administration building. The SMM saw that the side of a car parked in front of the building was partially covered with soot, assessed as caused by a smoke grenade. The Mission observed about 40-50 police officers present in the area. It did not observe further incidents.
The Mission continued monitoring in Kherson, Odessa, Lviv, Ivano-Frankivsk, Kharkiv, Dnipro and Chernivtsi.
*Restrictions of SMM's freedom of movement or other impediments to fulfilment of its mandate
The SMM's monitoring and freedom of movement are restricted by security hazards and threats, including risks posed by mines, unexploded ordnance (UXO) and other impediments – which vary from day to day. The SMM's mandate provides for safe and secure access throughout Ukraine. All signatories of the Package of Measures have agreed on the need for this safe and secure access, that restriction of the SMM's freedom of movement constitutes a violation, and on the need for rapid response to these violations. They have also agreed that the Joint Centre for Control and Co-ordination should contribute to such response and co-ordinate mine clearance. Nonetheless, the armed formations in parts of Donetsk and Luhansk regions frequently deny the SMM access to areas adjacent to Ukraine's border outside control of the Government (for example, see SMM Daily Report of 8 March 2019). The SMM's operations in Donetsk and Luhansk regions remain restricted following the fatal incident of 23 April 2017 near Pryshyb; these restrictions continued to limit the Mission's observations.
Denial of access:
On 9 March, an armed member of the armed formations at a checkpoint on the eastern edge of Novolaspa (non-government-controlled, 50km south of Donetsk) denied the SMM access to the village.
On 10 March, the SMM was unable to visit a heavy weapons holding area in a non-government-controlled area of Donetsk region. The site was locked and no personnel were present.
At a checkpoint north of Novoazovsk (non-government-controlled, 40km east of Mariupol), two armed members of the armed formations again denied the SMM access to the town, citing "ongoing special police operation".
Regular restrictions related to disengagement areas and mines/UXO:
The sides continued to deny the SMM full access to the three disengagement areas, as well as the ability to travel certain roads previously identified as important for effective monitoring by the Mission and for civilians' movement, through failure to conduct comprehensive clearance of mines and UXO.
Conditional access:
On 9 March, at a checkpoint near Staromykhailivka (non-government-controlled, 15km west of Donetsk), a member of the armed formations again allowed the SMM to pass only upon being escorted by a car of the armed formations during its patrol of the settlement.
Delay:
On 9 March, at a checkpoint 600m west of Verkhnoshyrokivske (formerly Oktiabr, non-government-controlled, 29km north-east of Mariupol), a member of the armed formations prevented the SMM from proceeding to Pikuzy (formerly Kominternove, non-government-controlled, 23km north-east of Mariupol) and Sakhanka (non-government-controlled, 24km north-east of Mariupol). After about 40 minutes, he allowed the SMM to pass through the checkpoint.
[1] Please see the annexed table for a complete breakdown of the ceasefire violations as well as a map of the Donetsk and Luhansk regions marked with locations featured in this report. During the reporting period, the SMM camera in Krasnohorivka was not operational.
[2]Disengagement is foreseen in the Framework Decision of the Trilateral Contact Group relating to disengagement of forces and hardware of 21 September 2016.
[3] Due to the presence of mines, including on a road between Bohdanivka and Petrivske, the SMM cannot access its camera in Petrivske, and thus the SMM has not been able to access observations from the camera since 22 June 2018.
[4] The SMM observed weapons that could not be verified as withdrawn, as their storage did not comply with the criteria set out in the 16 October 2015 notification from the SMM to the signatories of the Package of Measures on effective monitoring and verification of the withdrawal of heavy weapons. The SMM noted that two such sites continued to be abandoned.
[5] The hardware mentioned in this section is not proscribed by the provisions of the Minsk agreements on the withdrawal of weapons, unless otherwise specified.
Liudmyla Palamar
National Outreach Officer
26 Turhenievska Street
01054 Kyiv
Office: +380 44 392 0965
liudmyla.palamar@osce.org
smm-media@osce.org
Table of ceasefire violations as of 10 March 2019
Daily updates from the Special Monitoring Mission to Ukraine
Thematic Reports from the Special Monitoring Mission to Ukraine
OSCE Special Monitoring Mission to Ukraine -- The Facts
Sign up to receive OSCE Press Releases and Media Advisories by email as soon as they have been published.
The OSCE activity described on this page is helping to achieve the following Sustainable Development Goals


We are a leading Manufacturer of z shape modular bucket elevator from Pune, India.
We offer a precision engineered range of Modular Bucket Belt Conveyor that are incorporated with the features of heavy loads accumulation and drainage or air flow. Constructed from best quality components, these belt conveyors feature sturdy construction, optimum performance and long working life. Further, we also customize these belt conveyors as per the specifications laid down by our clients.
Looking for Z Bucket Elevator ?


text:
Ad Network Performs In-Browser Cryptojacking
An ad network provider is performing in-browser Coinhive cryptojacking on websites that use its service, 360 Netlab security researchers warn.
The practice has been ongoing since December 2017, several months after the ad network provider, a company called PopAds Publisher, started using domain generation algorithm (DGA) technology to bypass ad blockers, claiming it would allow customers to "monetize traffic that wasn't monetized before."
In mid-2017, the provider started to generate seemingly random domains that would ensure ads can reach end users. By the end of the year, however, these domains, which 360 Netlab refers to as DGA.popad, started participating in cryptojacking activities, all without end-users' acknowledgement.
Given that many people use ad blockers to prevent sites from displaying ads to them, ad networks often attempt to bypass blockers, and this provider decided to use DGA domains to host its advertisements. With these domains changing daily, it becomes difficult to block the ads, the researchers point out.
What's more, the ad network provider recently started using the DGA.popad domains to perform cryptojacking. These domains, the researchers discovered, have a strong connection with Coinhive family domains in DNS traffic and serve the coinhive.min.js web miner.
Some of the DGA.popad domains have a high ranking, with one of them found in the top 2000 sites on Alexa and several others in the top 3000 list.
Once a user accesses such a site, their computer's CPU starts being used to the full. According to 360 Netlab, the favicon.ico on the DGA.popad sites was found to run as a web miner. Most of the sites that would redirect users to DGA.popad domains are providing adult content and downloading services.
Because the impacted sites contain advertisements from this ad network, the cryptojacking activities are performed regardless of whether the user has an ad blocker installed or not.
Normally, users would be sent to a standard domain (serve.popads.net) hosted by the provider. If an ad blocker is used, the standard domain is blocked, and the visitor is sent to one of the DGA.popad domains. Regardless of the delivery mechanism, both the ad and the cryptojacking script are served.
"We are not able to make a detailed assessment. This is because only part of all the traffic passing DGA.popad and serve.popads.net will be inserted a web miner, but we are not sure which part will be selected, for now," the researchers note.
Related: Crypto-Mining Attack Targets Web Servers Globally
Related: Oracle WebLogic Server Flaw Exploited to Deliver Crypto-Miners
CISA, NSA, and MS-ISAC issued an alert on the malicious use of RMM software to steal money from bank accounts.
Chinese Hackers Adopting Open Source 'SparkRAT' Tool
Chinese threat actor DragonSpark has been using the SparkRAT open source backdoor in attacks targeting East Asian organizations.
Learning to Lie: AI Tools Adept at Creating Disinformation
Artificial intelligence is competing in another endeavor once limited to humans — creating propaganda and disinformation.


The Mining Hub is committed to protecting the privacy of personal information collected from registered users of our website.
When you submit your personal details, work experience, educational qualifications, and other personal information to an advertiser for a specific job application, that information will only go to the advertiser you have submitted to. The information you have submitted will only be retained and used by that employer for the purposes of contacting you regarding possible employment opportunities.
You may review, correct or delete your personal information at any time. Simply log onto your Account and make the necessary changes. If you change your mind about receiving newsletters and alerts, you may change your preferences by logging into your Account.
Our site links to other sites. When you click on one of these links, you are leaving our site and entering another site that we do not control and for which we cannot be responsible. You should carefully review the privacy statements or policies on any other site that you visit because those privacy statements or policies will apply to your visit to that site and may be different from our policy.
For Company Subscribers, the primary purpose for which we collect information from you is to provide you with the products and services you have requested. If you do not register as a Company Subscriber and complete all the fields requested, you will not be able to access features and services that are reserved for these subscribers only.
When you submit your bank account or credit card details to facilitate payment for a service offered on this site, appropriate security measures are in place to ensure the confidentiality, security and privacy of that information. Bank account or credit card details used to facilitate payment for services are not retained in any way on the website.
When purchasing from The Mining Hub, your financial details are passed through a secure server using the latest 128-bit SSL (secure pockets layer) encryption technology. 128-bit SSL encryption is approximated to take at least one trillion years to break and is the industry standard.
Unfortunately, no data transmission over the internet can be guaranteed as totally secure. Whilst we strive to protect such information, we do not warrant and cannot ensure the security of any information which you transmit to us. However, every reasonable step to preserve the security of all and any information you lodge with The Mining Hub is taken.
We may, from time to time, review and update this Policy, including taking into account new or amended laws, new technology and or changes to our operation. All personal information held by us will be governed by the most recent updated policy.
Contact us at info@themininghub.com if you have any questions or concerns about the privacy of your details or about our privacy practices.


So here they are to scale (except for the probe), and with 4 Vesta's current position at 23 million km compared to the previous distance at 28 million km near the beginning of September, shown in grey.
Arrival is now in 297 days, but since Dawn uses an ion engine this arrival will just involve the tiniest bit of energy used for orbital insertion, and long before that the probe will be nearby and will have almost matched orbits with Vesta. The chief engineer updates Dawn's progress at the end of each month so we'll see if he has anything to say on the subject (from when we'll be able to see 4 Vesta from Dawn) this time. I emailed them on it a few weeks ago on that as well.


text:
Q: Data Wrangling in r I'm trying to subset my data into new data frames in order to perform my analysis.
I have data frames on 134 samples which contain a lot of information however I'm only interested in the type, name, and expression columns.
How can I make a loop (or some way of iterating) through all of these samples to subset them into my columns of interest and then collate these into a single data frame?
I specifically would like to subset the data by type into separate tables (i,e one for proteins and one for mRNA). Then I would like to include the protein (or mRNA) expression of all of the samples into a single data frame.
Any help or direction to further resources would be much appreciated!
Thanks
Example code
sample1 <- data.frame(type = c("protein", "mRNA"),
                      name = c("DIABLO", "X1345"),
                      expression = c("1.23", "4.265"))

sample2 <- data.frame(type = c("protein", "mRNA"),
                      name = c("DIABLO", "X1345"),
                      expression = c("3.24", "5.33"))

sample3 <- data.frame(type = c("protein", "mRNA"),
                      name = c("DIABLO", "X1345"),
                      expression = c("2.56", "8.11"))

A: Combine all the data together in one dataframe and split on type column to get two separate dataframes.
library(dplyr)

bind_rows(mget(paste0('sample', 1:3)), .id = 'sample') %>%
  split(.$type) %>%
  list2env(.GlobalEnv)

mRNA

#   sample type  name expression
#2 sample1 mRNA X1345      4.265
#4 sample2 mRNA X1345       5.33
#6 sample3 mRNA X1345       8.11

protein
#   sample    type   name expression
#1 sample1 protein DIABLO       1.23
#3 sample2 protein DIABLO       3.24
#5 sample3 protein DIABLO       2.56

Although you might not need it but I added an extra column (sample) to identify from which dataframe each row is coming from.

A: First read all of the data files into a list and then extract the parts:
samples <- list(sample1, sample2, sample3)
# Alternate approach
# samples <- lapply(paste0("sample", 1:3), get)
proteins <- do.call(rbind, lapply(samples, function(x) x[1, ]))
proteins
#      type   name expression
# 1 protein DIABLO       1.23
# 2 protein DIABLO       3.24
# 3 protein DIABLO       2.56
mRNAs <- do.call(rbind, lapply(samples, function(x) x[2, ]))
mRNAs
#    type  name expression
# 2  mRNA X1345      4.265
# 21 mRNA X1345       5.33
# 22 mRNA X1345       8.11

redpajama_set_name: RedPajamaStackExchange


much more you can browse the search on this web site. When we discuss Does 2018 Toyota CH R offer Toyota Safety Sense TSS P after that we will certainly think of 2018 Toyota and also lots of things. But occasionally we should find out about to know much better. It is not far away with the extremely important . If you wish to open up the image gallery please click photo picture below. You can also download for your photo collection.