similar to: setting a number of values to NA over a data.frame.

A stupid question but I just cannot see how to do this. I have a loop that does some calculations and puts the results in a vector for each iteration, but I cannot see how to get the data out of the loop in such a way that I can use it. I can print it but how do I get it into a set of vectors or what ever. Any help gratefully received. Thanks Example cata <- c( 3,5,6,8,0, NA) catb <-

I am stuck on a simple problem where an example works fine but the real one does not. I have a data.frame where I wish to sum up some values across the rows and create a new data.frame with some of old data.frame variables and the new summed variable. It works fine in my simple example but I am doing something wrong in the real world. In the real world I am loading a labeled data.frame. The

I have a simple data base and I want to produce tables for each variable. I wrote a simple function fn1 <- function(x) {table(x)} where x is a matrix or data.frame. and used apply to produce a list of tables. Example below. How do I apply the colnames from the matrix or names from the data.frame to label the tables in the results in the list. I know that I can do this individually but

Dear R'sians, Would really appreciate if you could suggest a more efficient way to order the columns of a dataset. The column names of the dataset contain indices separated by a period. Following are examples of my code and the dataset. oC <- function(tg=x2) { lth <- length(grep("T",names(tg))) thix <-

hi I plot a series of observation data every minutes in a day as the attachment below plot(wnd,type='l',lty=1,col='red',lwd=1,xlab=xxlab,ylab=yylab,ylim=YY) In the figure, the x-axis tickment is the number of data How can I change it fore example 1h 2h 3h 4h and so on ? -- TANG Jie Email: totangjie at gmail.com Tel: 0086-2154896104 Shanghai Typhoon Institute,China

I want to calculate the "scat" etc, values for each level of id, output a vector of these values that is a <- c(scat, sdog, srat, sbat ) and do an rbind on them. Each level of id has a different value of rate. So far it it looks to me like I can do this by a series of repeated subsets of Df butI would think that I should be able to do this with a function but I am not having

Hello R-Users, I have the following problem with the subset function: See the following simple linear model. Here everything works fine: >germany<-lm(RENT~AGE1, in.mi01) However, if a use the same regression equation and only specify a subset, I get an error message: > berlin<-lm(RENT~AGE1, in.mi01, subset=C_X01=="Berlin") Error in lm.fit(x, y, offset

Hi Dennis, SOLVED!!! My thanks to both you, John, and others who chimed in. Took a little more digging before finally working, but it is working! Here's a little more of what I did and the ultimate resolutions: I installed the reshape2 library from CRAN. Executed > (.packages()) [1] "reshape2" "stats" "graphics" "grDevices"

Hi all, I'm having some trouble with rbind - this may be a bug or it may be my misunderstanding. If I do fileName <- paste(tempdir(),"test.txt",sep="/") file.create(fileName) x <- read.table(fileName, col.names=c("one","two","three")) I get a data frame with no rows, as documented. If I then try to rbind this with another data frame

#I have a data set that looks like this. A bit more complicated actually with # three factor levels but these calculations need to be done on one factor at a #I then have a set of different rates that are applied #to it. #dataset cata <- c( 1,1,6,1,1,2) catb <- c( 1,2,3,4,5,6) doga <- c(3,5,3,6,4, 0) data1 <- data.frame(cata, catb, doga) rm(cata,catb,doga) data1 # start rates #

## you may need to install HH install.packagess("HH") library(HH) hellisheidi <- read.table(text=" Component Sample1 Sample2 Sample3 CaO 45 52 48 SiO2 25 22 18 Al2O3 15 11 14 TiO2 6 5

Hello I have the following data frame which I read from an EXCEL file, and when i try to call one of its columns with a space in their names I am not being able to. For example if I do EURODOLLAR$ED1.Comdty Date I obtain the following error: Error: inesperado símbolo en "EURODOLLAR$ED1.Comdty Date" I have also tried using . or _ instead of the space and have obtained no succes. How do I

Dear Richard,? Thank you very much for your reply. I went through the code and it worked. I was also able to change the colours.? I was wondering if I can change the legend position; instead of being in the bottom to be on the left side.? I tried the following but without any success strip = FALSE strip.right = TRUE likert(t(hellisheidi), ReferenceZero=.5, xlab="X-lab",

Dear all,? I am trying to plot the following table in stacked barplot in percentages and also horizontal. Component? ? ? Sample 1? ? ? Sample 2? ? ? Sample 3CaO????????????????? ? 45????????????? ? 52????????????? ? 48SiO2????????????????? ?25????????????? ? 22????????????? ? 18Al2O3?????????????????15????????????? ? 11????????????? ? 14TiO2? ? ? ? ? ? ? ? ? ? ?6? ? ? ? ? ? ? ? ? 5? ? ? ? ? ? ? ?

Dear all here is my data called "rglp" structure(list(vzorek = structure(1:17, .Label = c("179/1/1", "179/2/1", "180/1", "181/1", "182/1", "183/1", "184/1", "185/1", "186/1", "187/1", "188/1", "189/1", "190/1", "191/1", "192/1",

R Users: I am trying to estimate a model of fertility behaviour using birth history data with maximum likelihood. My code works but is extremely slow (because of several for loops and my programming inefficiencies); when I use the genetic algorithm to optimize the likelihood function, it takes several days to complete (on a machine with Intel Core 2 processor [2.66GHz] and 2.99 GB RAM). Computing

Hello, I was trying to build some functions which I would like to integrate over an interval using the function 'integrate' from the 'stats' package. As an example, please consider the function h(u)=sin(pi*u) + sqrt(2)*sin(pi*2*u) + sqrt(3)*sin(pi*3*u) + 2*sin(pi*4*u) Two alternative ways to 'build' this function are as in f and g below: coeff<-sqrt(1:4)

https://bugs.freedesktop.org/show_bug.cgi?id=92306 Bug ID: 92306 Summary: GL Excess demo renders incorrectly on nv43 Product: Mesa Version: git Hardware: x86 (IA32) OS: Linux (All) Status: NEW Severity: normal Priority: medium Component: Drivers/DRI/nouveau Assignee: nouveau

As part of carrying out a complicated maximum likelihood estimation, I am trying to learn to program likelihoods in R. I started with a simple probit model but am unable to get the code to work. Any help or suggestions are most welcome. I give my code below: ************************************ mlogl <- function(mu, y, X) { n <- nrow(X) zeta <- X%*%mu llik <- 0 for (i in 1:n) { if

Hello, I am looking into description of surf.ls(spatial) and see under value $beta - the coefficients. When I use polynomial of degree 2 to fit surface I expect to get 4 coefficients: z = a_1 x^2 + a_2 xy + a_3 y^2 + a_4 What do beta really stand for and why do I get $beta vector of length 6? Thakns, Mark