[R] grep

[Steven Yen]

Now I've found another way to make it work. All I need is to pick up the 
names in the column (x.1.age...).

> v<-pr(goprobit.p); v

Maximum-Likelihood Estimates weighted = FALSE iterations = 5 logLik = 
-14160.75 finalHessian = TRUE Covariance matrix is Robust Number of 
parameters = 66 Sample size = 17922 est se t p g sig x.1.age 0.0341 
0.0138 2.4766 0.0133 -3.8835e-04 ** x.1.sleep -0.1108 0.0059 -18.6277 
0.0000 -4.4572e-04 *** x.1.primary -0.0694 0.0289 -2.4002 0.0164 
-9.9638e-06 ** x.1.middle -0.2909 0.0356 -8.1657 0.0000 -1.4913e-05 *** 
... > rownames(v) [1] "x.1.age" "x.1.sleep" "x.1.primary" "x.1.middle" 
[5] "x.1.high" "x.1.somewhath" "x.1.veryh" "x.1.somewhatm" ...Treating 
it as the first column does not work, because the first column contains 
the numbers. > v[,1] [1] 0.0341 -0.1108 -0.0694 -0.2909 -0.4267 -0.6188 
-0.7580 -0.3413 -0.3813 [10] -0.3101 -0.2977 -0.6310 -0.9132 0.1885 
-0.0887 -0.0850 0.0847 -0.1588

On 7/12/2024 11:57 PM, Steven Yen wrote:
> Thanks. First and second of the following worked, but the third (with 
> coef) did not. This may be because I programmed my own estimation 
> program. In short,
>
> names(goprobit.p$est) was recognized, but
> names(coef(goprobit)) was not
>
> > which(grepl("very|somewhat",names(goprobit.p$est)))
>  [1]  6  7  8  9 10 11 12 13 28 29 30 31 32 33 34 35 50 51 52 53 54 55 
> 56 57
> > jj<-grep("very|somewhat",names(goprobit.p$est)); length(jj)
> [1] 24
> > jj<-grep("very|somewhat",names(coef(goprobit))); length(jj)
> Error: object 'goprobit' not found
>
> On 7/12/2024 11:23 PM, Rui Barradas wrote:
>>
>> Hello,l
>>
>> Though the question is already answered, here is another answer to 
>> what is 'x'.
>> The output in the OP is not a lm or glm output but if your regression 
>> model was programmed according to recommended practices, there must 
>> be a 'coefficients' member in the list or object it returns and the 
>> following should work.
>>
>>
>> # this is 'x', a named character vector
>> coef(fit)
>> #
>> fit |> coef() |> names() |> grep("somewhat|very", x = _)
>>
>>
>> Hope this helps,
>>
>> Rui Barradas
>>
>> Às 10:26 de 12/07/2024, Steven Yen escreveu:
>>> Thanks. In this case below, what is "x"? I tried rownames(out) which 
>>> did not work.
>>>
>>> Sorry. Does this sound like homework to you?
>>>
>>> On 7/12/2024 5:09 PM, Uwe Ligges wrote:
>>>>
>>>>
>>>> On 12.07.2024 10:54, Steven Yen wrote:
>>>>> Below is part a regression printout. How can I use "grep" to identify
>>>>> rows headed by variables (first column) with a certain label. In this
>>>>> case, I like to find variables containing "somewhath",
>>>>> "veryh", "somewhatm", "verym", "somewhatc", "veryc","somewhatl",
>>>>> "veryl". The result should be an index 6:13 or 
>>>>> 6,7,8,9,10,11,12,13. Note
>>>>> that they all contain "somewhat" and "very". Thanks.
>>>>
>>>> Sounds like homework?
>>>>
>>>> which(grep("very|somewhat", x))
>>>>
>>>> Best,
>>>> Uwe Ligges
>>>>
>>>>
>>>>> est se t p g sig x.1.age 0.0341 0.0138 2.4766 0.0133 -3.8835e-04 **
>>>>> x.1.sleep -0.1108 0.0059 -18.6277 0.0000 -4.4572e-04 *** x.1.primary
>>>>> -0.0694 0.0289 -2.4002 0.0164 -9.9638e-06 ** x.1.middle -0.2909 
>>>>> 0.0356
>>>>> -8.1657 0.0000 -1.4913e-05 *** x.1.high -0.4267 0.0463 -9.2118 0.0000
>>>>> -3.6246e-05 *** x.1.somewhath -0.6188 0.0256 -24.1971 0.0000 
>>>>> -3.1337e-05
>>>>> *** x.1.veryh -0.7580 0.0331 -22.8695 0.0000 -2.9558e-05 ***
>>>>> x.1.somewhatm -0.3413 0.0426 -8.0112 0.0000 -1.8920e-05 *** x.1.verym
>>>>> -0.3813 0.0446 -8.5413 0.0000 -4.4029e-05 *** x.1.somewhatc -0.3101
>>>>> 0.0649 -4.7783 0.0000 -1.4353e-05 *** x.1.veryc -0.2977 0.0648 
>>>>> -4.5910
>>>>> 0.0000 -4.8986e-05 *** x.1.somewhatl -0.6310 0.0424 -14.8846 0.0000
>>>>> -1.9543e-05 *** x.1.veryl -0.9132 0.0462 -19.7525 0.0000 
>>>>> -4.4603e-05 ***
>>>>> ...
>>>>>
>>>>>     [[alternative HTML version deleted]]
>>>>>
>>>>> ______________________________________________
>>>>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>>>> PLEASE do read the posting guide 
>>>>> http://www.R-project.org/posting-guide.html
>>>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>> ______________________________________________
>>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide 
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
	[[alternative HTML version deleted]]