search for: bouldin

...designated x as a numeric vector, so what is going on here? Thanks. > x = as.vector(c(1:12));x [1] 1 2 3 4 5 6 7 8 9 10 11 12 > mode(x) [1] "numeric" > sample(x, 3) Error in sample(x, 3) : .Random.seed is not an integer vector but of type 'list' > Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...X Y V3 3 3 1 2.236068 4 4 1 6.324555 5 5 1 5.000000 6 1 2 8.062258 10 5 2 5.656854 11 1 3 2.236068 12 2 3 9.486833 17 2 4 2.236068 18 3 4 8.062258 24 4 5 5.385165 Clearly, this is not what I intend but I cannot figure out what I've done wrong. Any help appreciated. Thanks. Jim Bouldin

...were right, I'm not sure how I would then translate those quantile values into another apply function to get the mean, since they differ from row to row. I also tried: apply(matrix,1,mean,na.rm=T,trim=.05)) and the trim argument was simply ignored Stumped. Any help appreciated. Thanks. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...1] ~ 0 + x[,3] + x[,6]))) > test2[[1]] Call: lm(formula = x[, 1] ~ 0 + x[, 3] + x[, 6]) [omitted summary output] F-statistic: 40.94 on 2 and 7 DF, p-value: 0.0001371 It does not seem to be obtainable from anova(lm(model)) either, only the p values for the individual predictors. Stumped. Jim Bouldin Research Ecologist [[alternative HTML version deleted]]

...ot;model1","model2"); I want to assign formulas, such as lm(y~x[,1]) and lm(y~x[,2]), to the variables "model1" and "model2". There are of course, many more than 2 models involved, so brute force is the option of absolute last resort. Thanks for any help. -- Jim Bouldin, Research Ecologist [[alternative HTML version deleted]]

...not figure out how to convert such strings to corresponding object names within the R workspace itself, so as to give each d.f. a distinct name. The closest I got were various attempts with the as.name function, but couldn't get that to work either. Any help appreciated. Thanks. -- Jim Bouldin, PhD Research Ecologist

...ze 3, weighted by "weights"; the mean should be 7.50) [1] 7.406977 > 7.406977/7.5 [1] 0.987597 > b = mean(replicate(1000000,(sample(x, 3))));b # (1 million samples from x, of size 3, not weighted this time; the mean should be 6.50) [1] 6.501477 > 6.501477/6.5 [1] 1.000227 Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...11 NA I need to be able to consecutively number, in new columns, the non-NA values within each column (i.e. A[1,1] A[3,2] and A[3,3] would all be set to one, and subsequent values in those columns would increase by one, until the last non-NA value is reached, if any). Any ideas? Thanks Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...> apply(a[!is.na(a)],2,sum) give me this: "Error in apply(a[!is.na(a)], 2, sum) : dim(X) must have a positive length" when > dim(a) [1] 5 6 What is the trick to calculating summary values from rows or columns containing NAs? Drives me nuts. More nuts that is. Thanks. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...one a thorough search on this topic, but from the several instructional manuals and/or tutorials I've looked at, I don't see any mention of relational database capabilities in R? Have I missed something, and if so, can someone point me in the right direction to get started? Thanks! Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...0.40 1.01 0.67 0.17 1.36 1.84 1.06 0.56 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA etc IF I solve these problems, will I in fact be able to perform individual linear regressions on the (numerous) collections of 5 to 10 rows? Thanks as always for any insight. Jim Bouldin Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...ction(x) lm(a~a2)) ...I get exactly the same result. So, which is the appropriate apply-type function when two arrays (or d.f.'s?) are involved like this? Or none of them and some other approach (other than looping which I can do but which I assume is not optimal)? Thanks for any help. -- Jim Bouldin, PhD Research Ecologist [[alternative HTML version deleted]]

...448060....854818.. "-448060" "854818" c...446828....854736.. "-446828" "854736" So far so good, but columns 1:2 will not coerce to either numeric or integer, for unknown reasons. Thanks for any help (and/or suggestions on a better way to code this). Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...pw06b" match on the first five characters; their mean would be taken whereas any columns beginning with other than "mpw06" would be excluded. I need to compare every pair of columns in the frame, and in some cases, possibly three at a time. Thanks in advance for any ideas. Jim Bouldin Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

I can't seem to find info on how to unload packages that have been loaded. My goal in doing so is to gain access to functions that have been masked out by those packages. Or is there another way to do so? Thanks in advance. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...le but I'm having trouble subsetting vectors and matrices, for example extracting all values meeting a certain criterion, from a vector. Cannot seem to figure out the correct syntax and help page not very helpful. Or should I be using some other function than subset. Thanks for any help. Jim Bouldin

...st(beta = 4)) I get this error message: Error in nls(index1 ~ beta * (1/prl3), start = list(beta = 4)) : REAL() can only be applied to a 'numeric', not a 'logical' I've got no clue as to the REAL() to which this is referring. Any help appreciated. Thanks in advance. Jim Bouldin Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...lete list of functions, use ?library(help="base")?. and not the html help screen with full package or function description like I used to. Exceedingly problematic, and I can find nothing either in the FAQs or the R search sites on what to do. Solutions much appreciated, thanks. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740

...ea ~ constant + a.est * 2.7183^(b.est * ba.beg), start = list(constant = 4, : REAL() can only be applied to a 'numeric', not a 'logical'" I can't see where there are any "logicals" in this equation to cause this problem. Any help appreciated. Thank you. Jim Bouldin

...d data frame, and then tried merging this with the original data frame and selecting only those rows container the identifier. But this did not work no matter how the arguments were altered: all records were returned instead of the uniques. Completely stumped--any help appreciated. Thanks. Jim Bouldin, PhD Research Ecologist Department of Plant Sciences, UC Davis Davis CA, 95616 530-554-1740