[R] Non-conformable arrays - Error in var(betas)

Hello, If there is anyone who can offer me some guidance on how to resolve the error below, I would be very grateful. I am attempting to conduct a 'Resource Utilization Function' analysis using the 'ruf' package. My code is below: library(ruf) d78423 <- read.csv("RUF/78423_RelUse.csv") #

[R] predict() works with the design matrix but throws error with some rows of that matrix

Hello, I am having problems with predict() after a multinomial logit regression by multinom(). I generate a design matrix with model.matrix() and use it to estimate the model. Then, if I pass the entire design matrix to predict(), it returns the same output as fitted(), which is expected. But if I

Re: [R] is.na behavior

I don't think the question is about what an expression evaluates to, but what it is. I don't think it makes sense for an expression to be NA, but I am not a language designer. The expression NA tells the compiler that a literal NA_logical_ value is to be returned from evaluating that expressi

Re: [R] Change breaks x axis in plot

xx <- seq(1, 2, by = 0.1) yy <- sample(1:20, 11) plot(xx, yy, xaxt="n") axis(1, at=1: 2, labels=letters[1:2]) may be what you want. John Kane Kingston ON Canada > -Original Message- > From: luysgar...@gmail.com > Sent: Thu, 19 Nov 2015 16:18:16 -0300 > To: r-help@r-project.org > Subjec

Re: [R] Change breaks x axis in plot

Hello, Maybe something like this? plot(1:2, xaxp = c(1,2,1)) See ?par. Hope this helps, Rui Barradas   Citando Luis Fernando García : > I have the following graph: > my problem is the scale, I just need it to show the values for 1 and 2 and > remove the intermediate values (1.2,1.4,1.6,1.8).

Re: [R] is.na behavior

I am not sure I undestand the issue. But if the question is to decidedif an expression evaluates to NA, using eval should solve the problem. In fact, I do not really understand what an NA expression, and not an expression evaluating to NA, means. signature.asc Description: Message signed with O

[R] Change breaks x axis in plot

I have the following graph: my problem is the scale, I just need it to show the values for 1 and 2 and remove the intermediate values (1.2,1.4,1.6,1.8). How could I do it? I'm, attaching the script code for the plot. Thanks in advance ###AV=mean values ##SD=Sd values plot(D ~ ferre$Cantidad, t

Re: [R] [FORGED] How to find the likelihood, MLE and plot it?

> On 19 Nov 2015, at 18:12 , David Winsemius wrote: > >> >> On Nov 19, 2015, at 7:35 AM, C W wrote: >> >> ah. Let me fix that and get back to you. >> >> On a side note, why can't I put fn inside curve(), why do I have to use >> Vectorize()? >> >> Here's the code and the message: > >>

Re: [R] [FORGED] How to find the likelihood, MLE and plot it?

> On Nov 19, 2015, at 7:35 AM, C W wrote: > > ah. Let me fix that and get back to you. > > On a side note, why can't I put fn inside curve(), why do I have to use > Vectorize()? > > Here's the code and the message: > fn <- function(theta){ > > sum(0.5 * (xvec - rep(theta, 7)) ^

Re: [R] is.na behavior

Richard, I think the reason that this gives the warning is for the rest of us who don't think about asking about missing values in non-data objects. I could imagine someone choosing a poor name for a variable and doing something like: mean <- mean(x) is.na(mean) which would then tell them wheth

Re: [R] [FORGED] How to find the likelihood, MLE and plot it?

ah. Let me fix that and get back to you. On a side note, why can't I put fn inside curve(), why do I have to use Vectorize()? Here's the code and the message: > curve(fn, -5, 20) Error in curve(fn, -5, 20) : 'expr' did not evaluate to an object of length 'n' On Thu, Nov 19, 2015 at 10:29 AM

Re: [R] JAGS 4.x, rjags 4.x problems | Linux

Hi R-users, talking about the compilation of rjags (gnu 4.4.7 - blas and lapack routines) package in R I have solved the issue in the thread as follows: Before compiling JAGS I need to recompile BLAS and LAPACK libraries with -fPIC instruction: BLAS Makefile: [root@n14 library]# cat /home/

[R] knn - random result although use.all=TRUE

Dear all, I have this toy example to work with k-nn classification approach. (My data, code and results are at the end of the message) Working with knn function in library class and setting the parameter use.all=TRUE, I would not expect a random answer. Nevertheless I get a different answer ea

Re: [R] [FORGED] How to find the likelihood, MLE and plot it?

On 19 Nov 2015, at 16:17 , C W wrote: > Hi Rolf, > > I think the MLE should be 1.71, no? And yes, I am aware of the > maximum=TRUE argument. I still feel something is wrong here. > Just read more carefully what Rolf said: Your fn is MINUS the log-likelihood. So the graph is upside-down. -

Re: [R] [FORGED] How to find the likelihood, MLE and plot it?

Hi Rolf, I think the MLE should be 1.71, no? And yes, I am aware of the maximum=TRUE argument. I still feel something is wrong here. Thanks! On Wed, Nov 18, 2015 at 6:23 PM, Rolf Turner wrote: > On 19/11/15 11:31, C W wrote: > >> Dear R list, >> >> I am trying to find the MLE of the likeliho

[R] improve my ggplot look

Dear all,the following line of code print me a map of an area with the points I need. I only new two minor adjustments     ggmap(mp, darken = 0) + geom_point(aes(Longitude, Latitude, colour =Error), data = PlotPoints, size = 6)+ scale_colour_gradient2(low=muted("red"),mid="green", high=muted("bl

[R] Prediction data gaps by "neuralnet" package

Dear R users; I am trying to predict the gaps in my data set (colorofil.A) by using "neuralnet" package. But when I run the following codes it gives me error. Could anyone help me to fix the problem? Bests. > library("neuralnet") > head(ss1,2)[,c(5,24)] Season colorofil.A 1 Sp