similar to: question about the "solve" function in library "Matrix"

Hi All. I need a bit of help extracting the residual error variance from the VarCorr structure from lmer. #Here's a 2-way random effects model lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat) #Get the structure vc.fit <- VarCorr(lmer.1) #results in..... $person 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.7755392 $rater 1 x 1 Matrix

temp.ttt <- "ttt <- 1\nttt" conn.ttt <- textConnection(temp.ttt) source(conn.ttt, echo=TRUE) ## name of variable is echoed close(conn.ttt) cat(file="c:/temp/temp.R", temp.ttt) ## name of variable not echoed source("c:/temp/temp.R", echo=TRUE) temp.abc <- "abc <- 1\nabc" conn.abc <- textConnection(temp.abc) source(conn.abc, echo=TRUE)

You may also want to check this out: plot(ttt, type = "p") points(ttt, col = ifelse(ttt < 8, "black", "red")) Eivind K. Dovik Bergen, NO On Tue, 1 May 2018, Christopher W Ryan wrote: > Excellent! Worked like a charm. Thanks. > > --Chris Ryan > > On Tue, May 1, 2018 at 4:33 PM, William Dunlap <wdunlap at tibco.com> wrote: > >> The

Excellent! Worked like a charm. Thanks. --Chris Ryan On Tue, May 1, 2018 at 4:33 PM, William Dunlap <wdunlap at tibco.com> wrote: > The ts method for plot() is quirky. You can use the default method: > > plot(as.vector(time(ttt)), as.vector(ttt), type = "p", col=ifelse(ttt<8, > "black", "red")) > > > Bill Dunlap > TIBCO Software

How would I color points conditional on their value in a plot of a time series. Something like this: ## demonstration data ttt <- ts(rpois(12, lambda = 8), start = c(2000, 1), freq = 4) ttt plot(ttt, type = "p") ## doesn't work--all points the same color plot(ttt, type = "p", col = ifelse(ttt < 8, "black", "red")) ## also doesn't work--all

The ts method for plot() is quirky. You can use the default method: plot(as.vector(time(ttt)), as.vector(ttt), type = "p", col=ifelse(ttt<8, "black", "red")) Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, May 1, 2018 at 1:17 PM, Christopher W Ryan <cryan at binghamton.edu> wrote: > How would I color points conditional on their value in a plot

Hello, everyone, I don't quite understand the following message: > TTT <- t.test(1:10, y=c(7:20)) > class(TTT) [1] "htest" > TTT@p.value Error: Trying to get slot "p.value" from an object whose class ("htest") is not defined > TTT$p.value [1] 1.855282e-05 Why the message says the class of TTT is not defined while class(TTT) gets

Hi: I have two vectors of data, x and y and I want to get the "polynomial" expansion of (x+y)^p with any integer power p in R. Suppose p=2, then I want a matrix of five vectors, namely, x y x^2 y^2 x*y. The coefficient of the polynomial is not needed. I can write it manully if p is small. But I want it in the case of p=10 or even bigger, is there any function in R can do that

Hello, I am forwarding this from ESS mailing list, as it's a failure of internal R completion system: This fails: utils:::.assignLinebuffer('iris[iris$Spec') utils:::.assignEnd(15) utils:::.guessTokenFromLine() utils:::.completeToken() utils:::.retrieveCompletions() ## -> [1] "iris[iris$Spec" This works utils:::.assignLinebuffer('iris[ iris$Spec') # note the

At times we want to convert a two-sided formula to a one-sided formula. In S we can do this by dropping the second entry in the formula. In R that object no longer has a formula class. R> ttt <- score ~ age | Infant R> class(ttt) [1] "formula" R> length(ttt) [1] 3 R> ttt[-2] [[1]] ~ [[2]] age | Infant R> class(ttt[-2]) NULL R> do.call("~",

I have a timeseries object, ts, and want to get the first date in the series into a string so I can concatenate it with a SQL query. Input and output are shown below. I must be missing something very basic, but I can't seem to pry the data ("2008-07-01") into a string variable. Any suggestions would be appreciated. Thank you, Andrew =====script: class(ts) (ts[1,0]) #returns

We have frequent timeout problems after upgrading our imap servers from dovecot 1.1.16 to dovecot 1.2.11. One server acts as proxy only, and the other one is the "real" imap server". The credentials for the proxy service are stored in a remote MYSQL database. There were no trouble with dovecot 1.1.16. But now, with the most recent version, we get frequent login failures. It

I am trying to copy some digital images to a directory on my Redhat 7.3 Samba 2.2.3a-6 server and am getting extra files in the directory. I have included my SMB.Conf file contents below in case that is relevant. The files being copied reside on a Windows 2000 workstation and were created by Jasc Image robot. These files would have been accessed by both a Windows 2000 machine as well as a Windows

Hi, I encountered a problem of not being able to read in an excel spreadsheet using read.xls() in the xlsReadWrite package. can anyone help? Here is an example code write.xls(matrix(1:9,nrow=3),"ttt.xls") read.xls("ttt.xls") Error in read.xls("ttt.xls") : Unexpected error. Message: Can't find the file "ttt.xls" The "ttt.xls" file was

Hello, we are using the vfs module recycle with the following config: vfs objects = recycle recycle: repository = .Papierkorb recycle:directory_mode = 0777 recycle:subdir_mode = 0777 recycle: keeptree = Yes recycle: exclude = *.tmp, *.temp, *.log, *.ldb recycle: exclude_dir = tmp recycle:versions = Yes Problem is that everybody can see deleted documents in the recycle bin. But if I set

Why fails nls with "singular gradient" here? I post a minimal example on the bottom and would be very happy if someone could help me. Kind regards, ########### # define some constants smallc <- 0.0001 t <- seq(0,1,0.001) t0 <- 0.5 tau1 <- 0.02 # generate yy(t) yy <- 1/2 * ( 1- tanh((t - t0)/smallc) * exp(-t / tau1) ) + rnorm(length(t))*0.01 # show the curve

Dear all, I'd like to subset a df within a function, and use select for choosing the variable. Something like (simplified example): mydf <- data.frame(a= 0:9, b= 10:19) ttt <- function(vv) { tmpdf <- subset(mydf, select= vv) mean(tmpdf$vv) } ttt(mydf$b) But this is not the correct way. Any help? Thanks in advance Juli

I apparently have a dead FXS module. Is there any kind of test I can run on it (on a live system) to determine if its good or bad? The telephone has no dialtone, gets no calls. It has been working for several months, and just quit yesterday. Thanks! Found a Wildcard TDM: Wildcard TDM400P REV I (4 modules) Freshmaker version: 73 Freshmaker passed register test Module 0: Installed -- AUTO

Dear all, I have looked for an answer for a couple of days, but can't come with any solution. I have a set of functions, say: > t0 <- function(x) {1} > t1 <- function(x) {x} > t2 <- function(x) {x^2} > t3 <- function(x) {x^3} I would like to find a way to add up the previous 4 functions and obtain a new function: > rrr <- function(x) {1+x+x^2+x^3} without,

I have just installed a samba 4.7.4 from source. All right! wbinfo -u  and wbinfo -g return the list of users and groups from the domain But I can not use they in the DC filesystem. For instance: # touch ttt # chgrp 'domain admins' ttt chgrp: invalid group: ‘domain admins’ That should I do to ?