Dear Eric and William,
Why do the 1st and 2nd incantation of arima return sigma^2 as 5.233 vs
.?
The help for arima says ---> sigma2: the MLE of the innovations variance.
By that account the 1st result is incorrect. I am a little confused.
set.seed(123)
b <- arima.sim(list(order =
helloplease i have this matrixx<-rnorm(6,0,1)
aa<-matrix(x,nrow=6,ncol=6)
i have to extrat non diagonal value, i use this code
matrix(as.numeric(aa)[!as.numeric(aa) %in% diag(aa)],5,6)
but i didn't get the resultthank you
[[alternative HTML version deleted]]
helloi didn't obtaine the matrix after extrat non diagonalmy
programx<-rnorm(6,0,1)
aa<-matrix(x,nrow=6,ncol=6)
matrix(as.numeric(aa)[!as.numeric(aa) %in% diag(aa)],5,6)
nrow=5ncol=6thank you
[[alternative HTML version deleted]]
__
# construct the dataframe
`TK-QUADRANT` <- c(9161,9162,9163,9164,10152,10154,10161,10163)
LAT <-
c(55.07496,55.07496,55.02495,55.02496,54.97496,54.92495,54.97496,54.92496)
LON <-
c(8.37477,8.458109,8.37477,8.45811,8.291435,8.291437,8.374774,8.374774)
df <-
Hi Ashim,
Per the help page for arima(), it fits an ARIMA model to the specified time
series - but the caller has to specify the order - i.e. (p,d,q) - of the
model.
The default order is (0,0,0) (per the help page). Hence your two calls are
different. The first call is equivalent to order=c(0,0,0)
When that arrived it was a complete mess since you posted in HTML which
scrambles your code and you sent code which had syntax errors. Please
try again by posting in plain text and cut and paste your code. It would
also help if you stated exactly what you expected your output to consist of.
Hi.
You did not specify what do you want to do with the result.
functions
upper.tri, lower.tri and diag can manipulate parts of matrices.
Cheers
Petr
> -Original Message-
> From: R-help On Behalf Of malika yassa via R-
> help
> Sent: Wednesday, November 14, 2018 12:20 PM
> To: R-help
Hi
Your mail is mess due to HTML formating. Please use plain taxt mail.
You got an advice, did you try it?
With your code you just remove diagonal elements from your matrix. If this is
not your intention, you should specify more clearly what do you want to achieve
as the result.
Cheers
Petr
>
Another way:
> A <- matrix(1:9,3,3,
dimnames=list(Row=paste0("r",1:3),Col=paste0("c",1:3)))
> A
Col
Row c1 c2 c3
r1 1 4 7
r2 2 5 8
r3 3 6 9
> matrix( A[row(A)!=col(A)], nrow(A)-1, ncol(A), dimnames=list(NULL,
colnames(A)))
c1 c2 c3
[1,] 2 4 7
[2,] 3 6 8
Bill
An even better solution because it has fewer steps.
A <- matrix(1:9, 3, 3)
A
B <- A[-1, ]
B[upper.tri(B, diag=FALSE)] <- A[upper.tri(A)]
B
> A <- matrix(1:9, 3, 3)
> A
[,1] [,2] [,3]
[1,]147
[2,]258
[3,]369
> B <- A[-1, ]
> B[upper.tri(B, diag=FALSE)] <-
Dear Eric,
Many thanks for your reply.
Best Regards,
Ashim
On Wed, Nov 14, 2018 at 4:05 PM Eric Berger wrote:
> Hi Ashim,
> Per the help page for arima(), it fits an ARIMA model to the specified
> time series - but the caller has to specify the order - i.e. (p,d,q) - of
> the model.
> The
i) Your code creates w2 but references w1 to create aa.
So you needed
aa <- matrix(rep(c(0.4, 0.1, 0.2), 3), 3,3)
for a working example.
ii) This
> matrix(as.numeric(aa)[!as.numeric(aa) %in% diag(aa)],2,3)
removes any value that is present in the diagonal of aa. Look up ?"%in%" to see
what
> With your code you just remove diagonal elements from your matrix.
Worse; it removed _all_ elements from the matrix that match _anything_ in the
diagonal!
Which, in that example, was everything ...
***
This email and any
Sarah's answer is probably the best approach, but to do it using very basic R
methods that predate the very good spatial support that R now has, I would
likely do this:
## Thanks, Jim Lemon, for this step:
df1 <- read.table(text=
"latitude longitude Precip
45.5 110.5 3.2
On behalf of the American Statistical Association sections mentioned
below, I am forwarding the information that recently appeared on the
Statistical Computing Section web at
https://community.amstat.org/stat-computing/data-expo/data-expo-2019.
Old-timers will know that the ASA Data
> On 13 Nov 2018, at 16:19 , Paul Johnson wrote:
>
> Long ago, when R's t.test had var.equal=TRUE by default, I wrote some
> class notes showing that the result was equivalent to a one predictor
> regression model. Because t.test does not default to var.equal=TRUE
> these days, I'm curious
hello
for examplei have this matrix
w2<-c(0.1,0.2,0.4,0.2,0.4,0.1)aa<-matrix(w1,nrow=3,ncol=3)aa
[,1] [,2] [,3]
[1,] 0.4 0.4 0.4
[2,] 0.1 0.1 0.1
[3,] 0.2 0.2 0.2
if i use this code
matrix(as.numeric(aa)[!as.numeric(aa) %in% diag(aa)],2,3)
i will obtaine this matrix[,1] [,2] [,3]
Steve's method is very slick.
I think this is a bit easier to understand.
A <- matrix(1:9, 3, 3)
A
B <- matrix(nrow=2, ncol=3)
B[lower.tri(B, diag=TRUE)] <- A[lower.tri(A)]
B[upper.tri(B, diag=FALSE)] <- A[upper.tri(A)]
B
> A <- matrix(1:9, 3, 3)
> A
[,1] [,2] [,3]
[1,]147
[2,]
No entieno muy bien lo que preguntas. Además, los datos de Data no están
disponibles por lo que no puedo ver como quedaría el gráfico.
Sin embargo, si lo que quieres es que la linea de los mapas "tape" los
puntos que están cerca de los puntos, debes de copiar el mismo código para
visualizar el
Buenos días. Hago mis mapas con el código que os copio abajo, pero me
gustaría que me dibujase las líneas del mapa después de poner los
puntos para que se vean. Puedo hacer los puntos más pequeños, pero
entonces quedan separados y yo quiero que rellenen el mapa.
Gracias, como siempre,
Gracias Juan Diego, se trataba exactamente de eso, y claro, funcionó.
Como decía, ese es el código que utilizo para hacer los mapas, pero no
sé exactamente cómo funciona cada una de sus partes. Yo probé a copiar
esa misma línea más abajo pero me daba error porque no lo hacía
correctamente.
Estimada Dayana Muñoz Gil
Si ni comprendo mal, usted tiene en la base de datos una variable binomial,
digamos que es 0 o 1, ¿que posibilidad hay que en realidad en la base de
datos los registros sean 0, 1 y null? o 0 y 1 en el mismo registro, hay
bases de datos que permiten más de un chequeo y
22 matches
Mail list logo