similar to: persp() problem

Hi all, I am trying to create a function that will open connections to all files of one type within the working directory. I've got the function to open the connections, but I am having a bugger of a time trying to get these connections named as objects in the workspace. I am at the point where I can do it outside of the function, but not inside, using assign. I'm sure I'm

Hi all I am using the caret package and having difficulty in obtaining the results using regression, I used the glmnet to model and trying to get the coefficients and the model parameters I am trying to use the extractPrediction to obtain a confusion matrix and it seems to be giving me errors. x<-read.csv("x.csv", header=TRUE); y<-read.csv("y.csv", header=TRUE);

thanks for your attention. temporarily with the following script I'm looking at myself. #!/bin/bash INTERVAL=180 LOG=true function log { if $LOG; then echo $1 fi } while true; do testx=`ps -auxw | grep "smbd" | wc -l` testy=`netstat -ap | grep samba | grep sock | wc -l` if test "$testx" -gt "50" -a "$testy" -gt "50" ;then echo "samba

Readers, How to solve this unicode input error please? > postscript("~/tmp/test.eps",width=100/25.4,height=100/25.4,horizontal=FALSE,onefile=TRUE,paper="special") > testx<-seq(1:5) > testy<-seq(1:5) > plot(testy~testx) > mtext('text (O?)\n more text',side=3,line=1) Warning messages: 1: In mtext("text (O?)\n more text", side = 3, line =

I reported the problem two months ago . I knew this would happen . :(( 2016-09-17 22:42 GMT+03:00 Edson Tadeu Almeida da Silveira via samba < samba at lists.samba.org>: > Sorry, the problem came back, after 3 hours. > > att, > > 2016-09-17 13:11 GMT-03:00 Edson Tadeu Almeida da Silveira < > edson.tadeu at gmail.com>: > > > I don't know if it is a

Hi R users, I have a couple of questions about some problems that I am facing with regard to numerical integration and optimization of likelihood functions. Let me provide a little background information: I am trying to do maximum likelihood estimation of an econometric model that I have developed recently. I estimate the parameters of the model using the monthly US unemployment rate series

Dear R-helper, How could I count only some variable was exist after running sample (random) function. For example, > testx <- factor(c("Game","Paper","Internet","Time","Money")) > for(i in 1:2) { + x <- sample(testx,replace=TRUE) + print(x) + } [1] Money Money Time Internet Time Levels: Game Internet

Hello Please see the example below > class(testX) [1] "matrix" > class(testX[1,]) [1] "numeric" Why not matrix? What am I missing here? Is there a way to keep the same class? The reason for the question is that I want to implement a k-step ahead prediction for my own routines and R wrecks does not seem to like [1,] as shown below. >

Hi, I have a question about using lm on matrix, have to admit it is very trivial but I just couldn't find the answer after searched the mailing list and other online tutorial. It would be great if you could help. I have a matrix "trainx" of 492(rows) by 220(columns) that is my x, and trainy is 492 by 1. Also, I have the newdata testx which is 240 (rows) by 220 (columns). Here is

All, For my understanding, I wanted to see if I can get glmnet predictions using both the predict function and also by multiplying coefficients by the variable matrix. This is not worked out. Could anyone suggest where I am going wrong? I understand that I may not have the mean/intercept correct, but the scaling is also off, which suggests a bigger mistake. Thanks for your help. Juliet Hannah

I have gone over the examples and can't figure this out: rho<-.77 text(x=.05,y=.5,paste(expression(rho),rho)) I was hoping to get this to print a Greek rho with 0.77 beside it. Instead I get: rho 0.77 (i.e. Roman lettering) The help on expression() is quite opaque so I don't understand how it works. Thanks for any help. Bill Simpson

I am looking for examples of code that demonstrates the fine tuning of the strip panels in lattice graphics and uses plotmath characters. The code for the graphic is as follows: xyplot(lagy ~ n | rho1 * rho2, data= data, layout=c(2,6), span = 1, xlab = "Sample Size", ylab = "Bias in the Coefficient for the Lag of X", type = "o") rho1 is a four level factor

Using str() in a function. I am in the early phase of learning R, and I find I spend a lot of time trying to figure out what is actually in objects I have created or read in from a file. I'm trying to make a simple little function to display a couple of things about a object, let's say the summary() and the str(), sequentially, preferably without a bunch of surplus lines between them. I

Hi, I'm trying to write a simple function like case1 <- function (m, cov, Q, R) { theta <- (acos(R/sqrt(Q^3))) beta <- (-2)*sqrt(Q)*cos(theta/3)+m[1]/3 rho1 <- (-2)*sqrt(Q)*cos((theta+2*pi)/3)+m[1]/3 rho2 <- (-2)*sqrt(Q)*cos((theta-2*pi)/3)+m[1]/3 stderrb <- deltamethod( ~(-2)*sqrt(Q)*cos(theta/3)+x1/3,m,cov) stderrr1 <- deltamethod(

[Apologies in advance if this is too "statistics" and not enough "R".] I've got an experiment with two sets of treatments. Each subject either received all treatments from set A or all treatments from set B. I can compute the N pairwise correlations for all treatments in either set using cor(). If I take the mean of these N pairwise correlations, I see that the effects

Is it possible to run MS Word on a LInux box without Windows installed? If so, would someone be willing to help me get it running? I would like to use MS Word ver 7, but want nothing to do with Windoze. Word Perfect 8 for Linux can get testy at times. cedric

Dear all, I am confused with the behaviour of survfit with newdata option. I am using the latest version R-2-15-0. In the simple example below I am building a coxph model on 90 patients and trying to predict 10 patients. Unfortunately the survival curve at the end is for 90 patients. Could somebody please from the survival package confirm that this behaviour is as expected or not - because I

Dear folks: Let?s suppose I want a function to print return the name of the object passed to it. > myname <- function(object) {out<-deparse(substitute(object)); out} This works fine on a single object: > O1 <-c(1:4) > myname(O1) [1] "O1" However it does not work if you use lapply to pass it the same object from a list: > O2 <-c(1:4) > object.list <-

I've found a problem when using categorical variables in lmp() from package lmPerm According to help(lmp): "This function will behave identically to lm() if the following parameters are set: perm="", seq=TRUE, center=FALSE.") But not in the case of including categorical variables: require(lmPerm) set.seed(42) testx1 <- rnorm(100,10,5) testx2 <-

hey, my function doesn?t work. can somebody help me? the graphic doesn?t work and also the function. thnx a lot. N=10 n=100 p_0=c(1/5,1-1/5) power = function(p,m) { set.seed(1000) H=matrix(0,nrow=N,ncol=1) for(i in 1:N) { x <- matrix(rnorm(n, 0, 0.5), ncol = m) y <- matrix(rnorm(n, 0, 0.8), ncol = m) l <- diag(cor(x, y)) q_1 = qnorm(0.05, 0, 0.05) q_2 = qnorm(1 - 0.05, 0, 0.05)